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Drag  reduction  by  ejection  of  high  molecular  weight  polymers  on  free- 
running  bodies  of  revolution  has  been  demonstrated  repeatedly.  The  quanti- 
ties ol  polymer  required  have  made  the  gains  achieved  marginal  from  a volume 
utilization  tradeoff.  The  ejection  process  is  hypothesized  to  be  the  control 
J ing  factor.  Limited  data  obtained  in  pipe  flow  and  flat  plate  flow 
experiments  on  ejection  into  developing  boundary  layers  indicate  a drastic 
reduction  in  polymer  requirements  for  equivalent  percent  drag  reductions.  _ 
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Extension  to  axi  syinmet  r ic  I low  could  result  in  significant  achievable  gains 
in  volume  utilization. 

This  study  examines  these  polymer  ejection  processes  through  measure- 
ment of  wall  and  boundary  layer  concentration  prol iles  and  through  a photo- 
graphic study  of  Lhe  boundary  layer.  Tests  performed  with  fresh  water  ejec- 
tion and  solutions  of  the  drag  reducing  polymer.  Polyox  WSK-101,  lead  to  the 
hypothesis  that  optimal  ejection  lor  minimum  polymer  usage  requires  ejection 
into  a laminar  boundary  layer  prior  to  turbulent  flow  transition. 

Analytical  routines  are  developed  which  predict  boundary  layer  parameter,; 
and  polymer  wall  concentrations  for  this  postulated  optimal  ejection  process 
or  the  subopt imal  case.  Limited  veril icaLion  ol  the  model  is  made. 
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Drag  reduction  by  ejection  o(  high  molecular  weight  polymers 
on  free-running  bodies  of  revolution  has  been  demonstrated  repeatedly. 
The  quantities  of  polymer  required  have  made  the  gains  achieved 
marginal  from  a volume  utilization  tradeoff.  The  ejection  process 
is  hypothesized  to  be  the  controlling  factor.  Limited  data 
obtained  in  pipe  flow  and  flat  plate  flow  experiments  on  ejection 
into  developing  boundary  layers  indicate  a drastic  reduction  in 
polymer  requirements  for  equivalent  percent  drag  reductions. 

Extension  to  axisymmetric  flow  could  result  in  significant  achievable 
gains  in  volume  utilization. 

The  object  of  the  research  described  uerein  was  to  examine 
the  processes  described  above  through  measurement  of  wall  and  boun- 
dary layer  concentration  profiles  and  through  a photographic  study 
of  the  boundary  layer.  Tests  were  performed  with  fresh  water 
ejection  and  solutions  of  the  drag  reducing  polymer,  Polyox  WSR-301. 
Predictive  analytical  routines  were  developed  and  experimentally 
verified . 

The  experimental  apparatus  used  in  this  experiment  was  a drop 
tank  with  several  velocity  measurement  stations  and  photographic 
equipment.  Several  specially  designed  axisymmetric  bodies  were 
constructed.  Three  of  these  ejected  dye  by  aspiration  at  the 
minimum  pressure  point  near  the  nose.  The  fourth  had  a capability 
to  eject  through  a nose-screen  type  orifice  at  a constant  rate 
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and  withdraw  wail  or  boundary  layer  samples  at  four  stations  along 
the  body.  All  bodies  had  a maximum  diameter  of  three  inches  and 
applied  a half-body  forebodv.  Trie  three  dye  ejecting  bodies  incor- 
porated different  tail  configurations:  6°  cone  tail,  L/D  = 8.33; 

12°  cone  tail,  I./D  = 6.7;  and  spherical  tail,  L/D  = 5. 54.  The 
polymer  ejection  body  was  similar  In  configuration  to  the  6°  tail 
body . 

Tests  were  performed  in  a Reynolds  number  range  of  one  to  five 
million.  The  first  series  of  tests  were  performed  in  fresh  water 
and  in  a polymer  "ocean"  of  various  concentrations;  1.25,  2.5, 

5,  10,  20,  50,  and  60  WPPM.  Skin  friction  reduction  obtained  on 
all  bodies  approached  70  percent  at  polymer  concentrations  of 
20  WPPM,  agreeing  well  with  data  from  other  experimenters.  Total 
drag  reductions  of  33  percent  for  the  6°  body,  16  percent  for  the 
12°  body,  and  10  percent  for  the  spherical  tail  body  were  obtained, 
the  effect  of  higher  percentage  form  drag  being  evident. 

Photographic  studies  of  dyed  boundary  layers  in  these  tests 
and  in  ejection  tests  with  50  WPPM,  500  WPPM,  and  1000  WPPM 
displayed  several  interesting  characteristics.  Addition  of  small 
quantities  of  polymer,  2.5  WPPM,  eliminated  the  fine-scale  turbulent 
structure  leaving  only  coarse  turbulence.  Higher  concentrations 
resulted  in  suppression  of  the  coarse  turbulence  and  extreme 
thinning  of  the  boundary  layer.  In  the  ejection  tests,  concentra- 
tions of  500  WPPM  and  1000  WPPM  displayed  no  turbulent  structure  or 


mixing.  A dye  streaking  phenomenon  with  spacing  equivalent  to  that 
characteristic  oi  laminar  sublayer  streaks  was  displayed.  The 
number  of  streaks  was  also  approximately  equal  to  the  number  of 
ejection  holes,  approximately  700,  supporting  the  fact  that  mixing 
had  essentially  ceased. 

Concentration  profiles  were  measured  in  the  ejection  tests 
using  a fluorometric  method.  Tracer-contaminated  water,  5 WPPM, 

10  WPPM,  20  WPPM,  50  WPPM,  500  WPPM,  and  1000  WPPM  solutions, 
were  ejected.  Four  measurement  stations,  the  last  at  an  x/L 
of  .48,  were  sampled.  The  water  data  displayed  expected  concen- 
tration profiles,  agreeing  well  with  those  of  other  investigators. 

At  all  polymer  concentrations  tested,  unexpected  results  were 
achieved.  Wall  concentrations  remained  at  levels  predicted  by 
molecular  diffusion.  As  evidenced  from  the  photographs,  no  diffu- 
sion was  occurring.  An  initial  mixing  zone  neglected  by  most  inves- 
tigators was  the  controlling  zone  in  these  tests.  Only  for  the  5 
WPPM  case  did  a higher  diffusion  rate  begin  part  way  along  Lhe 
body.  The  data  for  this  axi symmetric  case  we  re  compared  with  data 
from  the  single  flat  plate  experiment  evidencing  the  same  phenomenon. 
A dimensionless  distance  developed  in  the  flat  plate  experiment  used 
to  predict  the  extent  of  this  initial  zone  was  modified  and  compared 
favorably  with  this  experiment.  It  is  hypothesized  from  these 
tests  that  optimal  ejection  for  minimum  polymer  usage  requires 
ejection  into  a laminar  boundary  layer  prior  to  transition  to 
turbulent  flow. 
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Analytical  routines  are  developed  which  predict  boundary 
layer  parameters  and  polymer  wall  eoncentrat Jons  for  this  postulated 
optimal  ejection  process  or  the  suboptimal  case.  Limited  verifica- 
tion of  the  model  is  made. 
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INTRODUCTION 


A.  General  Considerations 

High  speed  flow  of  both  gases  and  liquids  past  a fixed 
boundary  generates  a turbulent  boundary  layer  which  exerts  both 
a fluctuating  normal  stress  (flow  noise)  and  fluctuating  shear 
stress  (whose  average  is  drag)  on  the  body.  The  reduction  of 
resistance  to  the  turbulent  shear  flows  of  liquids  through  the 
addition  of  small  quantities  of  high  molecular  weight  polymers 
into  the  boundary  layer  of  the  flow  has  been  demonstrated 
repeatedly.  This  phenomenon  of  drag  reduction  has  far-reaching 
importance  in  the  reduction  of  power  required  for  pumping  of 
fluids  or  transport  of  bodies  through  liquids. 

Until  recently,  the  phenomenon  of  drag  reduction  was  con- 
fined essentially  to  characterizing  the  manner  in  which  the 
boundary  layer  velocity  profiles  were  affected,  thus  allowing 
predictions  of  boundary  layer  thickness  and  shear  drag.  Addi- 
tionally, the  early  work  primarily  dealt  with  well-developed  pipe 
flows  where  the  observation  of  drag  reduction  was  first  made. 

An  added  complication  results  from  the  fact  that  data  obtained 
with  commonly  used  instrumentation  have  been  found  to  be  affected 
by  the  polymers  resulting  in  erroneous  readings.  As  a result, 
the  mechanism  of  shear  stress  reduction  remains  undefined  although 


emperical  techniques  predicting  the  benefits  of  application  of  the 
phenomenon  have  been  developed. 


f ' '5g-!^^  jf  WfFWf  .•  f,  ,'.-f sps^TTV*  >w»v»  * 

- W *J~  - «1*»  . . — . .^-—  i_ 


Interest  In  the  problem  Imh  expanded  with  recent  Invest  Igntors 
examining  the  effect  of  polymers  on  external  flows,  as  over  flat 
plates.  Similar  results  as  in  the  pipe  flows  have  been  achieved. 
The  magnitude  of  the  turbulent  fluctuations  have  been  found  to 
decrease  with  a thickening  of  the  laminar  sublayer  and  a lowering 
of  the  shear  stress. 

Shear  stress  reduction  for  external  flows  over  bodies  can 
be  achieved  by  proper  ejection  of  polymer  solutions.  A limit  to 
the  amount  of  drag  reduction  achievable  with  polymer  use  has  been 
demonstrated  and  increases  in  shear  stress  has  been  obtained  when 
excess  polymer  is  used. 

It  is  evident  that  a better  understanding  of  the  mechanisms 
of  drag  reduction  in  developing  flows  with  polymer  ejection  is 
required  if  an  advantage  is  to  be  taken  of  the  phenomenon.  For 
most  cases,  of  interest  in  external  flows,  a carried  polymer  supply 
would  be  used  for  ejection  into  the  boundary  layer.  Limits  on 
diffusion  result  in  a far  too  high  concentration  initially  - with, 
possibly,  shear  drag  increases  - and  then,  further  downstream  on 
the  body,  too  little  polymer  resulting  in  insignificant  gains  in 
drag  reduction.  One  might  envision,  to  circumvent  the  problem, 
multiple  ejection  ports  to  optimize  the  polymer  wall  concentra- 
tion and  maximize  the  drag  reduction. 

Many  investigators  have  addressed  portions  of  the  problem 
outlined  above  as  discussed  in  the  following  chapter.  The 
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effect  of  long  chain  polymers  on  suppressing  boundary  layer  turbu- 
lence has  been  investigated  for  pipe  flows  and  external  1 lows  ; methods 
of  characterizing  the  boundary  layer  velocity  profiles  with  polymers 
present  have  been  developed  and  methods  for  prediction  of  flow 
separation  have  been  postulated. 

The  object  of  the  research  described  herein  was  to  develop 
predictive  methods  for  design  of  polymer  ejection  system  for  bodies 
of  revolution  and  to  perform  experiments  to  verify  the  predictive 
methods.  Additionally,  a study  of  the  effects  of  polymers  on 
boundary  layer  separation  was  made.  Several  "half  body"  configura- 
tions were  used  in  the  testing,  differing  in  their  conical  to 
hemispherical  tail  configurations.  The  bodies  were  dropped  in  a 
twenty-foot  long  cylindrical  tank  and  allowed  to  achieve  a steady 
state  velocity.  On  three  of  the  bodies,  dye  was  ejected  in  the 
boundary  layer  and  photographic  records  of  the  growth  and  separa- 
tion of  the  boundary  layer  at  several  velocities  were  made.  For  these 
tests,  the  ambient  fluid  was  water  and  a water  polymer  mix 
(polymer  ocean)  of  several  concentrations.  Another  series  of 
tests  were  made  with  a body  having  a significantly  different 
internal  configuration.  This  body  has  the  capability  to  store 
and  eject  concentrated  polymer  solution  from  its  nose  and  with- 
draw boundary  layer  samples  at  four  stations  along  its  length. 

The  boundary  layer  samples  in  turn  were  either  wall  samples  or 
samples  at  three  heights  within  the  boundary  layer.  The  stored 
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polymer  was  contaminated  with  a fluorescent  dye  which  was  later  used 
In  analyzing  the  concentrations  from  these  boundary  layer  samples. 
Thus,  it  was  possible  to  measure  the  change  in  wall  concentration 
and  concentration  profile  in  the  boundary  layer,  measure  steady 
state  body  velocity  and  compare  simultaneously  with  computed  results. 
This  research,  therefore,  concentrates  on  the  prediction  and  veri- 
fication of  the  drag  reduction  mechanism  on  bodies  of  interest  for 
underwater  application  to  gain  insight  which  will  help  in  the 
design  of  vehicle  ejection  systems. 

B.  Selection  of  Experimental  Apparatus 

The  intended  research  could  be  carried  out  in  water  tunnels, 
towing  tanks,  ring  channels,  in  buoyant  vehicle  experiments,  or  in 
drop  tank  tests. 

For  the  desired  research,  the  drop  tank  has  a number  of 
attractive  features.  Conceptually,  it  is  a relatively  simple 
experimental  tool.  The  ambient  turbulence  level  is  very  low, 
flow  visualization  is  easily  achieved  and  velocity  measurements 
readily  made.  Tests  over  a reasonable  range  of  Reynolds  numbers 
f are  achieved  by  weighting  the  body  being  dropped. 

One  of  the  chief  drawbacks  of  the  drop  tank  facility  as  com- 
pared to  a water  tunnel,  ring  channel  or  towing  tank  is  the 
inability  to  easily  measure  velocity  profiles  since  the  body 
being  measured  is  moving.  Since  flow  visualization  gives  an 
approximate  boundary  layer  thickness,  this  limitation  was  not  con- 
sidered to  be  an  important  drawback. 
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Fortunately,  the  author's  employer,  N'.tval  Underwater  Systems 
Center,  had  available  a drop  tank  with  suitable  dimensions  for  this 
research.  The  potential  for  meaningful  research  and  the  availability 
of  the  capital  equipment  led  to  the  decision  to  proceed  with 
ax isyinnio trie  boundary  layer  research  described  herein. 

C . Ri'su  1 t s 

Tests  were  performed  ir  a Reynolds  number  range  of  one  to  five 
million.  The  first  series  of  tc,;ts  were  performed  in  fresh  water 
and  In  a polymer  "ocean"  of  various  concentrations;  1.25,  2.5, 

5,  10,  20,  50,  and  60  WPPM.  Skin  friction  reduction  obtained  on 
all  bodies  approached  70  percent  at  polymer  concentrations  of  20 
WPPM  agreeing  well  with  data  from  other  experimenters.  Total  drag 
reductions  of  33  percent  for  the  6°  body,  16  percent  for  the  12° 
body,  and  10  percent  for  the  spherical  tall  body  were  obtained,  with 
the  effect  of  higher  percentage  form  drag  being  evident. 

Photographic  studies  of  dyed  boundary  layers  in  these  tests 
and  in  ejection  tests  with  50  WPPM,  500  WPPM,  and  1000  WPPM  displayed 
several  interesting  characteristics.  Addition  of  small  quantities 
of  polymer,  2.5  WPPM,  eliminated  the  fine-scale  turbulent  structure 
leaving  only  coarse  turbulence.  Higher  concentrations  resulted 
in  suppression  of  the  coarse  turbulence  and  extreme  thinning  of  the 
boundary  layer.  In  the  ejection  tests,  concentrations  of  500  WPPM 
and  1000  WPPM  displayed  no  turbulent  structure  or  mixing.  A dye 
streaking  phenomenon  with  spacing  equivalent  to  that  characteristic 
of  laminar  sublayer  streaks  was  displayed.  The  number  of  streaks 


was  also  approximately  equal  to  the  number  of  ejection  hole’s,  approxi 
mately  700,  supporting  the  fact  that  mixing  had  essentially  ceased. 

Concentration  profiles  were  measured  in  the  ejection  tests 
using  a fluorometric  method.  Tracer-contaminated  water,  5 WPPM, 

10  WP PM,  20  WPPM,  50  WPPM,  500  WPPM,  and  1000  WPPM  solutions, 
were  ejected.  Pour  measurement  stations,  the  last  at  an  x/L  of  .48, 
were  sampled.  The  water  data  displayed  expected  concentration  pro- 
files, agreeing  well  with  those  of  other  investigators.  At  all 
polymer  concentrations  tested,  unexpected  results  were  achieved. 

Wall  concentrations  remained  at  levels  predicted  by  molecular  diffu- 
sion. As  evidenced  from  the  photographs,  no  diffusion  was  occurring. 
An  initial  mixing  zone  neglected  by  most  investigators  was  the  con- 
trolling zone  in  these  tests.  The  data  for  this  axisymmetric  case 
were  compared  with  data  from  the  single  flat  plate  experiment  evi- 
dencing the  same  phenomenon.  A dimensionless  distance  developed 
in  the  flat  plate  experiment  used  to  predict  the  extent  of  this 
initial  zone  was  modified  and  compared  favorably  with  this  experi- 
ment. It  is  hypothesized  from  these  tests  that  optimal  ejection 
for  minimum  polymer  usage  requires  ejection  into  a laminar  boundary 
layer  prior  to  transition  to  turbulent  flow. 

Analytical  routines  are  developed  which  predict  boundary  layer 
parameters  and  polymer  wall  concentrations  for  this  postulated 
optimal  ejection  process  or  the  suboptimal  case.  Limited  verifi- 
cation of  the  model  is  made. 
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The  reduction  of  resistance  to  the  turbulent  shear  flows  of  liquids 
through  addition  ol  small  quantities  of  polymer  into  the  boundary  layer 
flowing  over  a body  has  received  considerable  attention.  The  first  re- 
portings occurring  independently  by  Mysel  (1949),  Toms  (1949),  and  Oldroyd 
(1949).  Although  extensive  theoretical  and  experimental  work  has  been 
carried  out  with  steady  pipe  flows  where  the  resistance  is  due  to  turbulent 
skin  friction  alone,  relatively  few  studies  have  been  conducted  on  the 
effect  of  additives  on  the  flow  about  axisymmetric  bodies  where  the  develop 
ment  of  the  boundarv  layer,  the  posit  ion  of  the  separation  point,  and 
the  unsteady  wake  significantly  contributes  to  the  total  resistance. 

Polymer  injection  in  developing  boundary  layers  has  only  recently  come 
under  study.  Polymer  ejection  type  data  for  external  flows  with  pres- 
sure gradients  are  virtually  non-existent.  Hoyt  (1972)  presents  an  en- 
cellent  summary  of  the  research  efforts  and  results  available  in  the 
field  of  drag  reduction. 

Literature  pertaining  to  the  current  study  will  be  reviewed  under 
four  separate  headings: 

1.  Bluff  body  experiments  with  polymer  solutions, 

2.  Turbulent  boundary  layer  theory, 

3.  Homogeneous  polymer  flow,  and 

4.  Polymer  ejection  studies  . 
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Bluff  Body  Experiments  W ith  Polymer  So  1 ut  i i>n 

The  earliest  known  experiment  with  bluff  bodies  in  drag  reducing 
polymers  was  carried  out  by  Crawford  and  Pruitt  (1963).  Additional 
experiments  with  different  sizes  and  shapes  of  bodies  (mostly  spheres 
cones  and  cylinders)  In  various  types  and  concent  rat  Ion  of  polymers  were 
conducted  by  others.  Work  performed  by  Lang  and  Patrick  (19bb)  is  re- 
presentative of  much  of  the  effort.  They  found,  as  did  many  of  the 
other  investigators,  that  the  drag  on  a sphere  was  considerably  reduced 
probably  due  to  the  rearward  movement  of  the  point  of  boundary  layer 
separation.  These  experiments  were  performed  in  a drop  tank  where  the 
fluid  ahead  of  the  ball  was  stationary.  A drag  reduction  of  69%  was 
attained  for  a 2-inch  diameter  steel  sphere  dropped  in  a 1000  WPPM  solu- 
tion of  Polyox  at  a water  Reynolds  number  of  1.4  x 10~\  It  was  also 
shown  that  the  polymer  produced  little  or  no  drag  reduction  on  stable 
bluff-based  bodies  whose  point  of  boundary  layer  separation  is  fixed. 

The  additives  produced  an  apparent  decrease  in  turbulent  mixing  in  the 
mean  wake  and  had  a tendency  to  display  a stringiness  in  the  wake  at 
high  polymer  concentration. 

Contradictory  results  were  achieved  by  Stow  and  Elliot,  1970,  who 
demonstrated  no  significant  drag  reduction  on  a tethered  ball  sus- 
pended in  a fully  developed  turbulent  pipe  flow.  This  is  not  surpris- 
ing, if  the  flow  about  the  sphere  exhibits  a turbulent  type  rearward 
movement  of  the  separation  point. 

The  performance  of  additives  in  reducing  drag  on  an  immersed  body 
is  complicated  by  a number  of  factors.  Certainly,  with  well  stream- 
lined bodies,  drag  reduction  is  effected  through  reduced  skin  friction 
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just  as  with  Internal  flow.  With  blunter  bodies,  however,  the  drag 
consists  primarily  of  form  dray,  and  Is  dominated  by  I lie  wake  sl/.e  and 
point  of  flow  separation,  skin  friction  frequently  being  negligible 
in  comparison.  On  typical  free-running  underwater  bodies  of  revolution, 
the  form  drag  can  represent  25%  of  the  total  drag.  The  two  investiga- 
tions previously  reported  demonstrate  the  importance  of  the  "free  stream" 
condition  on  the  movement  of  the  separation  point  with  polymers.  The 
effect  of  polymers  on  boundary  layer  separation  for  bodies  of  revolu- 
tion warrants  further  study  to  define  if  reduction  in  the  form  drag  as 
well  as  the  skin  friction  drag  can  be  affected  by  application  of  polymers. 

Turbulent  Boundary  Layer  Theory 

Since  the  concept  of  transit  ion  from  laminar  to  turbulent  flow  was 
demonstrated  by  Osborne  Reynolds,  in  his  classic  original  transition 
experiment  in  pipe  flow,  many  attempts  have  been  made  to  predict  the 
conditions  at  which  laminar  to  turbulent  flow  transition  will  take  place. 

Many  attempts  have  also  been  made  to  predict  velocity  profiles  and  wall 
shear.  So  many,  in  fact,  that  in  1968  the  Stanford  conference  was 
called  in  which  a total  of  29  methods  for  performing  turbulent  boundary 
layer  analysis  were  graded.  These  proceedings  have  been  edited  by 
Kline  (1968).  Such  general  analysis  can  be  divided  into  two  types: 

(1)  integral  methods  averaged  across  the  boundary  layer,  and  (2)  finite 
difference,  or  differential  methods  which  attempt  to  solve  the  full  par- 
tial differential  equations  of  the  boundary  layer.  A text  by  F.  M.  White 
(1974)  neatly  summarizes  and  discusses  many  of  the  competing  analysis  methods. 
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lis  considering  the  relative  Importance  el  visions  ami  t iirhu  leu! 

shear,  the  velocity  distribution  In  .1  turbulent  boundary  layer  follows 

a three- layer  concept.  The  layers  are: 

Inner  haver:  Viscous  shear  dominates 

Outer  Layer:  Turbulent  shear  dominates 

Overlap  Layer:  Both  types  of  shear  important 

The  mean  velocity  distribution  in  a two-dimensional  turbulent  boundary 

layer,  u(v),  depends  upon  four  local  parameters  — i , the  local  wall 

w 

shear  stress;  i , the  fluid  density;  p,  the  fluid  viscosity;  and  A,  the 
boundary  layer  thickness.  Brandt  1 (1913)  deduced  for  the  inner  law 
that  the  mean  velocity  did  not  depend  on  A resulting  in  a functional 
expression  for  the  Inner  law  of 

u = f (t  , p , u,  y),  ( 1 ) 

w 

Karman  (1930)  deduced  that  the  wall  acts  as  a source  of  retardation 
reducing  the  local  velocity,  u,  below  the  freestream  velocity  U^,  in 
a manner  independent  of  viscosity,  |i.  The  outer  or  velocity  defect 
relationship  is  then 

U-u=f('t,p,  y,  6).  (2) 

e w 

Coles  (1954)  performed  a dimensional  analysis  on  these  relationships 
resulting  in 

Inner  Law  — - = f ) (3) 

v*  v 

* U -u 

Outer  Law  6 - - = g(y/6)  (4) 

v* 

where  v*  is  a characteristic  velocity  called  the  wall  shear  velocity 
defined  as 

x L 

v*  = (-^)  \ (5) 

P 
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By  equating  equations  (3)  and  (4)  In  the  overlap  region,  one  may  obtain 
tlie  commonly  known  relation  for  the  "law  of  the  wall"  based  on  Inner 
variables 


-a-  • ' in  *>1*  + R 

V*  k v 


where  K and  B are  constants  determined  by  the  data  of  Nikuradse  (1930), 
to  be  .4  and  3.5,  respectively.  Coles  (1954)  displays  the  correlation 
of  the  dimensionless  velocity  profile,  u/v*  with  the  dimensionless  dis- 
tance for  the  wall  — — - up  to  a wake  of  about  300.  With  the  exception 

v 

of  separating  flows,  ail  the  data  nicely  collapse  into  regions  defined 
by  the  inner  and  logaritlunic  laws  having  dimensions 


Inner  Law 


Logaritlunic  I. aw 


- iL  = ^ 0 < .10 

v*  v v 


= 2.5  In  XV~  + 5.5  35  < ^ <300.  (8) 

M V 


For  regions  beyond  ^-V  = 300,  the  outer  region  or  wake  region. 

Coles  (195b)  postulates  a function  defined  as 

¥ W(6> 

be  added  to  the  equation  (6)  resulting  in 

u+  = i ln  y+  + B + ¥ w(i> 

where  II  is  related  to  Clauser's  (1954,  1955)  equilibrium  parameter 
x*  dP 

B=—-~.  0 

T dx 
w 

A reasonable  fit  to  data  is  given  by  the  relationship 

II  ; 0.8  (0  + 0. 5)°' 75  . (! 


Colt'S  ( 1 9 56)  proposed  the  following  curve  fit  to  the  wake  function  W. 

W#  = 2 (,2) 

Thus,  for  two-dimensional  flow,  empirical  relations  for  the  velocity 
profile  in  the  boundary  layer  are  available.  These  relations  will 
also  hold  for  thin  ax isymmet ric  boundary  layers  where  <5<<r0* 

Rao  (1967)  observed  that  in  the  viscous  sublayer  of  a thick  axi- 
syrametric  turbulent  boundary  layer  the  momentum  equation  reduces  to 

£ («£>•■ 

Integrating  with  u = 0 at  r = r , an  inner  law  for  the  sublayer  re- 


sults 


+ + . r 

u = r In  — 
o r 

o 


(14) 


where  r = r^  + y. 


When  JL  «i,  the  familiar  u+  = y+  relationship  returns  since  the  ln(l  + 


-^-)  approaches  y /rQ- 
o 


Rao  also  postulated  that  by  defining 


+ + . r (15) 

Y = r In  — 
o r 

o 

and  inserting  Y+  into  the  familiar  law  of  the  wall  (6)  in  place  of 
y+,  a correlation  with  the  two-dimensional  case  would  result.  The 
final  relations  are 


and 


„+  = v+ 


u+  = I ln  Y+  + B 
K 


Y+<5.0 


Y+>30. 


(16) 

(17) 
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Figure  1,  reproduced  Iron,  Ruols  paper,  displays  this  correlation  will 
data  obtained  by  Richmond  (1997)  and  Yu  (1998).  As  the  body  curvature 
effects  become  greater,  the  agreement  Improves  since  the  wake  tends 


to  disappear.  For  axisymmetric  cases  where  y/^  may  no.  be  considerably 


less  than  1,  equation  (19)  must  be  corrected  to  account  for  the  outer 
wake  region.  White,  et  a>  (1972),  accounts  for  the  wake  region  using 
an  inner  variable  approach.  The  basic  idea  of  the  method  is  that  in 


inner  variables,  the  velocity  profile  is  approximately  a logarithmic 


function  of  Y+  plus  a wake  which  depends  on  Y+  and  a single  dimension- 


less parameter  f,(x).  For  example, 


u+  * i in  y+  + B + f[  y+,  f.(x)  i 
K 


(18) 


for  the  two-dimensional  case.  The  derivation  of  the  functional  rela- 
tion for  the  axisymmetric  case  is  carried  tlnough  heri  , sinti  tlu. 
resulting  velocity  relationship  is  not  in  a lorm  readily  recognizable. 
Recalling  Prandtl's  mixing  length  theory: 

09) 


Du 

T : !'Ty’ 


Consistent  with  equation  (17)  for  the  axisymmetric  case,  f,  the  eddy 


/iscosity,  equals 
i Du  , 


t = p e?  I 


(20) 


3Y 


where  1 is  the  mixing  length. 


Now 


Du 

DY 


r Du 


(21) 


r 3y 
o 


substituting  in  (20) 


Fairly  near  the  wall  in  the  overlap  layer 
l = KY 

resulting  in  an  eddy  viscosity  which  is 


(23) 


. . a 


(2A) 


Finally,  upon  substitution  in  (19): 

2 


o '■>  r y 3u  \ 

, - p K>Y'~  (^> 

o 


(25) 


Equation  (25)  may  be  used  to  derive  an  axi symmetric  law  of  the 
wall,  using  a Couette  flow  approximation  for  the  near  wall  momentum. 


dP, 


(£)(")  ' r ^ = °' 


(26) 


Integrating  with 


r = r„  * rT  = Vw  ’ 


r = r , rt  = rt, 


1 dPe  y 7 

r t - r t + -x  (~j~" ) (r 
o w 2 ax 


*!>■ 


(27) 


Substituting  (25)  into  (27)  gives 


(28) 


The  pressure  gradient  parameter  a is  defined  as 


dp 

_ f ^ \ e 

a 'tv*  dx 
w 


(29) 


dp. 


Solving  (29)  for  and  substituting  in  equation  (28)  and  rearranging 

dx 


results  in 
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s'. 

r- 


r , , 1 v* 
t Ll  + ~ r 


r?  „3u  2 


w 


2 o v 


a (~t  - 1)1  = p KY2  ^ O 


w 


r2  '3y' 
o 


(30) 


Equation  (30)  may  be  placed  in  law  of  the  wall  variables  by  applying 
equations  (3),  (15),  and  (16)  and  equalities  from  Figure  2. 

Equation  (30)  becomes 
+ 


-du 


dY+  KY 


1 + 2Y+/r+ 

^ [l  + t ro  (e  ° - ^ 


(31) 


This  equation  reduces  in  the  limit  of  large  radius  to  the  law  of  the 
wall 


u+  (Y+,  0,  0)  = j In  (Y+)  + B 


, + + -KB 

where  u = 0 at  Y = e 
o 


= 0.1108  for  K = 0.4 
B = 5.5. 

-f*  *4"  / 

Integrating  across  the  boundary  layer  from  Y^  = 0.1108  to  Yq  (the 
edge  of  the  boundary  layer)  yields 
Y+ 

+ 1 


4-  -4- 

2Y/r  h 


“ -it  f -if  D - f r+  (1  - )]'  dY+. 

K 0.1108  Y 2 ° 


(32) 


Equation  (32)  represents  the  law  of  the  vjall  relation  for  the 
velocity  profile  accounting  for  pressure  gradient.  This  relation  holds 
throughout  the  boundary  layer  with  the  exception  of  the  viscous  sublayer 
region.  Figure  3 reproduced  from  White's  paper  shows  some  velocity 
profiles  obtained  by  integrating  equation  (32).  The  effect  of  positive 
(adverse  pressure  gradient)  raising  the  curve  above  the  incompressible 
log  law  and  the  effect  of  changes  in  r1  with  finite  a are  evident. 
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As  mentioned  previously,  many  methods  have  been  postulated  to 
predict  skin  friction.  White,  F.  M.  et  al  (1972,  1971)  extended  an 
earlier  integral  technique  developed  for  flat  plate  flows  (White,  F. 
M. , 1968)  to  the  case  of  thick  ax isymmetr ic  boundary  layer.  The 
method  concentrates  on  the  boundary  layer  equations  using  inner 
law  variables  and  is  interesting  in  that  boundary  layer  separation 
is  automatically  predicted.  The  integral  technique  provides  for 
relatively  simple  computer  solution  minimizing  required  computer 
time  and  has  been  demonstrated  by  White,  to  quite  accurately  pre- 
dict skin  friction.  A drawback  of  the  method,  as  with  most  integral 
methods,  is  that  thickness  calculations  are  not  very  accurate.  This 
is  not  a serious  drawback  for  the  intended  calculations  since  non- 
dimensionalized  boundary  layer  thickness  can  be  applied.  A priori 
specification  of  the  velocity  profile  is  required  with  this  integral 
technique,  as  with  others.  The  boundary  layer  continuity  and  momen- 
tum equations  for  turbulent  axisymmetric  flow  are  given  by 


(pur)  + — (pvr) 
3x  3y 


,9u.  . ,3u. 

pur  (— ) + pvr  (— ) 


- r (— ^)  + f-  (rr). 
dx  9y 


where  x and  y are  normal  to  the  body  surface  as  shown  in  Figure  2.  For 


the  portions  of  the  body  to  which  the  equations  are  applied,  the  angle 


0 between  the  normal  to  the  surface  and  the  normal  to  the  body  axis 
is  small  and  is,  therefore,  neglected.  Assuming  the  law  of  the  wall 


7V:-my' 


(32)  is  valid  across  the  entire  boundary  layer,  equations  (32)  through 
(34)  form  a closed  system  that  can  be  solved  for  the  skin  friction. 

The  oniy  additional  relations  required  for  the  present  study  are  those 
necessary  to  provide  a definition  of  the  polymer  wall  concentration. 


Appendix  A presents  a complete  development  of  the  equations  re- 
lated to  this  method  along  with  the  additional  terms  required  for 
solution  of  the  polymer  ejection  and  dilution. 

The  resulting  boundary  layer  equation,  after  considerable  alge- 
braic manipulation,  is 


+ r <i>  h - - v 


A dx*  1 


where 


p 2 2 

1 - 3 ■ & 1 • 


x*  * !■ 


V-IT’ 

O 


U L 

R _o_  , 

u V 


v dp 
_ / w . e 

a tv*  dx  ’ 
w 


and  the  definition  of  G , H,  and  I may  be  found  in  Appendix  A. 
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Other  integral  and  finite  difference  techniques  for  solution  of 
the  boundary  layer  problem  although  applicable  will  not  be  reviewed 
since  they  are  not  of  central  interest  to  the  work  being  presented. 

Homogeneous  Polymer  Flow 

Since  the  early  experiments  of  Toms  and  Olroyd  (1949),  a large 
number  of  high  molecular  weight  polymers  have  been  shown  to  be  effec- 
tive drag  reducers.  Polysaccharides  (Guar),  polyethylene  oxide, 
polyacrylamides,  and  sodium  carboxymethy 1 cellulose  have  received 
most  attention.  The  data  of  Hoyt  and  Fabula  (1964)  and  Virk  (1971) 
show  that  there  is  a maximum  drag  reduction  asymptote.  This  asymp- 
tote, for  a smooth  pipe,  corresponds  to  80%  of  the  friction  reduction 
that  would  be  attained  if  completely  laminar  flow  were  sustained  at  a 
given  Reynolds  number. 

The  efficiency  of  several  of  these  high  molecular  weight  polymers 
is  evidenced  in  Table  1 from  Hoyt  (1972). 


TABLE  1 

Concentrations  (WPPM)  of  Material  Required  to  Achieve  67-Percent 
Drag  Reduction  in  Pipe  Flow  at  Re=14xl03  (from  Hoyt  (1972)) 


Guar 

Gum  Karaya 
Polyox  WSR  301 


400 

850 

10 


Polyacrylamide,  Polyhall-250  20 
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Meyer  (1966)  and  Elata,  et  al  (1966)  have  shown  that  drag  re- 
duction in  pipes  is  due  to  a thickening  of  the  laminar  sublayer. 

It  was  shown  that  the  constant  B in  the  law  of  the  wall  (6)  to  remain  con- 
stant and  equal  to  the  Newtonian  value  until  a critical  threshold 

value  of  the  shear  velocity  was  reached  (v*)  after  which  B increased 

o 


logarithmically  with  v* 

B = 5. 5 + y In  (^£) 
v 


where  AB  = y In 


v* 


(Al) 


(A  2) 


White,  F.  M.  (1968)  has  plotted  data  from  several  investigators  and 
found  that 

1/2 


Y = 2.3  C 


w • 


(A3) 


Further,  the  data  indicate  a maximum  value  of  y of  approximately  1 1 

k 

and  a critical  shear  velocity,  v^,  for  the  onset  of  Jrag  reduction  of 
.08  ft/sec.  Many  authors  have  described  drag  reduction  as  a "nega- 
tive roughness  effect"  since  the  polymers  appear  to  thicken  the  sub- 

/ yv* 

layer  while  maintaining  the  same  slope  of  the  u/v*vs  In  curve 
in  the  overlap  region. 

In  the  outer  region,  where  the  wake  law  holds,  polymer  additives 
appear  to  have  no  effect  on  the  flow. 

Nadolink  (1968)  demonstrated  the  existence  of  the  thickened  sub- 
layer directly,  using  a high  speed  motion  picture  camera  and  a micro- 
scope . 
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Vlrk  (1966)  presents  evidence  Indicating  that 


the  onset  shear 


stress  Ls  inversely  proportional  to  polymer  molecular  r.id  ius  of  gyra- 


tion. Based  on  experimental  data,  the  critical  wall  shear  stress. 


t , which  must  be  exceeded  for  drag  redact  ion  to  occur  is  given  by 

w 


I = ,1(0.625  X 101'  p/R  ) = 

W ( i 


constant 


R£ 


(46) 


where  R is  the  rms  radius  of  gyration  of  the  molecule  as  deduced  from 


light  scattering  data. 

Fabula  et  al  (1969)  have  pointed  out  that  the  small,  value  of  the 
ratio  of  the  polymer  molecule  scale  to  the  scale  of  the  turbulent 
eddier  at  onset  indicates  that  individual  molecules  are  too  small 
by  several  orders  of  magnitude  to  Interfere  with  the  turbulence 
structure.  To  circumvent  this  problem  of  length  scales,  Fabula 
postulated  an  Interaction  between  the  time  scales  of  the  periodic 
molecular  deformation  in  the  viscous  sublayer,  given  by  y/2h  , where 
Y is  the  shear  rate,  and  a molecular  characteristic  relaxation  time 
tj.  t j is  determined  by  the  Zimm  or  Rouse  theories  which  relate  a 


characteristic  relaxation  time  of  the  solution,  t , to  the  solvent 


viscosity,  p , the  solution  viscosity,  p,  the  polymer  molecular 
s 


weight,  M,  temperature,  T,  and  the  concentration,  C 


t = a 
1 


(p-pg)M 


CRT 


(45) 


where  a is  a constant  having  a value  between  .4  and  .6. 

The  result  is  a criteria  for  the  critical  wall  stress  for  onset 
given  by 
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(4b  ■ 


2n(-'r/2n)  t ^ = 

or 

c ii  liRT 

i - - — — — . V* 

w 1 , aii  iiM 
* s 

Equation  (47)  may  be  rewritten  applying  a relation  between  intrinsic 
viscosity,  molecular  weight  and  molecular  dimensions  resulting  in 


c constant  ,,0, 

aw  : Rj~- 

(j 

Equation  (48)  of  Fabula  provides  a better  estimate  of  magnitude 
of  onset  shear  stress  whereas  Virk's  equation  (44)  gives  a better 
representation  when  the  constant  is  determined  by  a best  fit  of  the 
data . 

As  an  explanation  of  the  effectiveness  of  very  dilute  solutions, 
Fabula  postulated  that  entanglements  or  "blobs"  of  macromolecules, 
rather  than  individual  molecules,  are  responsible  for  drag  reduction. 
This  explanation  has  also  been  proposed  by  Kowalski  and  Brundret L 
(1974).  In  their  work,  a formula  has  been  developed  connecting  the 
size  of  the  entangled  molecules  with  the  size  of  a dissipative  eddy. 
The  macromolecules  entanglement  hypothesis  was  tested  to  predict  the 
so-called  onset  of  drag  reduction  in  pipe  flows  of  homogeneous  polymer 
solutions . 

Darby  (1972)  presents  a comprehensive  review  of  drag  reduction 
theories  in  which  comparison  of  molecular  hypotheses,  such  as  Virk's 
time  scale  hypotheses,  of  which  Fabula' s and  Kowalski's  are  examples; 
continuum  approaches  and  conventional  length  scale  boundary  layer 
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moil  i I ic.ii  ions  aro  made.  The  continuum  mec  han  i cs  approach  results  in 


tile  rone  1 us  ion  that  the  presence  of  elast  ii  properties  in  dilute  solu- 
tions being  a sufficient  criteria  for  drag  reduction.  This  is  accounted 
lor  hv  the  inclusion  of  a dimensionless  time  parameter  (the  Deborah 
number),  which  is  qualitatively  the  ratio  ot  a characteristic  time 
of  t ho  f luid  to  a characteristic  t inie  of  the  flow  system.  A presen- 
tation of  a dozen  different  forms  of  the  Deborah  number  is  made 
illustrating  the  problems  encountered  with  the  present  "state-of- 
the-art"  ot  this  approach. 

Transition  delay  from  laminar  to  turbulent  flow  has  been  reported 
by  White  and  McEligot  (1970)  and  is  found  to  depend  on  where  Lire  on- 
set shear  stress  is  reached.  If  tin*  onset  shear  stress  occurs  in 
the  Laminar  flow  region,  a delay  in  transition  to  turbulent  flow  can 


Extension  of  the  pipe  flow  data  to  provide  insight  into  How 
over  fiat  plates  has  been  successfully  performed  by  many  investiga- 
tors. (Iranvitle  (1971)  has  computed  the  maximum  possible  drag  reduc 
Lion  on  a flat  plate  of  80%  at  Reynolds  numbers  of  1 0! . Data  on 
flow  over  cylindrical  bodies  are  somewhat  sparse  hut  evidence  of 
35%  drag  reduction  by  Nadolink  ot  al  (1968)  for  a cylindrical  body 
dropped  in  a homogeneous  polymer  solution  is  representative  of  the 
results  to  be  expected. 


Polymer  Ejection  Studies 

Practical  applications  of  polymer  solutions  imply  the  necessity 
for  injection  whether  it  be  for  internal  or  external  flows.  Since 
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the  polymers  have  been  shown  to  effect  the  sublayer  resulting  in  drag 
reduction,  the  mechanism  of  injection,  the  diffusion  process  and  the 
quantities  of  polymers  required  as  compared  to  the  homogeneous  flow 
cases  treated  previously  are  of  extreme  practical  importance. 

Two  approximate  methods  have  been  used  to  treat  diffusion  patterns 
in  turbulent  shear  flows.  The  most  widely  known  method  employs  the 
eddy  diffusivity  model  which  assumes  that  the  flux  of  the  diffused 
matter  by  the  turbulent  fluctuations  Is  proportional  to  an  eddy 
diffusivity  term,  D , multiplied  by  the  local  concentration  gradient 

q = - D ^ . (49) 

y e dy 

Poreh  and  Hsu  (1971)  point  out  that  the  assumption  in  this  model  is 

that  D is  a function  of  the  flow  field.  The  fact  that  its  value  at 
e 

a point  can  be  specified  regardless  of  the  position  of  the  source 
only  holds  true  for  distances  from  the  source  that  are  large  compared 
to  the  lagrangian  integral  scale.  Measurements  made  by  Poreh  and 
Cermak  (1964)  indicate  that  this  characteristic  holds  for  turbulent 
shear  flows  and  the  Lagrangian  integral  scale  is  of  the  order  of  10 
boundary  layer  thicknesses. 

The  second  method  is  based  on  Batchelor's  (1957)  lagrangian 
similarity  hypothesis.  The  hypothesis  is  used  to  predict  the  turbu- 
lent motion  of  particles  in  steady,  self  preserving  shear  flows. 

Cermak  (1962)  applied  the  Lagrangian  similarity  hypothesis  to  predict 
diffusion  from  a continuous  point  and  line  source.  He  concluded  that 
results  from  application  of  the  Lagrangian  similarity  hypothesis  were 
significant  for  the  modeling  of  diffusion. 
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Hsu  (1971)  and  Poreh  and  Hsu  (1971)  applied  the  technique  ti 


predict  die  diffusion  boundary  layer  growth  in  the  intermediate. 


transition  and  final  zones  for  polymer  flows.  Applying  the  Lagran- 


gian  similarity  hypothesis,  the  following  equations  resulted  describing 


the  change  of  the  mean  vertical  position  y,  and  the  mean  longitudinal 


position  x,  for  an  ensemble  of  single  particle  releases 


where  b is  Batchelor's  constant,  and 


dx  _ 

d7  " u(7)' 


Combining  yields 


dy  _ bv* 
dx  u (Y)  * 


Ell ison  (1959)  estimated  Batchelor's  constant  to  be  given  by  b = K. 


His  analysis  further  suggests  that  y,  the  mean  position  of  particles 


at  a given  cross  section  x,  is  equal  to  Y,  the  mean  position  of  single 


particle  releases  when  x = x.  The  mean  position  of  particle  at  any 


cross  section  x may  be  defined  as 


l cy  dy 


l C dy 


By  replacing  y with  y/6,  where  6,  equals  the  value  of  y when  C/C 

d d max 


is  . 5 and  c by 


r~ = f(tK 

max  d 


Equation  (53)  may  be  Integrated  to  yield 


"i  -\i- 


(V.) 


Substituting  (54)  Into  (52)  results  In  an  expression  for  the  develop- 


ment of  the  diffusion  boundary  layer  with  distance  x 

d 6,  , . 

d bv* 


1 dx 


u (Y  ) 


(55) 


Hsu  (1971)  has  found  that  better  agreement  with  data  occurs  if  the 
constant  b in  (47)  is  replaced  with 

b = K(1  - |)  . 


(56) 


Figure  4 displays  a plot  showing  the  growth  of  the  diffusion  boundary 
layer  within  the  momentum  boundary  layer  and  depicting  the  several 
zones  of  diffusion. 

The  above  method  was  applied  assuming  similarity  in  concentration 
profiles  with  those  observed  by  Poreh  and  Cermak  (1964)  for  two-dimen- 
sional turbulent  mixing  of  ammonia  gas  from  a wall  line  source.  Porch 
and  Cermak  envisioned  a four-zone  process  which  is  worthy  of  descrip- 
tion since  the  concentration  profiles  have  come  under  much  study  as 
being  representative  for  polymer  diffusion.  The  four  zones  as  de- 
fined by  Poreh  and  Cermak  are:  (1)  the  initial  zone  - very  little 

reliable  data  were  obtained  in  this  region  due  to  very  large  velocity 
and  concentration  gradients.  The  extent  of  the  zone  (x/6)  was  not 
determined.  (2)  The  intermediate  zone  - the  diffusing  plume  is 
submerged  in  the  boundary  layer,  but  its  thickness  is  large  compared 
to  that  of  the  laminar  sublayer.  Diffusion  in  this  zone  depends  only 
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slightly  on  the  rate  of  the  boundary  layer  growth.  The  upper  limit 
of  the  zone  Is  about  x/4  = 18.  O)  The  transition  zone  - this  zone 
provides  lor  a decrease  in  the  rale  of  growth  of  the  diffusing  plane 
and  to  gradually  change  the  shape  ol  the  concentration  profile. 

X X 

Within  this  zone  18  < < 80.  Downstream  of  this  zone,  — > 60,  the 

6 6 

diffusion  plane  grows  at  the  same  rate  as  the  boundary  layer.  (6) 
The  final  zone  - diffusion  is  again  controlled  by  molecular  action. 

Morkovin  (1963)  described  the  data  taken  by  Poreh  and  Cermak 
as  shown  below: 
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Figure  5 displays  a plot  of  the  concentration  profiles  in  the  inter- 
mediate and  final  zones. 

Wells  (1968)  suggested  uniform  injection  through  a porous  wall 
since  it  raises  the  additive  concentration  to  the  drag  reducing  level 
in  the  wall  region  only.  Using  a Reynolds  - Prandtl  analogy  to 
analyze  the  diffusion  process,  he  calculated  that  distributed  ejec- 
tion would  require  40  to  140  times  less  additive  than  slot  ejection 
for  equivalent  drag  reductions. 

Walters  and  Wells  (1971)  conducted  tests  using  uniform  ejection 
of  polymer  solution  through  sintered  stainless  steel  smooth  cylin- 
ders. Fully  developed  turbulent  flow  was  achieved  in  the  test  sec- 
tion. Concentration  and  velocity  profiles  measurements  were  made  at 
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several  different  downstream  distances.  Pressure  drops  were  also 
measured  In  the  facility  to  give  Information  on  the  friction  reduc- 
tion. The  concent  rat  1 on  measurements  were  made  by  using  a f I uoromct r lc 
method.  Uranine  B was  used  as  the  fluorescent  dye  rather  than  several 
other  dyes,  since  the  laboratory  apparatus  could  be  easily  deconta- 
minated. The  dye  was  mixed  with  the  ejection  fluid  and  then  the 
tracer  sample  concentration  measured  with  a fluorometer.  Accuracies 
of  the  measurement  instrumentation  of  two  parts  per  billion  are  re- 
ported. The  conclusions  from  the  study  are  very  interesting.  The 
study  showed  that  the  porous  wall  approach  to  injection  required  lower 
quantities  of  additives  for  equivalent  drag  reduction  than  by  slot 
ejection.  For  certain  conditions  of  high  polymer  mass  flux,  a wall 
friction  increase  possibly  due  to  the  higher  viscosity  in  this  region 
was  noted.  As  compared  to  water  injection,  a one  or  two  order  of 
magnitude  reduction  in  total  diffusivity  in  the  ejection  region  was 
evidenced.  Downstream  of  the  ejection,  an  order  of  magnitude  reduc- 
tion of  total  diffusivity  was  noted  along  with  a significant  reduc- 
tion in  wall  friction. 

Latto  and  Shen  (1970)  performed  an  experiment  of  slot  ejection 
over  a flat  plate  positioned  in  a flume.  Using  hot  film  anemometry, 
it  was  found  that  the  momentum  diffusivity  was  less  than  for  pure 
water.  The  velocity  of  injection  was  also  found  to  be  important. 

It  was  found  to  be  desirable  to  keep  the  injection  velocity  as  low 
as  possible,  and  tangential  to  the  surface. 
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Wetzel  and  Ripkin  (1970)  experimentally  studied  polymer  injection 
Into  a developing  boundary  layer  in  a 9-ft-wide  open  channel.  Polymer 
solutions  were  injected  parallel  to  the  flow  near  the  bottom  wall. 

Pitot  tubes  were  located  at  positions  16,  28,  and  40  feet  downstream 
of  the  injection  slot.  Measurements  indicated  similar  velocity  pro- 
files at  each  station.  Several  methods  of  polymer  concentration  de- 
termination were  investigated.  The  fluorometric  method  was  found  to 
give  satisfactory  results.  The  resulting  concentration  profiles  for 
water,  1000  WPPM,  2000  WPPM,  and  3000  WPPM  were  found  by  the  author  to 
be  in  good  agreement  with  the  curves  developed  by  Morkovin  (1963) 
described  previously. 

A maximum  drag  reduction  of  35%  was  achieved  over  40-ft  boundary 
length.  Further,  at  a distance  of  16  ft  from  the  slot,  greater  drag 
reduction  was  attained  for  the  low  quantities  ejected  than  for  the 
higher.  At  further  downstream  distances,  the  reverse  was  true.  The 
behavior  was  attributed  to  more  complete  mixing.  Large  wavering 
parallel  streakes  were  evidence  shortly  after  injection  when  the 
polymer  was  color  dyed.  The  streaking  was  reported  to  be  a secondary 
three-dimensional  vortex  motion  superimposed  on  the  two-dimensional  flow. 

Fabula  and  Burns  (1970)  invoked  the  negative  roughness  analogy  to 
flat  plat  flow  with  polymers  so  that  the  outer  layer  mean  velocity 
similarity  is  unaffected  by  friction  reduction.  The  similarity  law 
of  mixing  with  polymeric  friction  is  predicted  to  be  the  same  as  with- 
out polymeric  friction  reduction.  Proceeding  with  this  assumption,  a 
relation  for  calculating  the  local  additive  concentration  at  downstream 
stations  along  the  wall  was  developed.  The  study  in< ludes  the  effects 
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of  turbulent  interraittency  and  the  similarity  profile  for  the  local 
additive  concentration  as  formulated  by  Foreh  and  Cermak. 


I.essmann  (1970)  extended  this  work  to  the  case  of  a body  of 
revolution  with  the  final  result: 


^ ^ ^ ( /1+v  cos  0,  , 

QiCi  = 27’  ro  > U C ( r } dy 

o o 


(59) 


where  the  terms  are  in  accordance  with  Figure  2. 

Wu  and  Tulin  (1970)  presented  experimental  data  obtained  by  in- 
jection of  various  polymer  concentrations  along  a smooth  and  rough 
flat  plate.  The  general  conclusions  were  that  the  slot  ejection 
angle  should  be  small  with  respect  to  the  flow  direction  and  the 
slot  width  should  be  comparable  with  the  thickness  of  the  viscous 
sublayer.  It  was  also  shown  that  a large  drag  reduction  was  obtained 
by  ejecting  the  additive  solution  at  a rate  comparable  to  the  normal 
viscous  discharge.  The  choice  of  additive  concentration  of  the  ejected 
solution  is  governed  by  the  length  of  the  boundary  and  its  roughness. 
These  findings  suggest  that  smaller  amounts  of  additives  are  needed 
for  injection  than  are  usually  estimated. 

Tullis  and  Rarau  (1973)  studied  the  characteristics  of  mean  tur- 
bulent flow  in  the  entrance  region  of  a rough  pipe  for  water  flow  and 
for  polymer  injection  into  a boundary  layer.  A 12-inch  diameter, 
200-ft-long  steel  pipe  was  used  for  the  study.  Polymer  was  injected 
through  a perforated  wall  pipe  section.  Drag  reduction  of  up  to 
80%  in  the  fully  developed  region  and  90%  in  the  inlet  region  were 
measured.  Comparison  between  water  and  dilute  polymeric  solution 
injection  showed  that  polymer  concentration  profiles  developed  slower 
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than  that  of  dye  in  water,  indicating  1 owe  i d i i t u:.  i v i t y . The  inlet 


length  needed  for  I low  to  tolly  develop  was  loom!  to  be  greater  for 
polvmer  injected  flows  than  for  the  case  of  no  Injection. 

Fruman  and  Tulin  (1974)  performed  a study  of  diffusion  of  a 
thin  tangential  jet  of  polymer  solution  injected  into  the  turbulent 
boundary  layer  of  a flat  plate  suspended  in  a high  speed  channel. 

Free -st ream  Reynolds  number  in  the  order  of  3.6  x 10;  were  achieved. 
Drag  measurements  by  reluctance  force  guages  and  wall  concentrations 
measured  by  a light  intensity  dyed  additive  method  were  taken.  The 
concentration  distribution  along  the  wall  was  found  to  be  represented 
by  two  regions.  Within  the  first  region  the  wall  concentration  is 
practically  constant  and  equal  to  the  Injected  polymer  concentration. 
In  the  second  region,  the  concentration  varied  Inversely  with  the 
distance  from  the  Injection  silt.  The  length  of  tills  first  zone  was 
some  15  to  20  times  that  of  water  injection.  The  length  of  this  first 
zone  appears  to  be  directly  related  to  the  thickening  of  the  viscous 
sublayer,  the  reduction  of  shear  stress  and  the  decrease  of  molecular 
diffusivity.  The  data  taken  correlates  well  when  formulated  using 
a heat  transfer  analogy  to  the  temperature  distribution  over  a flat 
plate.  The  distribution  of  the  wall  temperature,  T,  as  given  by 
Seban  (I960) 


T p . v . 1 . 2 -0.8 

— = 25.0  (-i— 2 -)  (-) 

T±  pu  s 


(60) 


Tulin's  data  for  p=p  yields 
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III.  KXI'FK  I MKNTAI.  Al’I’AKATl  'S 


The  experimental  apparatus  used  in  this  research  consists  of 
a clear  plastic  drop  tank,  a launcher,  several  test  models  and 
appropriate  instrumentation  and  photographic  equipment  to  measure 
model  velocity,  to  deduce  boundary  layer  thickness  and  to  measure 
boundary  layer  concentration  profiles.  Figures  6 and  7 show  the 
major  elements  of  the  experimental  apparatus.  The  characteristic 
of  each  of  the  above  elements  is  discussed  below. 

Drop  Tank 

The  drop  tank  facility  consists  of  a 20-ft  long  clear  plastic 
cylinder  2 ft  in  diameter.  Several  viewing  collars  are  available 
as  noted  in  Figure  7 which  may  be  positioned  where  desired.  These 
viewing  collars,  when  filled  with  water,  minimize  optical  distor- 
tion due  to  the  tank  wall  curvature.  The  liquid  in  the  tank  is 
continuously  filtered,  excepting  during  testing  periods,  to  main- 
tain high  clarity.  Located  to  either  side  of  the  tank  are  instru- 
mentation platforms  which  may  be  set  in  any  vertical  position 
desired.  A shock  pad  consisting  of  several  layers  of  dense  rubber 
topped  with  rubberized  horsehair  pads  totalling  approximately 
1-ft  thick  is  installed  at  the  bottom  to  absorb  the  energy  of  models 
being  dropped.  An  .062-inch  diameter  stainless  steel  wire  is 
suspended  down  the  center  of  the  tank  with  the  other  end  weighted 
to  approximately  75  lbs.  All  models  were  guided  down  the  tank  by 
this  wire.  The  drag  of  the  wire  on  the  bodies  has  been  determined 
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Figure  6.  Experimental  apparatus 


by  i till'  I. til  .splii'ii-  expcr  iment  with  simil.u  Miide  bearings 
Installed  in  the  sphere.  The  bearing  drag  force  determined  is 
approximately  .1  lbs. 
l.uuncher_  AssembJj/ 

A launcher  was  deemed  necessary  In  order  to  obtain  steady  state 
velocities  in  the  short  drop  tank  and  avoid  the  problem  of  pre- 
dicting boundary  layers  under  a transient  condition.  Figure  b 
schematically  shows  the  launcher  components  and  Figure  8 displays 
the  actual  launcher  used  in  many  of  the  tests.  As  can  be  seen  in 
Figure  8,  the  launcher  is  required  to  accelerate  the  test  bodies 
very  quickly  in  a very  short  distance  resulting  in  high  accelera- 
tions and,  to  the  dismay  of  occupants  of  the  building,  the  launcher 
ram  must  stop  in  a shorter  distance  with  even  higher  deaccelerations. 

The  launching  problem  consists  of  accelerating  a body  of  approxi- 
mately 10  lbs  weight,  plus  accelerating  launcher  components  weighing 
approximately  20  lbs  to  a desired  velocity  of  30  ft/sec,  the  maxi- 
mum velocity  considered  for  the  tests  to  be  performed,  in  a distance 
of  24  inches.  Applying  the  impulse  momentum  relation 

/ F dt  = / M dv  (62) 

where  F = the  force  required  (neglecting  drag  on  the  body) 
t = time  for  acceleration 
v = velocity 

M = mass  of  launcher  components  and  body  under  acceleration. 
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8.  Launcher  assembly  photograph 


The  time  for  an  acceleration,  a,  to  30  ft/sec  velocity  in  2 feet 


is  given  by 


t = 1 
a 


Determining  the  acceleration 
1 2 _ 1 v2 

X — T-r  3 V _ — 


Rearranging  and  solving 


a = — — = 225  ft/sec^  . 


Substituting  in  (63)  yields 


t = .1333  sec  . 


Integrating  (62)  with  the  established  limits  gives 


t = 0 v = 0 


t = .1333  sec  v = 30  ft/sec 


.1333  F - (30) 

32.2 


F = 210  lbs. 


A steady  force,  therefore,  of  210  lbs  would  do  the  job.  Figure  9 
is  an  assembly  drawing  of  the  launcher  constructed.  The  chief 
components  of  the  launcher  come  from  a high  pressure  hydraulic 
damper  cylinder  normally  operated  at  3000  psi.  For  safety  reasons, 
this  was  satisfying  since  the  pressures  that  the  cylinder  would 
be  subjected  to  in  this  work  were  low  but  the  shock  loads  high. 

The  overall  cylinder  length  was  about  2 feet  with  a working  piston 
diameter  of  4.550  in.  Fitted  to  the  end  of  the  piston  was  a conical 
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adapter  which  would  mate  with  the  h°  tailcone  ol  r he  prime  test 
model.  The  total  piston  assembly  Including  upper  guide,  piston 
and  lower  extension  pieces  was  approximately  62  indies.  All  pieces 
were  designed  to  minimize  weight.  A pressure  over  the  piston  of 
14  psi  should  then  provide  the  necessary  acceleration  force.  The 
problem  was  to  get  the  working  fluid  safely  to  the  piston  quickly 
enough . 

The  remainder  of  the  launcher  system.  Figure  4,  consisted  of 
a high  pressure  nitrogen  supply  (2200  psi)  regulated  down  to  150 
psi,  stored  in  a 4-cublc  ft  storage  flask  whose  contents  in  turn 
were  allowed  to  pass  through  a hand  operated  "fire"  valve  through 
a high  flow  low  pressure  regulator  to  the  launcher  Itself.  All 
pressure  lines  were  maintained  as  large  as  possible  ( // 1 2 AN)  where 
equipments  permlted.  The  launcher  chamber  pressure  was  recorded 
on  a Sanborn  recorder  so  Chat  launch  velocities  could  be  varied  by 
varying  the  low  pressure  regulator  output  pressure  in  a predeter- 
mined fashion. 

Through  tests,  it  was  found  that  a much  higher  final  pressure 
was  required  to  achieve  30  ft/sec  ejection  velocity.  A final 
pressure  of  about  60  psi  was  required  to  account  for  the  inertia 
of  the  low  pressure  regulator  and  the  losses  in  ihe  delivery  line. 
The  launcher  has  successfully  been  used  in  several  hundred 


launches . 
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Velocity  Measurement  Instrumentation 

Hie  velocity  measurement  Instrumentation  consists  of  three 
laser  photo-diode  stations  with  known  separation  as  shown  in 
Figures  6 and  7.  The  lasers  are  5 milliwatt  helium  neon  lasers 
manufactured  by  Spectra  Physics  (Model  120).  Lasers  were  used 
rather  than  a conventional  light  source  to  insure  that  a sharp 
cutoff  of  light  by  the  body  would  occur.  This  was  especially 
critical  where  measurements  across  the  body  length  dimensions 
were  made.  The  silicon  photo-diodes  used  to  detect  the  light 
cutoff  were  United  Detector  Technology,  Model  PIN-020A.  Their 
response  time  is  in  the  order  of  5 nanoseconds,  well  exceeding 

the  requirements  of  this  investigation. 

The  outputs  of  these  diodes  when  amplified  was  used  to  trig- 
ger two  Hewlett  Packard  counters  and  a special  purpose  dual  counter. 
These  may  oe  seen  In  Figure  6 In  the  Instrumentation  rack.  The 
bottom  instrument  was  only  used  for  troubleshooting  work  with  the 
equipment.  The  Instrumentation  was  so  arranged  to  allow  measure 
meat  of  the  velocities,  referring  again  to  Figure  6,  between  stations 
1 and  2,  2 and  3,  1 and  3,  and  across  the  body  (nose  to  tall)  at 
station  3.  Four  velocities  In  all  were  measured  across  a total 
span  of  5.7  feet  to  less  than  2 feet.  For  all  except  the  last 
measurement,  only  the  breaking  of  the  light  beam  by  the  nose  of 
the  model  was  used  to  start  the  counters.  This  is  rather  Impor- 
tant since  although  the  model  travelled  down  the  center  of  the 
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tank  guided  by  a wire,  small  lateral  movements , {'specially  with  tin' 

6°  tapered  tall  model,  would  result  In  rather  large'  unknown  dis- 
crepancies in  effective  length  of  tin.'  model  (where  the  laser  Inter- 
cepts the  model  on  entrance  and  exit)  and  as  a result,  in  the  velo- 
city measurements.  Lateral  motions  did  not  seriously  effect  measure- 
ments taken  from  the  nose  of  the  model  at  subsequent  stations  since 
the  curvature  of  the  nose  is  gradual,  introducing  less  error,  and 
since  all  lasers  were  positioned  on  one  side  of  the  wire,  at  about 
.25  inches  offset,  resulting  in  approximately  the  same  lateral 
offset  at  each  station  if  such  occurred. 

Bou ndarv  Layer  Thickness  Measurements 

The  title  of  this  section  is  a misnomer  of  sorts  since  actual 
boundary  layer  measurements  in  a drop  tank  experiment  of  this  type 
could  only  be  performed  with  extreme  difficulty.  Especially  in  the 
size  of  free  falling  model  being  tested  (3-tnch  diameter  x 25  inches 
long).  In  fact,  what  has  been  done  is  eject  from  the  test  bodies 
a visible  opaque  dye  and  photograph  the  body  and  Its  "boundary 
layer"  at  a given  position.  The  assumption  is  that  for  a dilute 
polymer  solution,  the  molecular  diffuslvity  and  viscosity  are  verv 
close  to  those  of  the  solvent.  As  concluded  by  Fabula  and  Burns 
(1970)  for  a far  downstream  case,  the  similarity  law  of  mixing 
with  polymeric  friction  reduction  is  predicted  to  be  the  same  as 
it  is  without  polymeric  friction  reduction.  Two  types  of  dyes 
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were  used;  the  first,  a food  coloring,  gave  excellent  results 
and  photographic  records  but  reduced  the  number  of  tests  to  3 to 
4 per  day  due  to  the  resulting  opaqueness  of  the  tank.  The  second, 
phenolphthalein  solution,  was  used  as  an  alternative.  The  solution 
PH  was  kept  sufficiently  low,  below  10,  so  as  to  not  effect  the 
polymers  used  (this  was  verified  by  drag  reduction  measurements  - 
no  velocity  changes  noted  with  or  without  the  dye) . The  dye  gave 
excellent  photographic  results  and  disappeared  within  seconds  due 
to  the  mixing  with  the  water  and  the  resultant  lowering  of  the  dye 
PH.  It  is  assumed  for  the  dye  as  with  polymers  that  the  diffu- 
sivity  would  be  equal  to  that  of  the  solvent. 

Figures  6 and  7 show  the  camera  station.  A graphic  camera 
and  Polaroid,  4"  x 5"  colored  film  was  used  for  the  photographic 
records.  The  shutter  of  the  camera  was  fixed  open  just  prior  to 
launch  and  a signal  from  the  first  laser  station  adjustably  delayed 
from  the  instrumentation  console  would  trigger  three  flash  units 
exposing  the  film.  The  adjustable  delay  allowed  for  correction 
for  velocity  changes  on  different  tests  or  to  photograph  different 
segments  of  the  body  on  subsequent  but  similar  tests.  Velocity  of 
the  vehicle  on  a run-to-run  basis  for  the  same  conditions  was  found 
to  be  in  the  order  of  2%  allowing  for  accurate  presetting  of  the 
delay  by  a calibrated  vernier. 
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Four  ax  lsymine  t r I c models  we  tv  constructed  lor  us*  in  thu  tost 
program.  Thruu  of  these  are  < al  led  the  dye  ejecting  bodies 
and  the  fourth,  Che  polymer  ejecting  body.  Their  purposes  are 
different  in  that  the  dye  ejecting  bodies  a.re  used  tor  drap, 
measurements  and  for  boundary  Laver  photographs  and  the  polymer 
ejecting  body  is  used  for  polvrner  ejection,  boundary  layer  con- 
centration measurement  studies.  This  latter  body  could  also  be 
used,  when  loaded  with  visible:  dye,  for  boundary  layer  photographic 
studies.  The  models  all  had  a maximum  diameter  of  three  inches  and 
an  identical  forebody  design.  In  fact,  the  dye  ejecting  models  all 
used  the  same  forebody  with  replaceable  tall  pieces.  The  forebody 
chosen  for  application  is  a classical  "half  body"  with  a length-to- 
diameter  ratio  of  3.88.  The  classical  half  body  shape  is  defined  by 
placing  a source,  of  strength  m,  in  a uniform  stream.  The  result- 
ing shape  is  defined  by  the  equation 

(67) 


rh  = Ro  4 


D 


where  RQ  = 4 


and  the  coordinates  arc  defined  from  the  source  as  shown  In 
Figure  10.  The  body  normally  extends  to  infinity  but  was  termi- 
nated at  11.65  inches  for  the  models  constructed.  The  half  body 
was  chosen  for  use  since  if  presented  a very  streamlined,  tractable 
shape  for  the  prediction  of  pressure  coefficients.  The  equation  for 
the  bodies,  in  x and  rn  coordinates  will  be  given  in  subsequent 
sections . 
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Dv  t ■ ['. j_c  c L i n g Bodies 


Three  Jlftor.nl  conf  Ir.urnt  I mm  of  dye  ""re 


Each  uses  the  half  body  forebody  dlscuasod 


prcv iousl y 


„y  followed  by  a three-inch  lone  spacer  section  of  constant 


, . f.  » { i n i o r o s won1  3SS cinb  1 G d t o 

three-inch  diameter.  The  various  tail  pieces  we 


this  assembly.  The  tall  pieces  constructed  were  a 6»  tapered  tall 


(,so  included  angle)  terminating  In  a .5-lnch  diameter  sphere. 


„ l2o  lnpered  tall  (240  included  angle)  terminating  in  a . 5-inch 


diameter  sphere  and  a he, sis ■»" 


Figures  11,  12  and  13  are  photographs  of  the  three  test  models. 


The  equations  for  the  external  shape  of  the  bodies  measured  from 


the  nose  as  shown  in  Figure  10  are  given  by . 


Half  Body  Section 


- = 079  + 4 . 346x  - .1022x2  + .122x3  5.29x 

L o 


for  0 < x < .981  Inches 


r„  ‘ -924 


+ 3296x  - .0659x2  + .000555x1  . ,000161x4  (69) 


for  .981  < x <_  11.^5  inches. 


itant  Diameter  Section 


r0  = 3 inches 


for  11.65  - x 114.65  inches  (70) 


6«  Tapered  Tail 


r = - .100511  x + 2.9686 


for  14.65  < x 1 24.58  inches.  (71) 


ing  model 


modi.'  1 


Equations  (68)  through  (72)  dcs<  ribed  the-  boundary  for  the  6° 


in  conLour  nL 


remove 


19.65  inches 


Cone 


for  19.65  < x < 20.15  inches 


Equations  (68),  (69),  (70),  (73)  and  (74)  describe  the  external 


boundary  for  the  12°  model 


Constant  Diameter  Section 


15.132  inches.  (75) 


Hemisnher ica L Tail 


x < 16.632  inches, 


Equations  (68)  , (69) , (75)  and  (76)  form  the  equation  for  the 


hemispherical  tail  model 


All  bodies  were  constructed  of  aluminum  and  allowed  for  the 


addition  of  weights  for  matching  of  velocities  when  testing  the 


•vjyry!}-—**,; 


models  under  increased  drag  conditions.  Figures  1A  and  15  describe 
the  body  components  and  features.  Referring  to  the  figures,  an 
explanation  of  the  operation  of  the  model  will  b . describe  the 
function  of  the  various  pieces.  The  model  is  made  ready  for  a 
launch  by  inserting  the  tank  guide  wire  through  the  body  guide 
tube,  © the  wire  bears  on  a nylon  bearing  at  either  end  of  the 


model  to  minimize  friction.  The  tail,  ^10j,  °f  the  model  is 
removed  (separated  sufficiently)  and  the  internal  cavity  filled 
with  dye.  The  dye  ejection  ports  being  previously  taped  over 
to  insure  no  loss  of  dye.  The  model  is  reassembled  and  the  vent 
screw  (unmarked)  in  the  tail  is  removed  and  the  model  placed  under 
the  water  level  in  the  tank  to  back  fill  the  cavity  totally  with 
water.  The  model  may  now  be  fired  upon  removal  of  the  tape  over 
the  ejection  holes.  When  fired,  stagnation  pressure  enters  the 
A stagnation  ports,  © is  transmitted  through  the  body  guide 
tube  and  pressurizes  the  bladder, ^A^,  forcing  the  dye  out  the 
four  dye  ejection  ports  located  at  the  minimum  pressure  point 
of  the  forebody  as  determined  by  a potential  flow  program.  After 
several  tests,  the  bladder  was  found  to  be  not  necessary  simpli- 
fying the  operation.  The  stagnation  ports  consisted  of  8,  .050 
inch  diameter  holes  located  on  a .625  inch  circle  at  the  nose 
of  the  body.  The  four  ejection  ports  were  also  .050  inch  diame- 
ter placed  on  a 1.125-inch  diameter  circle,  normal  to  the  surface 
of  the  body.  Table  2 summarizes  the  pertinent  information  about 
the  bodies. 
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STAGNATION  PRESSURE  INLET  PORTS 

BODY  WEIGHTS 

DYE  STORAGE  VOLUME 

BLADDER  FOR  DYE  EXPULSION 

REPLACEABLE  TAIL  CONFIGURATIONS 

DYE  EJECTION  PORTS 

HALF  BODY  CONTOUR 

SPACER  SECTIONS 

BOD.’  GUIDE  TUBE 

REPLACEABLE  TAIL  CONFIGURATIONS 


Figure  14.  Dye  ejecting  axisynunetric  body  description 
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TABLE  2 

Dye  Ejecting  Ax lsymmet r 1c  Body  Information 


6°  TAIL 

12°  TAIL 

HEMISPHERIC 

TAIL 

LENGTH  (FT) 

2. 0833 

1.679 

1.  386 

DIAMETER  (FT) 

.25 

.25 

.25 

L/D 

8.33 

6.7 

5.54 

FRONTAL  AREA  (FT2) 

.0491 

.0491 

.0491 

'WEIGHT  IN  WATER  (LBS) 
WEIGHTED 

4.43 

7.109 

3.95 

UNWEIGHTED 

3.47 

2.94 

Polymer  Ejecting  Body 

The  polymer  ejecting  body  is  similar  in  external  configuration 
to  the  6°  dye  ejecting  body.  The  equations  for  the  external  shape 
hold  for  this  body.  The  body  nas  been  designed  to  perform  two 
functions  while  traversing  the  length  of  the  drop  tank.  These 
functions  are: 

1.  Eject  polymer  or  other  fluids  at  a constant  rate 

2.  Withdraw  samples  from  the  boundary  layer  at 
four  axial  stations 

Figures  16  and  17  describe  the  various  elements  of  the  ejecting 
body.  These  elements  may  be  viewed  in  assembled  and  disassembled 
form  in  Figures  18  and  19,  respectively.  The  overall  dimensions 
of  the  model  are  25  inches  in  length  by  3 inches  in  diameter. 

Prior  to  discussing  the  design  values'.^!  the  various  elements,  a 


POLYMER  EJECTION  SCREEN 

BOUNDARY  LAYER  SAMPLE  TUBES 

BOUNDARY  LAYER  SAMPLE  COLLECTION  CHAMBER 

POLYMER  SOLUTION  STORAGE  VOLUME 

SUCTION  PISTON  FOR  BOUNDARY  LAYER  SAMPLING 

POLYMER  EJECTOR  PISTON  AMBIENT  PRESSURE  BACKFILL  VOLUME 

POLYMER  EJECTOR  PISTON  PRESSURIZATION  AREA 

C02  ORIFICE  AND  PUNCTURE:  NEEDLE 

CO2  SUPPLY 

HALF  BODY  CONTOUR 

VOLUME  AVAILABLE  FOR  WEIGHT 

POLYMER  EJECTOR  PISTON 

POLYMER  EJECTOR  PISTON  STROKE  MEASUREMENT  ROD 
C02  SUPPLY  ACTIVATION  DEVICE 
BODY  GUIDE  TUBE 


Figure  16.  Polymer  ejecting  axisyiumetrie  body  description 


59 


. 1.  . I'.-Yl'i-.n-,,-.,  -a-.  ~ ikl-'i-S..  » , V.  HM.  J. .1,  


c-nib  led 


I 

I?; 


I 


ft 


I 


, 

w 


I 


w 


I 


I 


K 


£2 


■»  */—  -:*■  '«Tts- . -~  .-.- 


deser  iption  of  the  operation  of  the  model  will  serve  to  identify 
all  the  major  components.  Referring  to  Figures  16  and  17,  the 
model  is  assembled  for  a launch  with  polymer  or  other  fluid  in 
the  storage  volume,  © and  a fresh  energy  source,  ^9^,  a CO2 
cartridge  installed.  The  model  is  placed  on  the  guide  wire  in 
the  drop  tank,  ^lT),  and  inserted  in  the  launcher.  A nylon  shear 
screw  is  attached  through  the  launcher  to  the  polymer  ejector  piston 
stroke  measurement  rod,  . This  rod  serves  to  hold  the  model 

in  the  launcher  and,  upon  impulsing  the  launcher,  triggers  the 
CO2  supply  activation  device,  © , puncturing  the  CO2  cartridge 


and  releasing  the  constant  pressure  (800  psi)  gas  supply.  The 
model  leaves  the  launcher  shearing  the  nylon  screw  allowing  the 
spring  loaded  stroke  measurement  rod,  0,  to  bear  against  and 


follow  the  polymer  ejector  piston,  in  its  travel.  The 

gas  from  the  CO2  supply  passes  through  an  orifice,  © and 
pressurizes  the  annular  piston  area,  © of  the  polymer  ejector 
piston  driving  it  forward.  This  action  results  in  ejecting  the 
polymer  from  the  storage  volume,  © through  the  ejection  screen, 
© and  providing  a suction  in  the  annular  area,  ^jT^.  This 
area  is  connected  to  four  equal  sampling  chambers,  © through 
tubes  located  within  30°  of  the  arbitrarily  chosen  top  of  the 
model.  The  suction  then  draws  either  a wall  sample  or  a samplev 
from  a probe,  © As  the  piston,  ^lT),  is  moving,  ambient  pressure, 
which  varies  from  atmospheric  at  the  surface  of  the  tank  to  near  9 
psi  above  atmospheric  at  the  bottom  of  the  tank,  enters  the  guide 
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( ubo , ©•  ■ mil  Is  transmitted  to  the  ambient  pres. •.lire  hark  I 111. 
volume,  (^6^.  Finally,  at  the  end  of  the  piston,  (^)>  travel,  the 
CO2  gfls  Is  vented  to  the  backfill  volume  chamber,  0 and  out  the 
guide  tube,  © . The  model  at  the  end  of  a test  is  thoroughly 


vented  of  gas  and  readilv  handled. 


The  model  was  designed  with  several  thoughts  in  mind.  First, 


it  would  be  desirable  to  eject  at  at;  low  a velocity  as  practical 


so  as  not  to  disrupt  the  flow.  Secondly,  ejection  in  the  nose 


stagnation  region  forward  of  the  minimum  pressure  point  would  be 


desired  to  minimize  the  possibility  of  early  tripping  of  the  boundary 


layer  to  turbulent  flow.  The  volume  ejection  rate  should  result 


in  an  annular  flow  of  about  the  same  dimensions  as  the  laminar 


boundary  layer  thickness  to  minimize  the  disturbance  on  the  flow 


field.  The  stroke  and  suction  process  should  be  at  a constant 


rate.  The  suction  process  should  be  at  a low  enough  rate  to  mini- 


mize the  disturbance  on  the  boundary  layer  and,  for  Lhe  wall 


suction  case,  to  insure  a sample  is  withdrawn  near  the  wall  and 


not  penetrating  the  boundary  layer  too  far.  Finally,  all  processes 


should  be  completed  prior  to  the  model  impacting  the  bottom  so 


as  not  to  be  either  ejecting  or  sampling  under  a zero  velocity 


condition.  Conversely,  it  is  desirable  to  delay  sampling  until 


the  boundary  layer  has  been  established  at  launch.  This  require- 


ment was  difficult  to  achieve  and  attempts  discarded  early  In  the 


design  process.  The  errors  introduced  are  less  than  8%  for  the 


worse  case. 


The  first  design  parameter  fixed  was  that  of  deciding  on  the 
required  polymer  ejector  volume.  After  deciding  to  keep  the  model 
as  light  as  possible,  constructing  with  aluminum,  a guessed  weight 
in  water,  based  on  some  calculations  for  a hemispherical  nose,  a 
cylindrical  section  and  a tail,  of  A. A3  lbs  resulted  for  a 2A-inch 
long  body  of  3-inch  diameter.  Assuming  drag  coefficients  of  .22 
and  .11,  the  velocity  of  the  body  was  determined  by  the  equation: 

U = (■  -W~B )1/2  (77) 

1/2  CDpA 

The  calculated  velocities  for  70°  water  were 
U = 20.6  ft/sec  for  Cp  = .22 

U = 29.0  ft/sec  for  Cp  = .11. 

The  Reynolds  numbers  for  these  velocities  were 
Re  = 3.9  x 106  0 U = 20.6  ft/sec 
Re  = 5.5  x 106  0 U = 29.0  ft/sec. 

The  laminar  boundary  layer  thickness  may  be  calculated  using  flat 
plate  approximations  and  assuming  a position  2 inches  back  along 
the  body  surface.  From  Schlichting  (1968) 


6_  = 5^2 

XS  (Re) 


1/2 


(78) 


Now  R„  = 3.25  x 10- 


RP  = A. 6 x 10' 


Therefore , 


6 = 1.52  x 10~3  ft 


U = 20.0  ft/sec 
U = 29.0  ft/sec. 

U = 20.0  ft/sec  . 
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A - l.:»7  X 10  ‘ ft  IJ  = 29.0  ft/s.*r. 


Now,  the  volume  flow  through  Lhe  annulus  ariunul  the  body  Is  given 


Q6  = . 6 1;  H D A 


where  . 6U  Is  used  as  an  approximation  to  an  average  velocity  In 


Lhe  region.  For  the  two  conditions  being  calculated 


Q = .0148  ft^/sec  = 25.5  in^/sc*c 


for  U = 20.6  ft/sec 


Q - .0174  ft-Vsec  = 29.9  in^/sec 


for  U = 29.0  ft/sec. 


The  usable  tank  length  for  model  tests  is  about  13  ft;  therefore. 


the  time  of  ejection  will  be  about 


ti  = — — — = .63  seconds 

1 20.6 


t2  = .448  seconds 


Required  volumes  for  storage  are 


U8  = 16.0  inJ 


U = 20.6  ft/sec 


Us  = 13.4  in- 


U = 29.0  f t/sec 


Comparing  these  volume  flux  rates  with  expected  velocities  from 


the  nose,  assuming  a one-inch  diameter  ejection  port  with  a 307. 


open  area  screen  installed 


(Ae) ( . 30) 


= 9 ft/sec  for  U = 20.6  ft/sec 


Qd 

(Ae) (-30) 


= 10.6  ft/sec 


for  U = 29.0  ft/sec. 


' - ' ‘ 


_____ 


Tliivsc1  velocities  are  high.  Velocities  In  the  okUt  ol  > Tt/. 
be  much  more  desirable.  Reduction  of  velocities  by  1/2  would  -suit 
in  a storage  volume  goal  of  about  8 cubic  inches.  This  reduction 
would  also  result  in  smaller  thicknesses  of  flow  in  the  boundary 

layer  as  compared  to  the  laminar  layer  thickness. 

The  effective  polymer  storage  volume  designed  into  the  model 

is  8.25  in-*  with  a stroke  of  4.67  in.  and  a stroking  time 
secs . 

The  ejector  screen  used  In  the  model.  Figure  18,  was  formed 
to  the  contour  of  the  forebody  nose.  The  screen  contained  517 
.024  inch  diameter  holes  per  square  Inch  with  a 247.  open  area. 

The  outside  diameter  of  the  screen  Is  1.37  inches  and  the  inside 
diameter  .3  Inches.  The  projected  area  Is  1.4  In2  with  an  open 
orifice  area  of  .337  in2.  This  ejection  area  resulted  In  a 5-ft 

per  second  ejection  velocity. 

, . Nn  n straight  screen  from  Harrington 
The  screen  used  is  a Mo.  u scraigi 

and  King  Perforating  Company. 

The  suction  time  was  next  to  be  determined.  Since  the 

single  piston  approach  was  decided  upon,  the  stroke  available 

and  the  stroking  time  for  use  was  equivalent  to  that  for  the 

ejector  piston,  4.67  inch  and  a time  of  .4  seconds.  The  suction 

rate  for  the  design  of  the  suction  piston  is  11.675  in2/sec. 

Applying  a method  defined  by  Dowdell  ,1973)  for  correction  of 

the  measured  "wall  concentrations"  due  to  the  sampling  process 
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1 

4 


drawing  further  out  in  L 1 1 » ■ boundary  layer,  an  estimation  of  a 
reasonable  sample  rate  and  sample  probe  or  If  Lee  size  were  made. 

It  was  first  decided  that  samples  would  be  made  at  stations 
3,  6,  9,  and  Id  incites  back  from  the  nose  of  the  body.  This  selec- 
tion was  somewhat  governed  by  the  need  in  the  tail  for  other 
machinery  space.  From  a measurement  standpoint,  it  was  desira- 


m 


bio  to  maximize  sample  size  but  from  a correction  standpoint, 
minimal  flow  velocity  is  desirable.  As  a result,  four  suction 


? I 


I 

*1 


ports  90°  apart,  were  located  at  each  measurement  plane.  Ad- 
dition,ally,  at  each  subsequent  plane  the  port  locations  were 
rotated  22°  '30'  so  as  not  to  be  affected  by  the  flow  entering 
the  ports  in  the  previous  pla.e,  Figure  18.  The  next  step  was 
to  define  the  approximate  boundary  layer  thickness  at  the  most 
forward  probe  station  for  the  highest  Reynolds  number  case  as 
used  before.  This  would  give  the  thinnest  boundary  layer  which 
is  the  worst  case.  Using  a flat  plate  turbulent  boundary  layer 
thickness  expression  from  Schlichting  (1968) 

A = 37  (8i 

x (Re)1/5  ' 

From  (80),  6 = .05  inches. 

A matrix  of  probe  sizes  (inside  diameter)  and  suction  volumes 


H I 


" I 


' 


were  examined  for  the  piston  rates  calculated,  body  velocities 
and  the  boundary  layer  thickness  calculated  above.  This  was 


■K 


also  performed  for  a boundary  layer  thickness  of  .025  inches, 
postulated  to  be  the  thickness  that  might  occur  when  polymers 
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were  added.  The  result  was  that  a .760  cubic  inch  suction  volume 


and  sampling  tube  inside  diameter  of  .030  inches  would  give 

acceptable  results.  A worst  case  correction  factor  of 
/•» 

measured  _ _ 54 

n 

w actual 

was  calculated  for  a .025-inch  thick  boundary  layer.  This  correc- 
tion reduced  to  about  .3  at  a boundary  layer  thickness  of  .05 
inches.  Although  a lower  correction  factor  would  be  desirable, 
an  ample  sample  volume  of  fluid  is  necessary  for  concentration 
determination.  This  suction  rate  and  flush  wall  sample  tube  size 
had  an  inflow  rate  per  tube  of  approximately  14  ft/sec. 

The  pitot  probes  to  be  used  in  sampling  the  boundary  layer 
away  from  the  wall  were  designed  to  minimize  the  effects  on  the 
flow.  Two  heights  were  chosen  for  measurement,  .025  inches  from 
the  wall  and  .055  inches.  These  selections  came  about,  in  part, 
due  to  physical  limitations.  With  the  physical  dimensions  of 
the  probes,  it  is  difficult  to  get  closer  to  the  wall  and  with 
the  projected  boundary  layer  thickness  for  the  first  station,  it 
is  not  desirable  to  get  any  further  from  the  wall.  These  heights 
would  allow  for  a reasonable  three-point  plot  of  boundary  layer 
concentrations.  Figure  20  displays  the  salient  features  of  the 
probes.  As  is  noted,  a thin  wall  tubing  is  used  which  allows  for 
flattening  of  the  tip  to  .01-inch  thickness.  The  resulting 
inner  dimension  for  sampling  is  .004  inches  high.  Sampling 
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Figure  20.  Special  boundary  layer  probe3 
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velocities  in  the  tube  arc  approximately  20  ft/sec.  Sixteen 
probes  were  used  in  each  test  replacing  the  sixteen  flush  wall 
probes . 

The  model  design,  as  can  be  seen  from  the  schematic,  Figure  17, 
contained  void  volumes  in  the  suction  lines  which  were  prefilled 
with  liquid.  The  void  volume  is  1.4316  cubic  Inches  resulting 
in  a sample  dilution  factor  of  .1327. 

The  stroking  piston  pressurizat ion  system  consists  of  the 
carbon  dioxide  (CO?)  cartridge,  a puncture  device  and  an  orifice. 
Since  the  goal  is  a constant  rate  ejection  process,  the  use  of 
a CO2  cartridge  as  the  gas  source  was  chosen  since  the  substance 
is  stored  in  liquid  form  with  a vapor  pressure  of  800  psi  at 
about  65°F.  Little  change  in  the  vapor  pressure  would  occur 
during  the  use  cycle  as  the  gas  is  drawn  off.  The  system  was 
designed  to  minimize  the  free  volumes  in  the  high  pressure  area 
to  insure  a sufficient  gas  supply  for  the  stroking  process,  Figure 
17.  The  CO2  cartridge  in  the  final  version  of  the  model  is 
inverted  from  that  shown  in  Figure  16,  with  the  puncturing  needle 
placed  on  the  opposite  end  also.  This  was  done  to  achieve  repro- 
ducibility in  the  stroking  of  the  piston.  Apparently,  either 
freezing  of  the  orifice  or  slight  contaminants  would  plug  the 
metering  orifice  resulting  in  extreme  variability  of  stroking 
times.  The  reversal  corrected  the  situations.  The  pressuriza- 
tion area  and  orifice  area  size  for  metering  the  CC>2  gas  were 


71 


.selected  to  provide  the  stroking  rate  previously  determined.  The 
area  was  somewhat  governed  hy  the  desire  to  maximize  r h « ■ drop  tank 
depth  pressure  balance  area,  0 Figure  17.  The  volume  In  this 
area,  besides  providing  for  balancing  of  the  depth  pressure  also 
minimizes  the  weight  change  during  a launch  by  water  backfilling 
the.  volume.  The  COt  piston  pressurization  area  designed  was  .785 
in^  with  a resultant  total  volume  of  3.68  in~*.  The  orifice  size 
was  determined  to  achieve  the  .4  second  stroking  time.  Consider- 
ing that  the  orifice  is  small  enough  such  that  choked  flow  condi- 
tions exist,  the  pressure  ratio  across  the  nozzle  is  given  by, 
Shapiro  (1953), 

k 

(81) 


P2  , 2 


PI 


(.  k+1 

k+1 


Which  for  k = 1.29  for  CO2  gas  ('■  68°F  yields 
P2 


Pi 


= 0.55. 


(82) 


Therefore,  flow  will  be  a maximum  for  downstream  pressures  less 
than  or  equal  to  440  psi  with  an  upstream  pressure  of  800  psi. 
The  f Low  rate  is  given  by,  Shapiro  (1953) 


(--) 

A1  max. 
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(83) 


where  p = .114  lbs/ft^  @ 68°F. 

O 

Substituting  in  (83): 

— - = 4.28  x 10“  lbs/ft^  sec.  (84) 

Al 
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Convert  inj>  to  match  the  volume  flow  requirements 


- J£—  - 37.54  x LO2  ft/sec.  (85) 

A]  PrA] 

To  achieve  the  9.2  tn-^/sec  volume  stroking  rate  for  the  pressurizing 
gas  volume,  Q must  equal  9.2  in2/sec  and  is  determined.  Sub- 
stituting in  (85) 

A]_  = 1.418  x 10-6  ft2  (86) 

= 2.04  x 10  ^ in^ 

Tlierefore,  the  orifice  diameter  must  be  .016  inches.  This  diameter 
orifice  has  been  used  in  the  model. 

The  important  characteristics  of  the  6°  ejecting  model  as 
finally  designed  and  constructed  are  displayed  in  Table  3. 

TABLE  3 

Characteristics  of  6°  Ejecting  Model 


Length  2.0833  ft 

D .25  ft 

L/D  8.33 

Weight  in  Water  4.46  lbs 

3 

Polymer  Solution  Storage  Volume  8.25  in 

Stroke  4.67  in 

Stroke  Time  .400  sec 

Suction  Volume  .760  in2 

Suction  Volume  Per  Sampling  Chamber 

(4  Chambers)  .19  in 

Sample  Dilution  Factor  .1327 

Ejection  Velocity  5 ft/sec 
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Boundary.  J..iyer  Goncent  rat  ion  Measurement  s 

Several  methods  exist  for  the  measurement  of  the  concent  ra- 
tion of  polymer  in  a sample  solution.  These  Include  pol arographic , 
turbidity,  simple  opaqueness  and  fluoromctric  methods.  A compari- 
son of  all  but  the  opaqueness  is  given  Ln  Wetzel  and  Ripkin  (1970). 
The  opaqueness  method,  which  relies  upon  a measurement  of  light 
transmission  through  a portion  of  the  boundary  layer,  cannot  be 
applied  to  the  moving  body  case  since  short  time  measurements  .and 
accurate  body  position  with  time  would  have  to  be  known.  The 
remaining  methods  are  all  applicable  to  the  conditions  of  this 
study  since  they  deal  with  a sample  of  solution.  The  method 
chosen  for  use  is  the  fluorometric  method.  This  method  has 
been  used  successfully  bv  Walters  and  Wells  (1971),  Wetzel  and  Rip- 
kin  (1970)  and  Tullis  and  Ramu  (1973).  The  method  consists  of 
injecting  a tracer  dye  into  the  fluid  being  analyzed  and  capturing 
a sample  of  the  fluid  for  analysis.  This  assumes,  of  course,  that 
the  diffusion  of  the  tracer  dye  is  identical  to  that  of  the 
boundary  layer  mixture.  A reasonable  assumption  as  data  comparing 
various  measurement  methods  by  Wetzel  and  Ripkin  (1970)  display. 

Tire  tracer  dye,  which  may  be  colorless,  fluoresces  when  radiated 
with  an  ultraviolet  source.  This  reradiation  is  measured  with  a 
photomultiplier  tube  and  is  proportional  to  the  concentration  of 
dye  in  the  sample.  The  fluorometer  used  in  these  tests  is  a 
G.  K.  Turner  Associates  Model  111  Fluorometer. 
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Several  tracer  dyes  may  be  used.  Two  prime  candidates  are 
Rhodamlne-B  and  Uranine-B.  Both  are  compatible  with  polymer 
solutions  in  that  the  mixture  was  stable  and  the  drag  charac- 
teristics were  unaffected  when  measured  in  a Turbulent  Flow 
Rheometer  similar  to  that  described  in  Hoyt  (1966).  The  Rhodamine- 
B gives  excellent  results  with  solution  concentration  of  1 part 
in  10  inches  on  a weight  basis  but  is  very  difficult  to  cleanse  from 
the  apparatus.  The  Uranine-B  dye  was  readily  adaptable  to  labora- 
tory use  and  gave  reproducible  readings  to  about  1 part  in  10^ 
on  a weight  basis.  This  was  sufficient  for  the  experiment  to  be 
performed  since  the  dye  concentration  in  the  solution  to  be 
ejected  could  be  adjusted  to  produce  samples  along  the  body 
sufficiently  above  the  background  levels  to  give  reproducible 
results.  Calibration  curves  were  experimentally  developed  for 
the  instrument  using  known  concentration  of  Uranine-B  dye  in  the 
tap  water  to  be  used  in  the  experiment.  The  Model  111  Fluoro- 
meter  has  four  sensitivity  settings  allowing  calibration  curves 
to  be  drawn  over  a range  of  about  1 part  per  10^  to  1 part  per 
10  , on  a weight  basis.  For  concentrations  greater  than  the 
lower  value,  dilution  prior  to  measurement  was  performed.  The 
calibration  charts  for  this  experiment  are  shown  in  Figures  21 
and  22.  Approximately  a 3.5  to  4.5  ml  sample  is  required  for 
measurement.  This  sample  is  placed  in  a cuvette  (a  small  test 
tube)  and  placed  in  the  instrument  and  the  measurement  taken. 
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FLUID  CONCENTRATION  URANINE-B  DYE 
OPENING  #1  (WPPM  X 10; 

OPENING  #10  (WPPM  X I02  ) 

Calibration  curves  for  model  111  fluorometer  - openings  1 and  10 


0 0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.S  0.9 


FLUID  CONCENTRATION  URANINE-B  DYE 
OPENING  *3  ( WPP  M X 10) 

OPENING  * 30  (WPPM  X 10  2) 


Figure  22,  Calibration  curves  for  model  111  fluorometer  - openings  3 and  30 
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C.u  i.  i i:. . ' t . 1 1-  ii  1 1-  w i j.i  Hi.'  i mi  i ■ . I 1 1 1. 1 ( i|.  in  . iii  . 1 1 '.m  prior 
to  liii-asui  .’in-  : 1 1 .-.Ini' i I I in>«  i p t 1 nt  ■;  do  effe.t  t ’ i ■ rendli.gs.  Wiping 

Wltli  .1  p.-.pi-r  l.’Wi-l  wo:  I no.mJ  to  111'  suili;  Irill  . All  S.lll.p  I os  WITH 
allowed  to  a.  Iiii-vo  conn  t ( mpc rat  nr v p*  lor  to  t.-idio;*. 

Model  (la  1 i bra r i on 

done t nt  r.i  t j on  IL-asnt  t-inc-nl  Cal  i ! > t it  ion 

Tin  ojivi  in;  ircdel,  a-,  design  .1  o.a't  ijn.  d tom  sample 
co  1 1 it  t i an  . ■ 1 ’ 1 1 ■ 1 1 . • i 1 d i aw  l op,  s imp  1 1 ..  I. h t mgh  l.n.r  1 I nsl . will  or 
tour  p i 1 1 t |ii  iilii'.,  pi  r t ample  r Ii.  mil  .01  I lies' • 1 > 1 loot  ion  < I iambi  ■ r s 
are  all  ionne<  ted  by  3 / 3 2 inch  Inside  diameter  tubing  to  the 
s uc  L i oi  piston.  The  transfer  tube.;  are  ot  unequal  lengths  result  inp, 

In  an  intetc-t  I r:>;  t.  r a ns  1 hi  t analysis  proMein  oi  deilnirig  t lie  suction 
rate  dlrlercti'  cs  p.t  stimple  chamber.  'Hie  prol  lent  was  circumvented 
by  calibrating  the  model  through  a series  of  experiments'.  Addi- 
tionally, a i ,t  1 i br.n  ion  was  desirable  since  it  was  difficult  to  deter- 
mine tin  ex  it  t v a id  /nliu.'.i-  per  sample  ciiarnh  r previously  stated  at 
l.a'Ut'  < e.l  ic  in  be. . which  rc  suited  in  a dilution  factor  oi'  .1  327. 

This  resulted  ir.  a i imp  It-  concentration  lower  titan  the  boundary 
layer  sample  rcmri'm  ration.  It  was  lirst  thought  to  leave  these 
void  spaces  air  tilled  eliminating  the  dilution  factor.  Early 
tests  showed  tost  wit.i  the  body  in  the  tank  in  rite  vertical 
position,  td,  v.iii-i  pressure  head  difference  would  result  in  a 
ventii.g  oi  (he  air  la  these  chambers,  since  the  / were  intercon- 
nected, oni  the  uppermost  chamber  and  a filling  with  tank  water 
through  die  lower  cli  lmbers.  I'refilling  with  water,  as  will  be. 
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discussed  later,  eliminated  the  problem.  To  calibrate  the  model, 
an  assembly  and  calibration  tank  was  constructed.  The  tank 
is  displayed  in  Figure  8 in  the  foreground  having  dimensions  of 
8-inch  diameter  by  3 feet  height.  The  model  is  installed  in  the 
tank  tail  first  as  shown  in  Figure  23  and  fired  by  pulling  on  a 
preattached  lanyard.  The  following  procedure  was  followed  in 
calibrati, .g  the  model.  The  model  was  assembled  excepting  for  the 
half  body  nose  and  placed  in  the  tank  after  connecting  the  lanyard 
to  the  firing  rod.  Tap  water  was  placed  in  the  tank  just  below  the 
level  of  the  polymer  ejection  volume  completely  filling  all  sample 
chamber  volumes.  The  nose  is  then  installed  and  a sample  back- 
ground reading  of  the  water  taken.  The  water  was  then  contami- 
nated with  the  Uranine-B  tracer  to  one  of  three  concentration  levels 
called  the  Ocean  Concentration.  These  were  approximately  .3  x 1C)~7, 
.8  x 10"7  and  2 x 10-7  WPPM,  which  were  determined  to  result  in 
an  appropriate  spread  of  chamber  concentration  levels.  The  model 
was  ther  fired.  The  samples  were  withdrawn  and  measured  in  the 
Fluorometer.  The  samples  were  withdrawn  by  removing  the  model 
from  the  calibration  tank,  wiping  the  outside  dry  while  holding 
in  a horizontal  position  with  the  top,  as  defined  by  the  position 
of  the  internal  suction  transfer  tubes,  being  in  the  uppermost 
position  and  placing  the  model  over  individual  collection  trays. 

The  samples  would  drain  into  these  four  separate  trays  for  analysis. 
It  is  necessary  to  keep  the  suction  transfer  tubes  in  the  uppermost 
position  to  insure  no  transfer  of  liquid  between  the  chambers.  The 
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procedure  was  similar  for  the  probe  case  excepting  that  an  upper- 
most row  of  probes  was  removed  and  then  a bottom  row  to  facilitate  the 
liquid  withdrawal.  A minimum  of  four  calibration  tests  were  run 
at  each  concentration  with  the  flush  probes  and  also  with  the.  pitot 
probes.  A least  squares  fit  was  performed  on  the  data  which  are 
plotted  for  each  sample  chamber  in  Figures  24  through  31  for  high 
and  low  chamber  concentrations.  Deviations  of  +5%  were  noted  in  the 
data.  Several  tests  were  performed  with  a 50  WPPM  polymer  solution 
in  the  tank  with  no  significant  change  in  the  data. 

Stroking  Time  Verification 

Verification  of  the  stroking  time  was  necessary  to 
insure  that  the  appropriate  volume  flow  rates  were  being  achieved. 

The  facility  shown  schematically  in  Figure  32  was  used  for  this 
purpose.  To  calibrate,  the  model  was  loaded  with  liquid  and 
clamped  to  the  bench.  A calibration  spacer  was  placed  on  the 
ejector  follower  and  the  rod  clamped  with  a nylon  shear  screw  to 
the  facility  firing  arm.  When  fired  by  moving  the  firing  arm  in 
a counterclockwise  direction  and  shearing  the  shear  screw,  the 
follower  rod  would  move  Into  the  model  following  the  ejector 
pistcn.  The  calibration  spacer  would  interrupt  the  laser  beam 
for  the  length  of  the  spacer,  starting  the  counters,  discussed 
previously,  giving  a time  over  the  known  distance  and  the 
velocity  of  the  rod.  It  was  found,  by  test,  that  the  velocity 
of  the  rod  was  approximately  5 times  faster  than  the  stroking 
time,  therefore  not  influencing  the  results.  A .5-inch  long 
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Figure  24.  Chamber  1 concentration  vs.  ocean  concentration 

-high  range 


OCEAN  CONCENTRATION  IWPPM  X 10  ) 

Chamber  1 concentration  vs.  ocean  concentration  -low  range 
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0.1  0.2  0.3  0,4  0.3  0.6  0.7  0.6  0.9 

OCEAN  CONCE  NTR  AT  ION  (WPPM  X 10) 

Figure  29.  Chamber  3 concentration  vs  ocean  concentration  - low  range 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 

OCEAN  CONCENTRATION  (WPPM  X 10) 

Figure  31.  Chamber  4 concentration  vs  ocean  concentration  -low  range 


calibration  spacer  was  used  in  the  test  and  placed  at  several  posi 


tions  on  the  follower  rod  to  detect  any  acceleration  of 
None  was  evident.  Nineteen  calibration  runs  were  made, 
rate  resulting  was  12.1  in/sec  +6%  as  compared  to  the  c 
value  of  11.7  in/sec. 


the  rod. 
The  strok 
a 1 culated 
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IV.  ANA!  Y I'  1 CAI.  CONS  I DERATIONS 
Boundary  Layer  Model 

This  .section  presents  the  analytical  np|)roach  for  the  analysis 
of  developing  turbulent  flow  over  an  axisymmetric  body  for  tno 
cases  of  no  polymer,  polymer  ocean  and  polymer  Injection.  A 
model  is  developed  for  predicting  growth  of  the  boundary  layer, 
diffusion  of  tracers  within  the  boundary  layer  and  prediction 
of  skin  friction  coefficients  and  total  drag. 

From  the  literature  review,  the  features  of  previous  studies 
for  turbulent  flows  can  be  summarized  as  follows: 

1.  A modified  law  of  the  wall  type  velocity  profile 
with  pressure  gradient  terms  was  developed  for  the 
thick  axisymmetric  boundary  layer  case. 

2.  Inner  variable  integral  equations  were  developed  for 
the  thick  axisymmetric  boundary  layer  and  verified 
for  certain  parameters. 

3.  Porous  wall  pipe  flow  studies  with  injection  of 
polvmers  have  suggested  significantly  reduced 
diffusion  characteristics.  This  representing  the 
ideal  case  of  introducing  polymers  into  the 
viscous  sublayer. 

4.  External  flow  studies  demonstrate  the  importance 
of  the  ejection  process  on  polymer  efficiency. 
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r> . Far  downs  I tram  dilution  laws  for  |><>  I vinn'  In  I low  those 
for  water  diffusion. 

6.  Similarity  concentration  profile  relations  have  been 
established  for  the  downstream  case. 

7.  The  extent  of  the  initial  mixing  zone  is  considerably 
extended  by  nearly  a factor  of  20  with  polymer  flows. 
Most  investigators  have  obtained  data  in  far  down- 
stream conditions  only. 

8.  Data  on  concentration  profiles  on  axisymmetric  bodies 
are  desirable  to  extend  and/or  verify  dilution  laws 
presently  available. 

In  the  analysis  that  follows,  the  velocity  profile,  continuity 
equation,  momentum  equation,  a Lagrangian  similarity  hypothesis 
for  the  turbulent  diffusion  boundary  layer  growth,  and  a conserva- 
tion of  polymer  equations  for  the  intermediate  and  final  zones  of 
diffusion  are  combined  to  yield  expressions  for  the  growth  of 
the  turbulent  boundary  layer  with  pol\  er  ejection.  Additionally, 
terms  are  added  to  predict  the  initial  zone  region  where 
molecular  diffusion  predominates  for  an  optional  ejection  pro- 
cess. An  alternate  analytical  model  is  postulated  which  would 
eliminate  the  need  for  the  concept  of  an  initial,  intermediate 
and  final  zone  concept. 
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Velocity  Profile  Relation 


Kquaclon  ( V)  provided  a relation  for  the  law  of  the  wall 
accounting  for  pressure  gradients.  Rewriting  equation  (32) 
accounting  for  the  Meyer  (1966)  correction  for  polymers,  AB, 
yields : 


Y+- 

+ l r° 

U = K 1 + 

n.nos  v 
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(1  - e 


2Y+/r+l1/2 


r+ll 

°) 


dY+  + 


AB 


(87) 


An  alternate  method  for  accounting  of  the  polymer  effect  would  be 
by  adjustment  of  the  mixing  length  constant,  K.  As  noted  by 
Virk  (1971),  a change  in  K by  a factor  of  5 is  possible  and  in 
full  agreement  with  data.  This  would  only  apply  for  the  case,  of 
maximum  drag  reduction  along  the  ultimate  asymptote  line,  however. 
Skin  Friction  Relations 

The  boundary  layer  continuity  and  momentum  equations,  equa- 
tions (33)  and  (34),  for  turbulent  axisymmetric  flow 

J_.  (p  u r)  + A-  (p  v r)  = 0 (33) 

9x  9 / 

and 

p u r (|H)  + „ v r (|H,  = - r + -j— (ri)  (34) 


are  applied  to  the  polymer  flow  case  studied  here. 

The  x derivatives  must  be  handled  by  the  chain  rule,  since  each 
of  the  parameters  (Y4,  a,  r^,  Cw)  in  the  law  of  the  wall  is  a 
function  of  x.  Thus,  we  substitute 
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(88) 
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As  may  be  noted  in  (88),  concentration  derivatives  are  added. 
It  is  assumed  throughout  that  pw  and  are  constant  in  this 
analysis . 

The  resulting  boundary  layer  equation,  after  considerable 
algebraic  manipulation  is: 
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where  G,  H,  I,  J are  given  in  the  appendix. 

The  prime  difference  between  equations  (89)  and  (35),  for 
the  polymer  case,  is  the  last  term  accounting  for  the  change  of 
concentration  with  x. 

Appendix  A presents  a complete  development  of  the  equation 
related  to  this  method. 

Prior  to  solution  of  equations  (87)  and  (89),  it  is 
necessary  to  determine  a means  of  calculating  C^.  Test  (1974) 
applied  techniques  developed  by  Lessmann  (1970) , and  Fabula  and 
Burns  (1970)  to  provide  a solution  for  C^.  This  technique 
satisfied  the  case  of  final  zone  of  mixing  but  did  not  account 
for  the  development  of  a diffusion  boundary  layer  in  a similar 
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fashion  to  the  development  of  a hydrodynamic  boundary  layer. 
Nevertheless,  the  results  of  the  work  are  significant.  With  the 
application  of  the  hypothesis  suggested  by  Batchelor  (1957), 
describing  the  development  of  the  diffusion  boundary  layer  and  data 
by  Poreh  and  Cermak  (1964)  and  Poteh  and  Hsu  (1971) , empirical 
relations  are  developed  which  provide  the  values  necessary  to 
determine  the  friction  coefficients  for  the  intermediate  and 
final  zone. 

For  the  case  of  no  polymer  present,  equations  (87)  and  (89) 
may  be  solved  setting  all  terms  containing  equal  to  zero. 

It  is  only  necessary,  as  pointed  out  in  White  (1972),  to  assume  an 
initial  X0  and  prescribe  the  flow  condition  (Ue,  rQ)  as  a function 
of  x.  Additionally,  for  the  case  of  operation  in  a polymer  ocean, 
that  is  Cy  is  everywhere  constant,  all  derivatives  of  Cy  and  with 
respect  to  C„  are  set  equal  to  zero  and  again  a solution  can  be 
obtained  as  for  the  case  of  Cy  equal  to  zero  but  instead  inserting 
a cunstant  value  of  Cw  into  equation  (87). 

Many  of  the  previous  investigators  have  neglected  the  contri- 
bution of  the  initial  zone  in  the  process  of  diffusion  and  resulted 
in  good  agreement  with  data.  The  suggestion  being  that  the  diffu- 
sing process  with  polymers  is  similar  to  that  with  water.  This  is 
true  in  the  intermediate  and  final  region  of  the  diffusion  process 
where  molecular  diffusion  is  not  the  controlling  factor.  Experimen- 
tal methods  and  apparatus  used  significantly  control  the  results. 

For  the  external  flow  case  studied  here,  ejection  in  a laminar 
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flow  region  prior  to  transition  to  turbulent  flow  may  h ignif i cant ly 
influence  the  process.  Data  by  Walters  and  Wells  (1971)  and  Fruman 
and  Tulin  (1974)  suggest  the  initial  zone  becomes  very  large.  It 
is  hypothesized,  here,  that  by  injection  prior  to  transition,  the 
turbulence  intensity  levels  are  extremely  subdued  resulting  in 
elimination  of  the  high  level  mixing  process  associated  with  the 
turbulent  eddies  penetrating  the  viscous  sublayer.  When  injecting 
into  established  flows,  as  performed  by  many  investigators,  the  damp- 
ening of  the  turbulent  intensity  does  occur,  but  not  before  the 
polymer  is  dispersed  at  a rather  normal  diffusion  rate  consistent 
with  other  fluids  (water).  This  results  in  low  wall  concentration , 
rather  quickly,  and  therefore  low  drag  reduction  efficiency.  The 
sections  that  follow  describe  an  analytical  model  which  includes 
an  initial  zone.  Figure  33  schematically  describes  the  model. 

Considering  first  the  initial  zone  where  molecular  diffusion 
is  presumed  to  predominate,  an  approximation  to  the  concentration 
entering  the  turbulent  transition  zone  is  made.  Writing  a sim- 
plified polymer  conservation  relation: 

QiCi  = (Qt  + q2)  Ct  (90) 

where  Q*-  and  Ct  are  the  flow  rate  and  concentration  at  the  transi- 
tion point , respectively , and  Q2  is  an  initial  condition  correction 
factor.  The  volume  flow  rate  Qt  is  given  by  the  mean  local  laminar 
flow  velocity,  .6Ue,  at  the  transition  position  radius  roC  and  a 
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Figure  33.  Schematic  representation  of  diffusion  model 


laminar  layer  thickness,  6t,  determined  from  equation  (78) 

Qt  = .6Ue  27rr0  6t.  (91) 

The  initial  condition  to  the  laminar  flow  boundary  layer  must  be 
applied  if  the  ejected  volume  of  fluid  exceeds  the  fluid  in  the 
laminar  boundary  layer  at  the  periphery  of  the  ejector.  In  equation 
form 

Qi  “ Qe  = Q2 

where  QG  is  the  volume  flux  through  the  laminar  thickness  calcu- 
lated at  the  ejector  periphery  in  a similar  fashion  to  equation 
(91),  and  Q2  is  the  volume  flux  difference  which  provides  an 
"initial  condition"  and  must  be  added  to  the  boundary  thickness 
at  each  station.  The  assumption  is  made  that  the  polymer  solution 
is  uniformly  mixed  within  this  laminar  layer.  The  local  concen- 
tration at  the  transition  point  is  then  determined  from  (90).  The 
volume  flow  rate  into  the  viscous  sublayer  is  determined  by  simple 
ratio  of  the  thickness  of  the  viscous  sublayer  height  to  the  laminar 
layer  height  at  transition. 

Qv  = (Qt  + Q2)  • (92) 

Approximating  in  this  manner  would  tend  to  overpredict  the  flow 
into  the  viscous  sublayer  while  underpredicting  the  flow  in  the 
region  above  the  sublayer.  The  flow  rate  into  the  viscous  region, 
Qv,  and  the  concentration,  Qt,  are  the  parameters  applied  to  the 
drag  reduction  and  boundary  layer  calculations.  This  approxima- 
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mat  ion  allows  for  an  excess  quantity  of  fluid  being  Introduced  above 
the  viscous  sublayer.  The  portion  of  flow  injected  into  the  inter- 
mediate layer  of  the  turbulent  flow  is  given  by 

(Qt  + Q2>  - Qv  = Q3'  ^) 


This  fluid  will  undergo  a rapid  diffusion  characteristic  of  inter- 
mediate and  far  downstream  diffusion  as  discussed  earlier.  The 
concentration  of  polymer  at  the  transition  point  is  as  given 

_ Qi  Ci 


previously,  Ct 


Qt  + Q2 


Considering  equation  (59),  the  relation  for  this  diffusion 
process  is  given  by 


Q3Ct  = 2 tt  (rQ  + yv)  Ct  / u ~ (1  + — )dy  (94) 

v o yv 

;u 

where,  for  the  assumption  of  uniform  concentration  in  the  sublayer, 
Ct  = Cy  on  the  right  side  of  (94).  For  yv  <<  r0  equation  (94) 
becomes 

OO 

Q3Ct  = 2nr0  Cy  / u (1  + f-)dy.  (95) 

0 ^ 0 

Equation  (95)  may  be  integrated  across  the  boundary  layer  at  each  x 
station  resulting  in  a measure  of  the  concentration  profile  above 
the  viscous  sublayer.  It  is  suggested  that  this  relation,  (95), 
be  solved  as  described  in  subsequent  sections,  as  an  independent 
solution  to  the  two-zone  diffusion  process  with  flow  and  concen- 
tration inputs  of  Q3  and  Ct , respectively,  but  other  boundary  layer 
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and  polymer  wall  concentration  parameters  governed  by  the  basic 
initial  zone  calculations.  This  will  avoid  additional  computer 
programming  complexity.  This  calculation  is  only  required  if  a 
measure  of  boundary  layer  concentration  at  particular  heights  in 
the  boundary  layer  is  desired.  The  extent  of  the  initial  zone  has 
yet  to  be  determined.  Fruraan  and  Tulin  ( 19 7 A)  postulate  a model 
for  prediction  cf  the  extent  of  the  initial  zone  for  flat  plates, 
based  on  an  analogy  to  heat  transfer.  They  show,  for  a near 
optimal  ejection  process,  that  th*>  initial  zone  terminates  for 
values  of  dimensionless  distance  greater  than  8.  The  dimen- 
sionless distance  used  to  correlate  these  data  is: 


]UJ>) 

1.5  * 

1 

\ V;jl ' 

s 

Ci_ 

where  s is  the  slot  width.  For  the  axisymmetric  case  with  nose 
ejection,  an  effective  slot  width  may  be  defined  as  the  ejection 
orifice  area  divided  by  the  ejection  circumference.  For  the 
model  being  applied  in  this  study,  the  effective  slot  width  is 
.07837  inches  and  C-p  is  replaced  with  Ct.  No  correction  will 
be  made  to  this  relation  for  application  to  the  axisymmetric 
case  pending  results  of  the  test  program.  It  should  be  reiterated 
that  the  initial  zone  length  will  normally  be  very  small  for 
dilute  solutions  of  polymer  as  would  occur  if  ejection  were  ini- 
tiated in  a highly  turbulent  region.  For  this  type  of  ejection 
process,  an  intermediate  and  final  zone  model  would  be  quite 
adequate. 
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Upon  leaving  the  initial  zone,  relations  for  the  prediction 
of  wall  concentration  are  required.  Application  of  the  Lagragian 
similarity  hypothesis  and  the  concentration  similarity  profiles 
discussed  earlier  coupled  with  the  axisymmetr ic  boundary  layer  rela- 
tions and  velocity  profile  provide  the  necessary  relations. 

A relation  for  solution  of  in  the  highly  turbulent  flow 
field  may  be  extracted  by  application  of  equations  (57),  (58)  and 
(59) 


QiCi  = 2Trr0  / u (1  + f-)  dy  (97) 

0 S?  ro 

where,  for  the  case  of  polymer  flows  with  an  initial  zone, 
would  be  replaced  with  QvCt.  Changing  to  law  of  wall  variables 
and  dropping  the  overbars  since  we  are  dealing  with  average 
values  throughout  the  turbulent  flow  field  analysis 

. 2Y+/r+ 

QiCj  - 2nr0  « r / u+  ~ e ° dY+.  (98) 

0 

Solving  for 

^ = Qici 

Y+  2Y+/r+  (99) 

e 0 

2nr„  v / u+  dY+ 

0 °w 


Equation  (99)  for  polymer  concentration  is  analogous  to  equation 
(34)  for  momentum.  It  is  necessary  to  define  a relation  for  C 
to  solve  equation  (99)  as  it  was  necessary  to  develop  an 
expression  for  u to  solve  equation  (34).  The  similarity  relations 
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given  by  equations  (57)  and  (58)  provide  the  needed  concentration 
expressions 


= -0.693  (y/6d) 


1.5 


Intermediate  Zone 


and 


2.15  Final  Zone 


(57) 

(58) 


C/Cw  . g-0-693  (y/6d)‘ 

Within  the  final  zone,  for  zero  pressure  gradient  flows,  6d/Yp 
remains  constant  at  0.64.  If  there  is  a pressure  gradient,  data 
by  Mellor  (1966)  indicate  that  the  value  of  6d/Ye  changes  to 
.53  in  this  zone. 

The  expression  for  the  growth  of  the  diffusion  boundary  layer, 


<5d,  equation  (55),  completes  the  analysis 
d6d  b v* 


al 


dx 


(55) 


u(Y) 


It  is  only  necessary  to  determine  the  constant,  a^,  in  equation 
(55).  Applying  equation  (48),  changing  variables  and  inserting  a 
similarity  concentration  profile  with  variable  exponent 


Y 


“ -0.693(y/6d)K3 

6d  f °w  e y/<5d  d(y/<$d) 

o 


7 -0.693(y/6d)K3 

J cw  e 
0 


d(y/<Sd) 


(100) 


which  results  in 

Y = K2  6d  = ax  6d  . (101) 

Numerically  integrating  (100)  for  various  values  of  K3  between 
.5  and  5 yields  values  of  K2  which  may  be  applied  to  (101) 

Figure  34  displays  the  results.  Substituting  (101)  into  (55) 
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VALUE  OF 


0.1  1.0  1 0 

EXPONENT  K 3 

Figure  34.  Value  of  K2  for  exponent  Kg 
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yields  a relation  for  6^ 


K2  111  „ R v* 
2 dx. 


(102) 


•a  U1 

where  u-^  is  the  value  of  u at  K26  and  xs , and  K is  the  Karmen 
constant.  Virk  (1971)  has  shown  that  the  constant,  K,  used  in 
the  law  of  the  wall  varies  with  polymer  addition  resulting  in  an 
ultimate  asymptote  for  the  velocity  profile  in  drag  reduction 
cases.  An  effective  value  of  Karman's  constant  may  be  determined. 
This  new  "constant,"  which  is  a function  of  AB,  will  be  called 
K5.  The  manner  in  which  the  modification  of  K to  K5  is  determined 
relates  to  the  basic  equations  for  the  velocity  profiles  with 
and  without  polymers,  equation  (87)  and  (32),  respectively. 
Equating  these  equations 


Y 

3 e 
K / 

0 


1 - 


a + 

2 r0  a 


- e 


2Y+/r^1  1/2 


+AB  = 


k5  * [x  - ! ro (1  - e 


2Y+/r+l  1/2 


0 


+ - 

o) 


dY+. 


Solving  for  yields 

f 


K5  = 


Y+  1 r + 2Y+/r+1l/2 

Ye  fl  - | r+  (1  - e )l  dY+ 

0 ^ J 


iY+,  r , 2Y+/rJii/4 

k / h-  L1  - f r„  (1  - e >J 


(103) 


(104) 


) dY  + AB 
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The  effect  is  apparent.  K5  is  reduced  dependent  on  the  value  ofAB. 
In  effect,  reduction  of  the  eddy  viscosity  through  the  mixing  length 
is  occurring.  Depending  on  the  value  of  u+  andAB  chosen,  K5 
values  approaching  those  of  Virk's  (1971)  ultimate  asymptote  can 
be  achieved.  The  effect  of  this  change  in  K to  some  smaller  value 
= f(AB)  is  to  lengthen  the  intermediate  mixing  zone. 

Equations  (102),  and  (57)  or  (58)  provide  relation  for  the 
solution  of  C/Cy  which,  in  turn,  may  be  substituted  into  (99) 
yielding  an  expression  for  C„.  Equations  (871  and  (89),  the 
boundary  layer  velocity  profile  and  momentum  equations,  respec- 
tively, may  then  be  solved.  These  relations  apply  for  the  inter- 
mediate and  final  zones  of  the  diffusion  process.  Within  the 
initial  zone,  equations  (90)  and  (92)  define  the  polymer  wall 
concentration  and  flow,  respectively.  The  equations  for  the  inter- 
mediate and  final  zone  are  summarized  below: 


k2 


(102) 


-0.693(Y/6d)  2>i5 

— = e For  xc/Ye  > 60  (58) 

Gw  S 


-0.693(Y/6d)1<5 

— = e For  xP/Y  < 60  (57) 

Cw  " e 
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(99) 


Cw  = 


Qici 

Y+  2Y+/r+ 

2vr0  v fe  u+  e ° d^ 

0 ^ 


1 


K 


Y+ 

o 


f 

0.1108 


2Y+/r+ 

' °)J  1/2 


dY+  + 


AB  (87) 


dX_ 

dx* 


(3  a H - Gj_) 


X Gx 


X2ro  , 2Y+/r+ 
— (e 


1) 


(89) 


It 

H 


RLV 


drT 
o 

dx* 


I - 


X^wj 

dx* 


where  the  terms  are  as  defined  previously. 

An  alternate  approach  may  be  applied  to  attack  the  problem 
for  developing  flow.  It  is  assumed  the  similarity  expression  for 
the  concentration  profile  holds  throughout  the  boundary  layer. 

Lacking  data  for  the  effect  of  turbulence  intensity  changes  in 
developing  flow  with  polymer  ejection,  it  is  hypothesized  that  the 
similarity  relation  is  functionally  related  to  the  boundary  layer 
thickness,  <5,  the  diffusion  boundary  layer  thickness,  6^,  the 
Karman  constant,  K,  and  a distance  from  the  inception  of  turbu- 
lence to  the  ejection  slot,  xe,  (set  equal  to  1 for  cases  of  the 
ejection  being  at  the  transition  to  turbulent  flow  point  or 
before).  The  exponent  in  the  similarity  expression  is  the  controlling 
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parameter  since  It  indicates  the  diffusion  rate.  The  expression 
for  the  similarity  profile  may  be  written  as 


C_  = e -0.693  (Y/6d)K3 
Cw  " P 


where  K.3  = f(6,  6d,  R,  xe)„ 

Poreh  and  Cerinak  (1964)^  as  discussed,  have  empirically  assigned 
specific  zones  relating  downstream  distance  and  the  boundary 
layer  thickness  to  establish  the  exponent  value.  It  seems 
apparent,  that  in  the  intermediate  zone,  the  value  <•  c the 
exponent  should  not  be  constant.  Concentration  profiles  as  a 
function  of  downstream  distance  by  Walters  and  Wells  (1971)  bear 
this  out.  Rewritting  (105)  in  terms  of  K5 

K3  = f (6 , 6d,  K5,  xe).  (106) 

Assuming  a nondimensional  relationship  which  allows  the  exponent 
to  grow  as  the  diffusion  boundary  layer  grows  results  in 

K3  = C K5  r • (107) 


For  final  zone  of  diffusion  where  the  characteristic  diffusion 
rates  of  polymer  and  water  have  been  shown  to  be  the  same, 


K5  = K = .4 


.64  for  zero  pressure  gradient  flows 

1 since  ejection  preceeds  the  turbulent  transition 
point  and  K3  = 2.15, 
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Solving  for  C in  (107)  yields 

C - 8.4  (108) 

From  whx"h  (107)  is  defined  as 

K3  =>  8.4  K5  . (109) 

For  developing  polymer  flows  with  the  ejector  at  or  preceeding  the 
point  of  transition  to  turbulent  flow  (using  the  test  model 
described  earlier) , 

K5  = .2 

6^  = .0008  ft  (for  the  model  flow  rates) 
and  6 = .003  ft  (equation  (80)) 

the  value  of  K3  becomes 
K3  = .45 

and  the  value  of  the  moment,  K2,  may  be  determined  from  Figure  34. 
The  constant  K3  * .45  implies  an  extremely  suppressed  diffusion 
process,  which  in  fact,  on  computation,  occurs.  Application  of 
(105)  with  the  Lagrangian  similarity  hypothesis  for  predicting  the 
diffusion  boundary  layer  growth  would  eliminate  the  need  for  a 
multi-zone.  The  exponent,  K3 , would  be  self-adjusting  throughout 
the  boundary  layer.  The  simple  form  chosen  for  K3  may  be  incorrect. 
Sufficient  data  would  be  required  in  all  diffusing  zones  to  properly 
define  whether  the  similarity  profile  extends  to  the  boundary 
layer  and  to  define  the  proper  functional  form.  Application  of  this 
method  would  require  integration  of  equation  (100)  rather  than 
use  of  Figure  34.  This  is  required  since  equation  (100),  although 
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integrated  between  0 and  03  in  reality  reaches  final  values  for 
values  nearer  4 or  5.  For  cases  where  may  be  small,  as  for 
much  reduced  diffusion,  the  value  may  become  quite  large  and 
should  be  bounded  by  the  viscous  sublayer  thickness,  for  6^  , 

and  the  boundary  layer  thickness  for  6max.  This  will  prevent  the 
moment,  K2 , from  becoming  very  Irrge  for  low  values  of  K3  as  implied 
by  Figure  34.  Extensive  effort  in  pursuing  this  approach  was  not 
expended  since  the  molecular  diffusion  coupled  with  the  two  zone 
model  is  believed  to  represent  a sounder  approach. 

The  equations  may  be  solved  for  the  skin  friction  distribu- 
tion Cf(x),  the  boundary  layer  thickness,  6,  and  the  diffusion 
boundary  layer  thickness,  <5j.  Additionally,  the  concentration  at 
the  wall,  C^,  and  at  any  position  in  the  boundary  layer,  C,  for 
a given  polymer  or  tracer  flow  and  ejection  concentration  and 
ejection  rate  can  be  determined.  The  initial  conditions  required 
are  an  estimation  of  Aq,  and  the  known  flow  condition  Ue  and  r0 
as  functions  of  x.  Inputs  of  velocity,  temperature  and  related 
parameters,  ejection  rate  and  polymer  concentration  are  required. 

It  has  been  assumed  that  the  flow  becomes  turbulent  at  the  minimum 
pressure  point  and  all  initial  conditions  are  determined  at  this 
point.  From  these  input  conditions,  6^  may  be  determined  from 
a mass  balance  at  r0(x)  and  60  may  be  calculated  from  (80). 

AB0  may  be  calculated  from  (42).  An  estimate  for  the  skin 
friction  term  AQ  comes  from  the  flat  plate  relation. 


Cf  = 0.0592  R, 


(110) 


-1/5  = 2_ 

***  X2  ' 

Figures  35,  36  and  37  present  the  pressure  coefficients  vs  nor- 
malized axial  distance  for  the  three  bodies  discussed  earlier.  The 
pressure  coefficients  and  other  flow  conditions,  Ue/U0,  have  been 
determined  using  the  Douglas-Neumann  Potential  Flow  Program  (1958) 
and  are  required  for  solving  the  previous  system  of  equations. 

The  total  drag  on  the  body  has  been  determined  by  integrating  the 
calculated  skin  friction  and  by  integrating  the  pressure  distribu- 
tion over  the  body  to  the  point  of  separation.  As  pointed  out  by 
White  (1972),  equation  (89)  provides  a distinct  test  for  separation 
when  the  coefficient  (3aH-Gi)  vanishes  which  causes  X to  approach 
°“  and  Cf  approaches  0.  At  this  point,  the  pressure  is  determined 
and  multiplied  by  the  projected  area  at  the  radius,  rQ(x)  of 
separation.  Thus,  both  friction  and  form  drag  are  accounted  for. 

The  form  drag  calculated  in  this  manner  may  exhibit  large  errors 
due  to  the  steepness  of  the  calculated  pressure  coefficient  in  the 
tail  region  as  seen  in  Figures  36  through  37.  Data  on  pressure 
coefficients  for  similar  bodies,  Nadolink  (1958),  do  not  exhibit 
this  steepness.  The  calculation  is  simple,  however,  and  does 
provide  an  indication  of  the  total  drag. 

The  equations  have  been  programmed  on  a Control  Data  Corpora- 
tion Computer,  CDC  3300.  All  numerical  integrations  use  the  stan- 
dard RUNGE-KUTTA  subroutine.  Appendix  B contains  a sample  of  the  com- 
puter program  and  solution  for  the  combined  model.  A difficulty  does 
occur  in  the  formulation  when  a favorable  pressure  is  encountered 
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Figure  37.  Calculated  pressure  distribution  over  axisymmetric  body  - spherical  tail  configuration 


(negative  “>) . Equation  (89)  diverges  resulting  in  no  solution. 

The  problem  is  resolved  by  not  allowing  a.  to  become  negative. 

This  does  introduce  a discontinuity  into  the  solution  but  recovery 

from  the  discontinuity  is  rapid.  Two  methods  of  solution  are 

possible  for  the  determination  of  X in  equation  (89).  The  full 

equation  may  be  used,  which  for  small  polymer  wall  concentration 
dC„ 

gradients  , provides  for  stable  solution.  For  higher  concen- 

dx 

tration  gradients,  an  alternate  approach  is  recommended.  A 
solution  for  X may  be  obtained  neglecting  the  last  term  in 
equation  (89).  The  value  of  X is  then  corrected  at  each  calcula- 
tion interval  as  noted  below: 

X ■ X calculated  “ ‘ 4B>  <1U> 

where  AB0  is  the  previous  calculated  value  of  AB.  In  this  manner, 

the  new  calculated  value  of  X Is  corrected  for  an  updated  AB 
resulting  in  the  appropriate  X.  The  full  equation  (89)  does  pro- 
vide a solution  at  the  higher  concentration  gradients  but  the  solu- 
tion tends  to  become  oscillatory  requiring  a large  number  of  calcu- 
lation points  to  minimize  the  oscillations. 

The  analytical  model  has  been  exercised  in  several  different 
ways  to  test  the  results  that  may  be  obtained.  For  purposes  of 
ease  in  discussing  these  models,  they  will  be  given  specific  names. 
First,  the  molecular  diffusion  model  only  contains  the  molecular 
diffusion  portion  of  the  analytical  program  and  never  switches 
to  a higher  diffusion  two-zone  model.  The  second  model  is  called 
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the  two  zone  mode  1 , the  Initial  zone  is  vorv  small  and  not  Ineluded. 
llils  model  Is  applicable  to  water  diffusion  or  to  polymer  diffusion 
If  the  ejector  were  placed  In  an  established  turbulent  flow.  The 
third  model  will  be  railed  the  combined  model  and  represents  the 
near  optimal  external  flow  process  for  polymer  ejection,  ejecting 
into  the  laminar  region  prior  to  turbulent  flow  transition.  The 
postulated  model,  with  variable  K3  is  tested,  in  a single  case, 
for  comparison.  'Die  total  model  is  shown  schematically  in  Figure 
33. 

The  calculations  were  all  made  using  the  6°  tail  model,  a 
velocity  of  27  ft/sec,  an  ejection  rate  of  20.6  in^/sec,  unless 
otherwise  noted,  and  a medium  temperature  of  60°F.  The  length 
Reynolds  number  is  4 . 5 x 10^.  The  primary  emphasis  in  these 
calculations  has  been  on  the  concentration  profiles  and,  more  spe- 
cifically, wall  concentration  profile.  It  should  be  remembered, 
however,  that  the  computer  routine  and  analytical  model  allow 
determination  of  concentration  throughout  the  boundary  layer, 
the  growth  of  the  diffusion  boundary  layer  within  the  hydrodynamic 
boundary  layer  and  other  boundary  layer  parameters.  Skin  friction 
and  total  drag  is  also  determined  and  a comparison  of  the  drag 
reduction  resulting  from  the  different  diffusion  processes  as 
the  three  models  (really  one  applied  differently  excepting  the 
variable  K3  model)  predict  are  compared  in  Table  4. 


Figure  38  displays  polymer  wall  concentration  ratios  vs  axial 
distance  calculated  using  the  molecular  diffusion  model  and  the 
two-zone  model.  For  comparison  purposes  in  the  drag  reduction 
calculations,  the  water  case  for  the  two  zone  model  is  considered 
to  be  the  baseline.  This  model,  for  water,  is  accurate  as  will 
be  shown  in  subsequent  sections.  The  molecular  diffusion  model 
presents  the  optimal  case  for  drag  reduction  on  external  flows. 

The  viscous  sublayer  growth  is  rather  small  over  the  body  length, 
a growth  of  less  than  .001  inches  over  the  x/L  distance  displayed, 
resu lting  in  the  flat  concentration  profile.  The  prime  reason  for 
the  reduced  concentration  to  that  ejected  is  due  to  the  diameter 
ratio  from  the  ejection  diameter  to  the  near  maximum  body  diameter. 
Since,  for  molecular  diffusion,  the  concentration  at  the  wall  does 
not  change  appreciably,  the  single  line  shown  applies  to  all 
polymer  concentration  so  long  as  Y equals  11.5.  For  values 
of  polymer  concentration  less  than  50  WPPM,  AB  in  the  model  would 
change  resulting  in  a wall  concentration  variation.  For  this 
case,  only  v*  affects  the  calculations.  The  two  zone  model  for 
50  and  500  WPPM  shows  rapid  dispersion  of  the  polymer  resulting  in 
less  than  1/10  the  ejected  concentration  at  about  1/2  the  body 
length.  This  type  of  diffusion  process  would  predict  high 
polymer  concentration  or  ejection  rates  required  for  reasonable 
drag  reduction.  The  model  would  be  reasonably  accurate  for  long 
bodies  where  the  initial  zone  is  small  compared  to  the  body  length. 
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I'lio  elicit  of  concentration  !h  evident  with  a maximum  wall  concen- 
tration, over  the  portion  of  the  body  shown  exceeding,  that  required 
for  maximum  drag  reduction  for  ejected  concentration  above  approxi- 
mately 200  WPPM.  The  model,  as  Implemented,  does  not  include 
changes  in  viscosity  for  high  polymer  concentrations.  The  effects 
must  be  accounted  for  if  large  concentrations  are  applied. 

Figure  39  displays  the  two-zone  model  compared  to  a data 
adjusfed  model  based  on  data  generated  in  the  test  program.  Essen- 
tially, the  two-zone  model  Is  reduced  to  a single  zone  model  with 
the  exponent  value  in  tlie  similarity  expression  (58)  being  that 
determined  from  the  experimental  data,  .75.  The  effect  of  suppressed 
diffusion  is  apparent  with  concentration  levels  being  much  higher. 
Care  must  be  taken  in  interpreting  this  figure,  as  will  be  discussed 
in  the  experimental  section  where  the  data  is  qualified.  For  the 
moment,  the  effect  of  a change  in  the  exponent  is  the  important 
factor  shown  here. 

Figure  40  displays  the  variable  K3  model  predicting  very  low 
diffusion  and  therefore  high  wall  concentrations’ as  compared  to 
the  two-zone  model.  The  results  are  interesting  warranting  further 
development  of  the  functional  form  of  the  K3  function  (106). 

Figure  41  displays  the  results  from  the  combined  model  com- 
pared to  the  two-zone  model.  The  transition  from  the  molecular 
model  to  the  two-zone  model  is  clearly  evident  as  controlled  by 
the  effect  of  concentration  on  the  initial  zone  length  in  equation 
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(96)  and  data  from  Fruman  and  Tulin  (1974).  This  model  should 
adequately  predict  the  polymer  diffusion  process  for  the  initial, 
intermedia i_e ,•  and  final  zone  of  diffusion  subject  only  to 
verification  of  the  abrupt  transition  from  the  initial  to  inter- 
mediate zone  similarity  profile.  Quite  possibly  a more  gradual 
transition,  as  suggested  by  the  variable  K3  model,  would  be  appro- 
priate. Figure  42  displays  the  effect  of  the  ejectior  rate,  for 
a fixed  50  WPPM  polymer  concentration,  on  the  wall  concentration 
ratio.  The  shift  in  the  transition  from  the  initial  zone  to  the 
intermediate  zone  is  apparent  and  governed  by  the  changed  ejection 
velocity  in  equation  (96).  Figure  43  taken  from  Fruman  and  Tulin 
(1974)  displays  wall  concentration  data  calculated  from  the  com- 
bined model  plotted  against  the  reduced  dimensionless  distance, 
equation  (96).  As  may  be  seen,  the  combined  model,  applying  the 
Lagrangian  similarity  hypothesis,  results  in  good  agreement  with 
data  obtained  by  Tulin. 

The  percent  total  drag  reduction  and  percent  skin  friction 
reduction  have  been  calculated  for  the  cases  discussed  using  the 
relations 

% TDR  = 100  (1  - -ft—--)  (112) 

uw 

and 

% SFR  = 100  (1  - ~---LY)  (113) 

SDW 
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NORMALIZED  AXIAL  LENGTH  X /L 

Figure  42.  Effect  of  polymer  ejection  rate  on  wall  concentration  ratio  - 
combined  model  - 5G  WPPM  polymer  concentration 
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Figure  43.  Dimensionless  concentration  vs  reduced  dimensionless  distance 
(from  Fruman  and  Tulin  (1974)  ) 


where 


DpQLY  = total  calculated  drag  with  polymer  ejection 
Dw  = total  calculated  drag  with  water  ejection 
SDpQLY  = shear  drag  calculated  with  polymer  ejection 
SDty  = shear  drag  calculated  with  water  ejection. 

Table  4 lists  these  results.  A discussion  of  the  table  will  not 
be  given  since  the  results  are  as  would  be  expected  for  the  wall 
concentration  shown  in  the  previous  figures. 

The  system  of  equations  presented  herein  is  believed  to  have 
several  unique  qualities.  First,  presuming  that  a reasonable 
ejection  process  occurs,  the  equations  account  for  a change  in  the 
mixing  length  constant  resulting  in  a significantly  reduced  diffu- 
tion  process.  This  should  represent  a near-optimum  case  for  develop- 
ing flow.  Secondly,  by  maintaining  the  mixing  length  constant  equal 
to  that  of  water,  the  model  represents  the  ejection  of  polymers  in 
a fully  developed  turbulent  flow  where  the  similarity  concentration 
profiles  are  equivalent  to  those  of  water.  Thirdly,  the  model 
accounts  for  the  initial  zone  of  diffusion  as  well  as  the  subsequent 
zones,  which  should  provide  an  accurate  prediction  of  the  wall  con- 
centration. These  beliefs  will  be  tested  in  subsequent  sections. 
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TABLE  4 


Calculated  Draft  Reductions 


Model  Applied 

Conditions 

Pol v, nor  CONC/Q1 (inJ/secl 

% TDR 

% SFR 

2 zone 

Water/ 20. 6 

0 

0 

Molecular 

Diffusion 

50,500/20.6 

32.6 

65 

Variable  K.3 

50/20.6 

31.2 

62.3 

Data  Adjusted 

50/20.6 

500/20.6 

27 

32 

55 

65 

2 zone 

50/20.6 

500/20.6 

22.6 

32.6 

44.6 

65 

Combined 

50/41 . 2 
50/20.6 
50/10.3 
100/20.6 
500/20.6 

30.3 
24.6 
17.8 
30.1 

32.3 

60.6 

50.3 
37.8 

60.4 
64 . 2 
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Two  typos  of  tests  were  performed  in  this  research.  These  were: 

(1)  polymer  ocean,  drag  reduction  and  flow  visualization  tests  and 

(2)  polymer  ejection  boundary  layer  sampling  tests.  Both  series 
made  use  of  the  full  experimental  apparatus  but  in  slightly  different 
ways.  The  experimental  procedure  for  the  polymer  ocean  tests  will 

be  discussed  first. 


A . Polymer  Ocean  Tests 

General 

The.  polymer  ocean  tests  made  use  of  the  three  dye  ejecting 
bodies  described  previously  and  displayed  in  Figures  11,  12  and  13. 
The  test  series  is  called  polymer  ocean  since  for  these  tests  the 
entire  drop  tank  is  filled  witli  either  fresh  water  or  filled  with 
a polymer  solution.  The  model,  therefore,  is  operating  under 
conditions  of  uniform  polymer  concentration.  A visible  dye  was 
ejected,  by  aspiration  in  many  of  the  tests. 

Polymer  Addition 

The  polymer  used  in  these  tests  was  polyethylene  oxide 
(Polyox-WSR-301 , Union  Carbide  Corp.).  Polyox  is  a long,  3 million 
molecular  weight,  water  soluble  polyether.  Due  to  its  considerable 
size  and  the  sensitive  ether  linkage,  the  mixing  of  highly  concen- 
trated solutions  must  be  undertaken  with  great  care.  A 25%  by  weight 
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polymer  slurry  was  made  up  for  mixture  with  the  water  in  the  tank. 
'ITie  test  program  called  for  polymer  ocean  concentrations  of  fiO  Wl’PM 
to  1.25  Wl’PM.  Since  the  water  tank  weight  of  water  equaled  3452. 5 
lbs,  3452.5  x ICT6  lbs  of  polymer  addition  were  required  for  each 
WPPM  of  polymer  concentration  desired  (or  13810  x 10-(l  -bs  of  25% 
polymer  slurry  by  weight).  The  tank  was  first  brought  to  60  WPPM 
concentration  level,  tests  performed,  and  then  drained  to  a pre- 
determined level  and  refilled  with  fresh  water  to  obtain  50  WPPM. 
After  testing  at  this  level,  the  procedure  was  repeated  to  obtain 
20  Wl’PM  and  5 WPPM.  The  tests  requiring  concentrations  down  to 
1.25  WPPM  and  2.5  WPPM  were  begun  with  a 20  WPPM  concentration. 

The  mixing  procedure  consisted  of  premixing  the  preweighed 
25%  polymer  slurry  in  a 25-gallon  drum  by  slowly  pouring  the 
slurry  against  a wall  of  the  drum,  and  spraying  a high  intensity 
fan  shape  spray  of  water  on  the  polymer.  Gentle  stirring  followed 
for  approximately  1/2  hour  and  the  mixture  left  to  stand  for 
2 to  4 hours.  The  clear  solution  was  then  poured  into  the  par- 
tially empty  drop  tank  and  the  tank  refilled.  The  tank  was  gently 
mixed  by  raising  and  lowering  the  retrieval  net  to  insure  a 
homogeneous  mixture.  Samples  were  taken  at  each  concentration  and 
measured  in  a Hydrat  apparatus  similar  to  that  described  in  Hoyt 
(1966). 
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Model  Preparation 

The  models  for  testing  were  configured  in  the  6°  tail,  12° 
tail  or  spherical  tail  configuration  as  discussed  previously.  The 
model  preparation  for  test  was  as  described  previously.  It  is 
repeated  here  for  completeness.  Referring  to  Figures  14  and  15, 
the  model  is  made  ready  for  a launch  by  inserting  the  tank  guide 


(j^)»  cl*e  w 


ire  bears  on  a nylon 


wire  through  the  body  guide  tube, 
bearing  at  either  end  of  the  model  to  minimize  friction.  The 
tail,  (io),  of  the  model  is  removed  (separated  sufficiently)  and 


the  internal  cavity  filled  with  dye.  The  dye  ejection  ports 
being  previously  taped  over  to  insure  no  loss  of  dye.  The  model 
is  reassembled  and  the  vent  screw  (unmarked)  in  the  tail  is 
removed  and  the  model  placed  under  the  water  level  in  the  tank 
to  back  fill  the  cavity  totally  with  water.  The  model  may  now 
be  fired  upon  removal  of  the  tape  over  the  ejection  holes.  When 


fired,  stagnation  pressure  enters  the  4 stagnation  ports 


,0 


is  transmitted  through  the  body  guide  tube  and  pressurizes  the 


,0 


bladder,  ( 4 j,  forcing  the  dye  cut  the  four  dye  ejection  ports 
located  at  the  minimum  pressure  point  of  the  forebody  as  deter- 
mined by  a potential  flow  program.  After  several  tests,  the 
bladder  was  found  to  be  not  necessary  simplifying  the  operation. 
The  tail  of  the  model  was  then  inserted  into  the  launcher  in 
preparation  for  launching.  The  launcher  receptacle,  as  designed, 
mates  with  the  6°  tail  only.  The  other  models  are  butted  to  the 
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launcher.  The  models  are  held  in  place  by  flshllne  that  is  tied  to 
one  side  of  the  launcher,  wrapped  around  the  wire  In  front  of  the 
nose  of  the  model  and  tied  to  the  other  side  of  the  launcher.  At 
launch,  the  string  is  severed  by  the  forward  motion.  When  in  the 
ready  to  fire  condition,  all  except  the  tail  is  underwater.  After 
firing,  the  model  is  maintained  on  the  wire  in  the  tank,  the  ejection 
holes  retaped  underwater,  the  tail  removed  and  the  internal  chamber 
drained  and  refilled  with  dye.  The  model  is  reassembled  as  before 
and  made  ready  for  a launch. 

Launch  Procedures 

The  launcher  operation  was  quite  simple.  The  low  pressure 
regulator  on  the  launcher  was  preset  to  some  known  pressure  level. 

The  hand  firing  valve  was  opened, f iri ng  the  launcher.  The  launcher 
chamber  pressure  was  recorded  during  this  operation.  The  model 
velocity  was  measured  at  the  several  stations  and  noted  for 
acceleration  or  deacceleration.  The  launcher  pressure  was  then 
adjusted  to  obtain  steady  state  velocities  during  the  model  travel 
down  the  tank.  After  fire,  the  firing  valve  is  closed,  the 
launcher  piston  chamber  vented,  and  the  firing  piston  retracted 
in  preparation  for  the  next  test. 

Velocity  and  Photographic  Measurements 

Little  needs  to  be  said  regarding  procedures  here.  It  was 
only  necessary  to  reset  the  counter  and  set  the  camera  flash  guns 
delay  circuitry  to  the  appropriate  value.  At  launch,  the  camera 
shutter  was  held  open.  The  appropriately  delayed  flash  would 
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expose  the  film.  Transit  times  through  the  various  stations  were 
automatically  recorded. 

Polymer  Ocean  Test  Series 

Two  types  of  tests  were  performed.  Launch  tests  in  which  the 
launcher  was  used  and  free  drop  tests,  in  which  the  model  was  held 
just  underwater  and  allowed  to  gravity  drop  down  the  tank.  The 
launch  tests  were  all  performed  with  the  6°  tail  model.  The  free 
drop  tests  were  performed  with  the  12°  tail  and  spherical  tail 
model.  This  type  of  test  was  initiated  since  gross  instabilities 
in  trajectory  of  the  models  (12°  and  spherical  tail)  occurred  when 
launchings  were  attempted.  Undoubtedly  due  to  the  launcher 
receptacle  having  a 6°  taper  design  to  fit  that  model. 

The  test  series  are  summarized  in  Table  5.  For  the  launch  tests, 
only  those  tests  at  which  terminal  velocity  was  achieved  are  noted.  A 
total  of  491  tests  were  run  of  which  135  tests  are  reported. 
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TABLE  5 


Summary  of  Tests  Performed 
Polymer  Ocean  Series-Nu iber  of  Tests 


POLYMER  CONCENTRATION 
PPMW 


FRESH  WATER 
1.25 
2.5 
5 

10 

20 

50 

60 


LAUNCH  TESTS 
6o  TAIL 


4 

10 

7 

2 

12 


FREE  DROP  TESTS 
12°  TAIL  SPHERICAL 
TAIL 


7 

8 

6 

9 

6 

6 

6 


7 


11 


1.9 


B . Polymer  Ejection  Tests 

General 

The  polymer  ejection  tests  made  use  of  the  polymer  ejecting 
body  described  previously  and  displayed  in  Figure  18.  In  these 
tests,  varied  polymer  concentrations  as  well  as  pure  water  all 
dyed  with  a fluorescent  dye,  Uranine-B,  were  employed.  Simul- 
taneously, either  wall  samples  or  samples  at  two  differem.  dis- 
tances from  the  wall  in  the  boundary  layer  were  taken  at  four  axial 
positions.  Visible  dye  was  also  mixed  with  the  polymer  on  specific 
tests  for  boundary  layer  photographs. 
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The  polymer  used  in  this  series  of  tests  was  also  polyethylene 
oxide  (rolyox-WSR-301 ) . 

The  test  matrix  called  for  solutions  in  the  50  to  1000  PPMW 
range.  A master  solution  was  prepared  at  2000  ppmw  and  diluted  to 
the  required  concentration. 

The  master  solution  was  prepared  by  sifting  the  premeasured 
polyox  powder  onto  the  surface  of  the  carefully  weighed  waLer  which 
was  being  slowly  stirred  by  a magnetic  mixing  bar.  After  a clear 
solution  resulted  (approximately  1 hr.),  the  solutions  were  left 
to  stand  for  approximately  40  hrs.  to  minimize  the  high  viscoelastic 
effects  of  freshly  mixed  solutions,  and  to  assure  homogeneity. 

The  master  solution  was  kept  in  a dark,  cool  place  after  initial 
mixing  to  reduce  the  auto-oxidation  problem  with  polyethers.  When 
required,  the  proper  amount  was  pipetted  out  of  the  master  solution 
and  diluted  to  the  required  concentration.  These  dilute  solutions 
were  used  within  24  hrs.  of  the  make-up  time. 

The  uranine  master  solution  was  prepared  in  the  same  manner. 
Thus,  when  a test  solution  was  required,  the  proper  amount  of 
each  constituent  was  added  to  enough  water  to  make  one  liter  of 
solution.  This  was  poured  back  and  forth  into  the  container  ten 
times  to  assure  complete  mixing,  and  then  let  stand  for  one  hour 


to  achieve  homegeneity. 


A series  of  calibration  tests  with  the  fluorometer  indicate 


that  this  procedure  is  adequate  to  obtain  accurate  results. 

The  phenolphthalein  solution  was  also  mixed  with  the  polymer 
solution  in  the  same  manner  as  the  uranine  to  provide  a visible 
dye  for  photographic  record  purposes.  Fluorometric  readings  were 
not  taken  during  these  tests  since  the  slight  opaqueness  of  the 
phenolphthalein  resulted  in  erroneous  fluoressence  readings  were 
attempted . 

The  test  series  conducted  required  polymer  solutions  of  5 PPMW, 
10  PPMW,  20  PPMW,  50  PPMW,  500  PPMW  and  1000  PPMW. 

Model  Preparation  and  Data  Retrieval 

The  ejecting  6°  tail  model  has  been  previously  described.  The 
assembly  for  a test  and  the  subsequent  withdrawal  of  samples  is 
discussed  below.  The  model  is  assembled,  as  described  previously, 
excepting  for  the  half-body  nose  and  placed  in  the  assembly  tank, 
Figure  23.  The  appropriate  flush-mounted  probes  or  boundary  layers 
pitot  tubes  are  installed.  The  installation  procedure  for  the  flush 
probes  consists  of  partially  inserting  the  flush  probes,  placing 
a small  amount  of  wax  around  their  periphery,  heating  lightly 
with  a propane  torch  and  pressing  the  probes  flush  with  the 
surface.  The  pitot-like  probes  were  installed  in  a similar  fashion 
excepting  that  a feeler  gauge  set  at  either  .020  inches  or  .050 
inches,  depending  on  probe  set,  was  placed  under  the  flattened  tip  of 
the  probe.  Figure  20.  The  probe  area  at  the  surface  of  the  body 
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was  filled  with  wax  and  f laired  in  with  a razor  blade.  The  height 


could  be  reasonably  set  within  .003  inches  by  this  method.  The 
axial  orientation  was  set  by  eye,  but  with  the  extent  of  the  probes, 
.03  inches,  an  ofi.et  of  several  degrees  was  easily  noticed  and 
corrected. 

The  assembly  tank  was  now  filled  with  water  to  just  below 
the  polymer  storage  volume  level.  A background  sample  of  the 
assembly  tank  water  and  drop  tank  water  was  taken  at  this  time 
to  correct  sample  fluorence  readings.  The  solution  to  be  ejected 
was  then  poured  carefully  into  the  polymer  storage  volume,  Figure 
17,  and  the  halfbody  nose  carefully  installed.  'The  model  now 
resides  in  the  assembly  tank,  nose  up,  and  must  be  transferred  to 
the  drop  tank,  nose  down.  The  ejector  screen  was  then  taped  with 
electrical  tape  to  insure  no  leakage  occurred  on  transfer  to  the 
drop  tank.  The  guide  wire  was  threaded  through  the  model  and  the 
model  transferred  to  the  drop  tank.  The  model  is  now  inserted 
into  the  launcher  receptacle  and  a nylon  shear  screw  affixed  through 
the  launcher  to  the  firing  rod  on  the  model.  This  holds  the  model 
in  place  and  fires  the  CC>2  cartridge  ejecting  the  polymer  and  with- 
drawing boundary  samples  at  launch.  After  firing  and  the  launcher 
retracted,  the  model  is  kept  underwater  while  preparations  are 
made  to  quickly  wipe  it  dry  on  removal  from  the  tank  prior  to 
placing  on  four  separate  collection  trays.  The  model  is  removed 
f>-om  the  wire,  wiped  dry  and  placed  on  a stand,  horizontally,  over 


four  collection  trays,  (keeping  the  top  uppermost  for  reasons 
explained  under  the  calibration  section).  After  the  samples  have 
drained  into  the  collection  trays,  the  model  was  totally  dis- 
assembled and  washed  prior  to  the  next  test.  The  chamber  samples 
were  then  tested  applying  the  fluorometric  techniques  described 
previously.  The  fluorometer  readings  were  transferred  to  fluorescent 
dye  concentration  measurements  in  PPMW  using  the  appropriate  cali- 
bration curve,  Figure  21  or  22,  and  then  converted  to  undiluted 
values  applying  the  external  to  internal  calibration  curves, 

Figures  24  through  31.  These  values,  when  nondimensionalized  with 
the  ejected  dye  concentration,  represent  the  normalized  concen- 
tration of  any  ejected  tracer,  or  the  polymer,  at:  that  particular 
station  or  boundary  layer  height  under  the  assumption  of  similar 
diffusion  qualities. 

Launch  Procedure 

The  launch  procedure  is  the  same  as  for  the  previous  tests. 

No  attempt  was  made  to  achieve  terminal  velocity  since  a signifi- 
cant number  of  tests  are  required  to  do  this.  The  ejector  type 
tests  took  approximately  one  hour  per  complete  test  cycle  as  com- 
pared to  13  minutes  for  the  polymer  ocean  tests  (when  photographic 
records  were  desired). 

Velocity  and  Photographic  Records 

The  procedures  for  these  tests  were  as  described  for  the 
polymer  ocean  tests. 
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Polymer  Ejection  Test  Series 


The  polymer  ejection  tests  were  performed  with  fresh  water 
and  six  polymer  concentrations.  As  previously  mentioned,  no 
attempt  was  made  to  control  velocity.  The  test  series  is  summa- 
rized in  Table  6.  The  ejection  rate  for  all  tests  was  held  con- 
stant at  20.625  in-Vsec. 


TABLU  6 

Summary  of  Tests  Performed 
Fnlymer  Ejection  Series-Number  of  Tests 


POLYMER 

CONCENTRATION 

PPMW 

FLUSH 

PROBE  HEIGHTS 
.020  IN. 

.050  IN. 

PICTURE 

TESTS 

FRESH  WATER 

4 

3 

7 

5 

50 

4 

4 

4 

1 

500 

6 

4 

3 

1 

1000 

3 

- 

- 

2 

20 

1 

- 

- 

- 

10 

2 

- 

- 

- 

5 

2 

— 

_ 
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VI.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


A . General 

The  test  series  were  directed  towards  providing  information  in 
several  areas.  These  were,,  for  the  axisymmetric  turbulent  flow 
case  (1)  does  the  polymer  affect  the  point  of  separation  thereby 
affecting  the  form  drag,  (2)  qualitatively,  is  there  any  change  in 
the  boundary  layer  characterization  and  lastly,  (3)  what  is  the 
wall  concentration  necessary  to  achieve  a maximum  skin  friction 
reduction  and  what  is  the  diffusion  process  for  an  ejector  ejecting 
into  the  laminar  nose  region.  A discussion  of  the  effects  of  poly- 
mers on  the  separation  region  and  the  effects  of  polymer  on  the 
boundary  layer  characteristics  will  be  performed  simultaneously 
since  they  relate  to  the  flow  visualization  portion  of  the  study. 

B . Boundary  Layer  Characterization 
and  Separation  Tests 

Both  the  simple  dye  ejecting  body  and  the  polymer  ejecting 
body  were  used  in  this  series  of  tests.  The  mechanism  of  main- 
taining a high  wall  concentration  was  different  in  these  tests 
The  dye  ejecting  body  was  operated  in  the  drop  tank  with  a uniform 
polymer  concentration  throughout  the  tank.  The  polymer  ejecting 
body  depended  on  a diffusion  process  for  wall  concentration  control. 
Direct  comparison  can  be  made  which  clearly  displays  the  effect 
of  the  polymers  on  the  character  of  the  boundary  layer. 
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Brady  (1973)  summarized  the  work  of  Kline  et  al  (1963)  quite 
nicely.  He  states: 

"Kline  et  al  found  the  'laminar  sublayer'  to  be 
made  up  of  a regular  structure  of  low  and  high 
velocity  longitudinal  streaks  which  meander 
transverse  and  normal  to  the  wall.  A dimension- 
less streak  spacing  \+  = \v*  ° 76.5  was  found 
for  zero  pressure  gradient.  These  streaks 
either  break  up  or  randomly  burst  from  the  sub- 
layer into  the  fully  turbulent  region. 

'Hie  fully  turbulent  region  is  one  of  intense 
mixing  and  high  dissipation  of  energy.  Pro- 
truding from  this  region  are  intermittent 
large  eddies  — visualized  as  peninsulas  of 
turbulence . 

The  peninsulas  of  turbulence  which  extend 
into  the  third  boundary  layer  zone  — the 
outer  turbulent  region  — gives  it  a charac- 
teristic not  unlike  the  wake  behind  a cylin- 
der. It  has  rather  large  lumps  of  turbulence 
at  intermittent  spaclngs. 

It  is  clear  that  the  driving  force  for  the 
entire  turbulent  boundary  layer  is  the  gene- 
ration of  velocity  streaks  in  the  sublayer 
and  their  subsequent  bursting  outwards. 

Kline  presents  the  results  of  many  other 
investigators,  as  well  as  his  own,  in 
support  of  this  hypothesis." 

Figure  44  is  a pictorial  sketch  of  the  boundary  layer  flow  taken 
from  Kline  (1963)  displaying  the  process  quite  nicely.  The 
subsequent  test  results  bear  this  process  out  quite  well  and 
display  the  effects  of  polymers  on  the  process. 

The  first  series  of  tests  were  made  with  the  dye  ejecting 
bodies.  In  this  series,  one  or  more  of  the  four  dye  ejecting 
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Figure  44.  Pictorial  sketch  of  boundary  layer  flow  (from  Kline,  (19631  ) 


holes  was  plugged  to  attempt  to  minimize  tank  contamination.  The 
tests  were  of  the  launcher  or  free  drop  type  depending  on  the  tail 
conf  lgur.it  ion. 

Figure  45  displays  the  boundary  I aver  character  for  the  (5° 
tail  configuration  operating  fn  water  at  a Reynold's  number  of 
4.94  x 10*\  The  highly  turbulent  character  with  the  ejected 
eddies  is  clearly  evident.  Measurements  of  the  boundary  layer 
taken  from  this  picture  and  others  not  displayed,  at  a mean 
position  of  the  ejected  eddies,  is  shown  in  Figure  46.  The  analyti- 
cal program  previously  discussed  was  exercised  to  predict  the 
boundary  layer  thickness.  This  is  shown  in  Figure  46  also.  Good 
agreement  occurs,  within  10%  for  this  case,  if  one  considers  that 
the  measured  position  was  the  mean  of  the  eddies  and  not  the  tips 
as  theorized  by  Kline  (1963).  Figure  47  presents  a normalized 
version  of  this  same  information.  both  the  calculated  and  measured 
values  were  normalized  by  their  respective  values  at  x/L  = .5. 

As  is  shown  in  the  figure,  an  excellent  prediction  of  the  boundary 
layer  shape  as  a function  of  distance  down  the  body  results  for 
this  body  by  the  analytical  routine. 

A second  test  was  performed  with  the  12°  tail  configuration 
at  a length  Reynold's  number  of  3 x 10^  and  a polymer  concentration 
of  2.5  WPPM.  The  forebody  portion,  shown  in  Figure  48  displays 
the  same  character  as  the  no-polymer  case  seen  previously.  The 
ejected  eddy  height  seems  somewhat  larger  In  this  case  probably 
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Figure  45. 


Photograph  of  dve  ejection  into 
6°  tail  - no  polymer  1h>  = 4.94 


the  bound. ire 
x I0h 


layer  - 


NORMALIZED 

BOUNDARY  LAYER  THICKNESS 


Figure  46.  Boundary  layer  thickness  vs  normalized  axial  position  - 
6°  tail  configuration  water  case 


0.2  0.4  0.6  0.8  1.0  1.2 

NORMALIZED  AXIAL  POSITION  X/L 

Figure  47.  Normalized  boundary  layer  thickness  vs  normalized  axial 
position  - 6°  tail  configuration  water  case 
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due  to  the  suppression  of  the  smaller  scale  turbulence  as  displayed 
by  White,  A.  Figure  49  displays  the  measured  and  calculated  boundary 
layer  thickness  again  using  the  mean  eddy  height  for  measurement 
purposes.  The  thickness  measured  in  this  test  is  believed  to  be 
in  error  (the  body  rotated)  since  the  first  visible  position  of 
the  dye  is  too  far  for  the  nose.  The  shape  of  the  profile  should 
be  correct,  however.  The  normalized  values  are  displayed  in  Figure 
50  again  showing  reasonable  shape  agreement. 

The  effect  of  polymer  suppression  of  the  small  scale  turbu- 
lence was  much  more  apparent  in  the  tests  with  the  polymer  ejecting 
body.  The  polymer  solutions  used  in  these  tests  were  all  made  up 
within  a 24-hour  period.  Table  7 lists  the  pertinent  data  for 
this  test  series. 

Figures  51,  52  and  53  display  photographs  of  the  polymer  eject- 
ing body  ejecting  visible  dye  and  water. 

Table  7 

Polymer  Ejecting  Body  - Boundary  Layer  Visualiza- 
tion and  Measurement  Tests  - Test  Parameters 

POLYMER  CONCENTRATION 
(WPPM) 

Water 

50 

500 


VELOCITY 

LENGTH 

(FT/SEC) 

REYNOLDS  NUMBER 

26.9 

5.33  x 106 

28.33 

4'.  67  x 106 

29.1 

5.31  x 106 

29.2 

4.9  x 106 

1000 


NORMALIZE 

BOUNDARY  LAYER  THICKNESS 


Figure  49.  Boundary  layer  thickness  vs  normalized  axial  position  - 
12°  tail,  2.5  WPPM  polymer 


Figure  50.  Normalized  boundary  layer  thickness  vs  normalized  axial  position 
12°  tail  configuration,  2.5  WPPM 
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These  views  clearly  show  the  coarse  and  f ine  scale  turbulent  type 
flow.  Measurements  of  the  boundary  layer  are  plotted  in  Figure  60. 
Also  to  be  noted  in  Figure  51  is  the  uniform  ejection  process  which 
apparently  is  not  seriously  effecting  the  flow.  From  these  views, 
the  flow  is  apparently  fully  turbulent  near  the  nose  of  the  body. 

The  tests  were  repeated,  this  time  ejecting  a 50  WPPM  polymer 
solution.  Figures  54  and  55  display  the  results.  The  fine  scale 
eddy  structure  is  apparently  missing  leaving  only  the  coarse  struc- 
ture. More  interestingly,  the  boundary  layer  thickness  has  decreased 
a disproportionate  amount  in  the  midbody  region  as  compared  to  the 
tail  suggesting,  as  it  should  be,  higher  effective  wall  concentrations 
forward  on  the  body  which  result  in  thinner  boundary  layers  in  this 
region  and  a rapid  growth  in  the  tail  region.  Since  the  boundary 
layer  is  beginning  to  get  rather  thin  in  the  forebody  region,  the 
measurement  method  used  should  be  described.  Slides  have  been 
made  of  the  figures  presented  in  this  study,  these  were  projected 
on  a screen  such  that  the  body  diameter  (3  inches)  projected  to 
about  2 feet  in  diameter.  The  resultant  boundary  layer  thicknesses 
near  the  nose  were  a projected  .2  inch.  An  overlay  slide  with  one 
inch  axial  marks  was  used  as  an  overlay  (with  2 projectors  to 
insure  proper  scale).  This  was  required  since  the  body  does  not 
reach  3 inch  diameter  for  some  distance  form  the  nose. 

The  last  two  groups  of  tests  in  this  series  were  with  500 
WPPM  and  1000  WPPM.  Figures  56  and  57  show  the  forebody  view 
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Figure  54.  Photograph  of  polymer  ejecting  body  - midbody  view 
50  WPPM  polymer  ejection  Re  = 4.67  x 10 
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Figure  55.  Photograph  of  polymer  ejecting  body  t<ul  view 
50  WPPM  polymer  ejection  Re  - 4.65  x I.U 
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and  tail  view  for  the  500  WPPM  case  and  Figures  58  and  59  show 
similar  views  for  the  1000  WPPM  case.  Apparently,  a limit  of 
drag  reduction  is  being  reached,  presuming  this  may  be  inferred 
from  the  fact  that  there  is  little  difference  in  the  boundary  layer 
reduction  between  these  cat.es  as  seen  in  Figure  60.  Moreover, 
an  additional  phenomenon  is  occurring.  It  is  first  noticeable  in 
Figure  57  and  readily  apparent  at  the  1000  WPPM  level  of  Figures 
58  and  59.  A significant  number  of  streaks  are  evident  in  the 
flow.  The  approximate  streak  spacing  is  about  .015  inches.  The 
number  of  ejector  holes  is  approximately  700.  For  the  body  diameter 
of  3 inches,  the  number  of  streaks  would  be  about  620.  Additionally, 
applying  the  dimensionless  streak  spacing  found  by  Kline  et  al 
(1963)  in  the  laminar  sublayer,  for  a v*  of  approximately  .6 
ft/sec  and  viscosity  of  about  that  of  water,  results  in  a "laminar 
sublayer"  streak  spacing  of  the  order  of  .015  inches.  Apparently, 
at  these  concentration  levels,  the  turbulence  intensity  is  so  damped 
that  the  streaks  become  well  ordered  and  have  difficulty  in  bursting 
out  into  the  other  layers  of  the  boundary  layer.  It  is  not  until 
the  tail  region  that  a very  filamentary  bursting  is  occuring 
suggestive  of  the  long  stringiness  associated  with  elastic  polymers. 
Even  the  angle  of  the  streaks,  essentially  in  the  axial  direction, 
are  suggestive  of  very  low  energy  levels  bursting  the  streaks  out- 
ward as  compared  to  the  water  case.  The  implication  of  these 
photographs  is  that  turbulent  diffusion  has  ceased  and  a molecular 
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Photograph  of  polymer  ejecting  body 
1000  WPPM  polymer  ejection  Re  = 4. 


Figure  58 


Figure  59 


Photograph  of  polymer  ejecting  oody  - tail  view 
1000  WPPM  polymer  eieetion  Rc  = 4.9  x 10^ 


diffusion  only  in  the  sublayer  is  occurring.  This  would  suggest 
extremely  long  mixing  lengths.  Tests  by  Nadolink  (1968)  displayed 
the  same  characteristic  but  for  polymer  ocean  type  tests.  Figure  61 
displays  a normalized  plot  of  the  boundary  layer  growth  for  these 
tests.  The  polymer  cases  clearly  show  a delayed  growth  trend 
indicative  of  suppressed  diffusion. 

The  last  series  of  visual  observations  centered  around  the 
question  of  whether  application  of  polymers  effected  the  separa- 
tion point  at  the  tail.  Several  tests  were  run  with  the.  dye 
ejecting  body  with  the  hemispherical  tail  in  fresh  water  and  in  a 
20  WPPM  polymer  ocean.  The  length  Reynolds  number  for  both 
cases  was  1.877  x 10^. 

Figure  62,  the  test  with  plain  water,  displays  a thick 
boundary  layer  near  the  tail  and  a separation  angle  measured  from 
the  vertical  of  approximately  62°  to  65°.  Figure  63  displays 
the  results  of  a similar  test  with  20  WPPM  polymer  concentration  in 
the  tank.  The  boundary  layer  is  characteristically  thin  but  the 
separation  point  resulting  in  an  angular  measurement  of  80.6°  for 

either  polymer  case  or  water  case.  On  a nondimensional  basis,  the 
computer  program  predicted  separation  at  x/L  = .925  while  the 
data  indicate  x/L  = .948  for  separation.  Nonetheless,  the 
polymer,  for  the  case  tested,  did  not  apparently  change  the  sepa- 
ration point.  Similar  tests  with  the  12°  tail  configuration 
also  showed  no  effect  on  the  separation  point.  Since  sphere 
tests  conducted  by  Lang  and  Patrick  (1967)  indicated  that  sphere 
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NORMALIZED 

BOUNDARY  LAYER  THICKNESS 


Figure  GO.  Boundary  layer  thickness  vs  normalized  axial  length  - 
polymer  ejecting  body 


Figure  61.  Normalized  boundary  layer  thickness  vs  normalized  axial  length 
polymer  ejecting  body 
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1’botograpl)  of  dye  ejecting  body  - spherical  tail 
tail  view  - water  case  Re  = l-877  x 11)1 
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drag  was  reduced  by  a movement  of  the  separation  point  from  the 
laminar  separation  point  to  a further  back  turbulent  separation 
point.  These  results  are  not  surprising  for  the  case  of  a turbu- 
lent separation  point  initially.,  These  results  are  consistent 
with  those  of  Stone  and  Elliott  (1970).  Unfortunately,  this 
small  change  in  predicted  separation  point  places  the  separation 
at  the  extreme  of  the  negative  pressure  coefficient,  Figure  37, 
resulting  in  a very  high  calculated  form  drag. 


C.  Drag  Reduction  Tests 

The  next  series  of  tests  were  directed  towards  determining 
drag  reduction  efficiency  for  various  polymer  concentrations  and 
tail  configurations  and  to  test  the  analytical  model  for  the  6° 
tail  case.  The  form  drag  calculation  for  this  case  are  simpler  due 
to  the  geometry  and  pressure  coefficient  profile.  The  test  matrix 
run  was  as  outlined  in  Table  5.  A total  drag  coefficient  was 
calculated  for  each  series  of  tests  using  the  relation 


CD  = ~ w - B ~ U14) 

1/2  p Af  U02 

where  the  frontal  area,  Af  = 0.0491  ft2.  From  the  drag  coefficient, 
Cq,  the  percentage  total  drag  reduction  was  then  calculated  for 
each  polymer  concentration  tested.  The  percentage  total  drag 
reduction  is  given  by 


% TDR  = 100  (1  - 


Cdpoly-, 

C% 


(115) 
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where 


Cn  = Drag  coefficient  for  polymer  tests 
uPOLY 

Cj^  = Drag  coefficient  for  water  tests. 

A calculation  of  percent  skin  friction  reduction  was  also  made  for 
each  test.  The  relation  for  percent  skin  friction  reduction  are. 

% SFR  = 100  (1  - CfP0LY)  (116) 

% 

where . 

Df 

cf  = CD  ^ 

and  Op  = measured  drag  coefficient 

Dp  = total  shear  drag  analytically  calculated 
Dt  = W-B  = body  weight  m water. 

It  is  noted  here  that  Dt,  the  actual  drag,  is  used  in  the  / SFR 
calculation  and  not  the  total  calculated  drag.  This  is  done 
since  the  skin  friction  as  calculated  by  this  method  has  been 
shown  by  White  (1972)  to  be  accurate.  The  addition  of  the  term 
for  calculating  form  drag  is  only  considered  to  be  an  indicator 
and  may  be  considerably  in  error  due  to  the  pressure  coefficient 
inaccuracies.  The  6°  tail  tests  were  the  only  tests  where  the  launcher 
could  be  used.  These  tests,  therefore,  were  the  only  ones  in  this 
series  where  terminal  velocities  were  achieved. 

Extreme  Instabilities  were  noted  when  launch  tests  were  attempted 
with  the  12°  tail  and  spherical  tail  models.  This  is  attributed  to 
the  launcher  receptacle  being  designed  to  accept  the  6°  tail.  As 
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a result,  the  12°  Lail  and  spherical  tail  were  not  properly 
guided  during  the  launch  phase.  Even  with  the  weighted  g tide 
wire,  photographs  of  the  test  models  travelling  at  angles  of  near 
10°  relative  to  the  wire  were  obtained.  Drag  data  obtained  in 
this  manner  were  not  of  any  use.  Additionally,  velocity  measure- 
ments were  never  achieved  when  this  occurred  since  the  guide  wire 
would  generally  trip  the  lasers  during  the  launch  transient. 

The  free  drop  12°  and  spherical  tail  tests  were  projected  to 
terminal  velocity  using  a procedure  outlined  in  Nadolink  (1968). 

The  procedure  solves  the  equation  of  motion  for  an  accelerating 
body  having  as  knowns  the  weight  in  air,  the  weight  in  water,  the 
velocity  at  a particular  distance  and  an  entrained  mass  coeffi- 
cient emperically  determined  and  plotted  in  the  report  as  a function 
of  L/D.  The  procedure,  because  of  the  form  of  the  final  relation 
is  iteratively  solved  for  terminal  drag  coefficient,  which  will 
predict  the  velocity  at  the  station  which  is  being  measured.  From 
this,  equation  (114)  is  used  to  calculate  the  terminal  velocity. 

Tables  8,  9 and  10  present  the  data  obtained  in  this  series  of 
tests  as  well  as  the  percent  total  drag  reduction  and  skin  fric- 
tion reduction  calculated.  Other  pertinent  information  noted  is 
the  total  drag  calculated  from  the  analytical  routine  for  the  iden- 
tical test  conditions.  As  may  be  noted  from  Table  8,  the  total 
drag  predicted  by  the  analytical  ri  utine  is  within  10%  for  the 
6°  tail  case.  Additionally,  the  total  skin  friction  determined 
from  equation  (89)  results  in  a friction  drag  coefficient  based 
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W-B  =3.47  LBS.  POLYMER  TESTS  (WEIGHT  IN  WATER) 


TABLE  9 

DATA  SUMMARY  12°  Tail  Tests 
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W-B  = 7.109  LBS.  (WEIGHT  IN  WATER) 


TABLE  10 
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on  surface  area  of  Cf  = .00311  which  compares  very  favorably  with 
that  determined  from  the  accepted  formulas  for  flat  plates  of  .00336 
at  the  same  length  Reynolds  number.  The  calculated  skin  friction 
coefficient  is  plotted  in  Figure  64  taken  from  Granville  (1971). 

For  the  tS/r0  being  worked  with  in  this  study  and  for  a mild  pressure 
gradient,  a comparison  is  considered  reasonable. 

Figure  65  presents  a plot  of  the  percent  total  drag  reduction 
for  the  various  tail  configurations.  As  noted  in  the  figure,  peak 
drag  reductions  are  occurring  for  the  6°  and  12°  tail  case  at 
polymer  concentrations  of  between  15  and  20  WPFM.  The  effect  of 
increased  form  drag  with  the  12°  and  spherical  tail  models  results 
in  a decreased  total  drag  reduction  since  the  polymer  only  reduces 
the  skin  friction.  Figure  66  is  a plot  of  percent  skin  friction 
reduction.  This  plot  clearly  displays  the  achievement  of  equal 
skin  friction  reduction  on  all  the  bodies  tested  as  would  be 
expected.  Calculation  of  the  skin  friction  coefficient  based  on 
surface  area  for  the  20  WPPM  case  with  the  6°  tail  results  in  a Cf 
of  .0011  at  a length  Reynolds  number  of  4. 85  x 10^.  This  point  is 
plotted  on  Figure  64  and  is  in  good  agreement  with  the  maximum 
drag  reduction  value  given  by  Granville.  Plotted  on  Figure  64 
also  are  the  skin  friction  coefficients  calculated  for  the  12° 
and  spherical  tail.  These  also  are  in  reasonable  agreement  with 
the  maximum  drag  reduction  line  as  proposed  by  Granville. 
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SPHERICAL  TAIL 


D . Polymer  Concentr ation  Profile  Tests 


The  polymer  ejecting  body  was  used  throughout  this  entire 
series  of  tests.  The  primary  emphasis  in  these  tests  was  to 
obtain  polymer  concentration  measurements  at  the  wall  and  at  two 
heights  in  the  boundary  layer.  These  heights  were  selected  to  be 
.025  inches  from  the  body  to  the  centerline  of  the  probes  and  .055 
inches  from  the  body  to  the  probe  centerline  of  the  probes  for 
reasons  previously  discussed.  The  probes  are  as  described  pre- 
viously and  ohown  in  Figure  20.  Only  one  stroke  rate  was  used 
in  these  tests  although  other  stroke  rates  are  possible  with  the 
model  through  an  orifice  change.  The  stroke  rate  selected  resulted 
in  a volume  rate  of  flow  of  20.6  in^/sec.  Polymer  concentrations 
of  50  and  500  WPPM  of  Polyox-WSR-30  were  ejected  and  samples  taken 
at  the  wall  and  the  two  probe  heights.  Testing  with  5 WPPM,  10 
WPPM,  20  WPPM  and  1000  WPPM  was  also  performed  but  only  with  the 
flush  probes  installed.  Tests  were  also  performed  with  fresh 
water  dyed  with  the  fluorescent  dye  to  determine  what  differences 
in  diffusion  exist  between  the  polymer  contaminated  solution  being 
ejected  and  the  water.  All  tests  were  run  using  the  launcher  and 
no  attempt  was  made  to  attain  terminal  velocity.  The  tests  were 
all  made  at  a high  enough  initial  velocity  that  near  terminal 
velocities  existed  and  transient  characteristics  were  deemed  to 
be  unimportant.  The  tests  could  not  be  used,  however,  to  determine 
drag  reduction  or  drag  coefficients.  The  tests  went  very  smoothly 
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with  only  one  or  two  failures  in  the  entire  test  series  when  a 


launch  was  made.  Velocity  variations  for  similar  tests  were 
within  one  percent.  Unfortunately,  concentration  measurements 
were  more  diverse  with  10  to  15%  variations  not  uncommon  although 
the  bulk  of  the  data  fell  in  a narrowband.  The  test  operation  was 
very  procedure  oriented  with  deviations  generally  resulting  in  an 
aborted  attempt  at  a launch.  The  long  hours  spent  in  calibrating 
the  model  was  key  to  the  smooth  test  series  since  carefully  followed 
procedures  resulted.  The  test  matrix  run  is  that  shown  in  Table  6. 

In  general,  each  test  series,  a series  being  either  no  polymer  or  a 
specific  polymer  concentration,  began  with  the  flush  probes  installed. 
The  .025  inch  probes  were  next  installed  and  finally  the  .055  inch 
probes.  No  intermixing  of  probe  heights  occurred,  even  on  a per 
station  basis.  The  test  model  functioned  perfectly,  as  designed, 
with  over  several  hundred  firings  being  accounted  for  considering 
stroke  rate  calibration  tests,  polymer  concentration  tests,  photo- 
graphs firings  and  actual  data  gathering  firings.  The  only  malfunc- 
tion occurring  was  in  the  firing  rod.  One  of  the  functions  of  this 
rod  was  to  follow  the  polymer  ejection  piston  during  the  ejection 
process.  A photograph  of  the  tail  region  of  the  model  could  then 
be  analyzed  to  reveal  the  position  of  the  firing  rod  still  external 
to  the  body.  This  information  combined  with  body  velocity  and 
distance  to  the  camera  station  would  allow  a continual  check  on 
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t lie  stroking  rate  to  be  made.  Unfortunately,  the  launch  process 
proceeded  to  distort  the  rod  slightly  on  each  launch  resulting 
in  it  becoming  stuck  and  not  properly  following  the  ejection  piston. 

By  inserting  a spacer  in  the  ejector  piston  volume  to  limit 
its  stroke,  an  independent  check  on  stroke  rate  was  made.  With 
the  appropriate  spacer  size  for  the  body  velocity,  camera  station 
and  stroke  rate  installed,  a test  was  made  with  visible  dye  and  a 
photograph  taken.  The  photograph  displayed  several  inches  of  the 
forebody  free  from  dye  while  the  remainder  of  the  body  uniformly 
covered  w'th  dye.  This  verified  that  the  proper  ejection  rate  was 
occurring  under  test  conditions. 

The  effects  of  the  probes  on  the  boundary  layer  flow  was  also 
of  concern.  The  rotation  of  subsequent  sampling  stations  on  the 
body  was  a result  of  this  concern.  Figure  67  displays  a photograph 
of  a .025  inch  probe  within  a dyed  boundary  layer  for  a water  test 
case.  No  significant  distortion  of  the  boundary  layer  is  apparent. 
The  velocity  in  this  case  was  20.5  ft/sec. 

The  results  of  the  test  series  are  tabulated  in  Tables  11,  12, 
13,  14,  15,  16  and  17.  Although  polymer  concentrations  were  not 
directly  measured,  the  assumption  that  the  Uranine-B  in  solution 
reacts  and  diffuses  in  the  same  manner,  as  previously  discussed, 
allows  direct  translation  to  polymer  concentration  to  be  made. 

Only  for  the  water  case  is  the  tracer  concentration  discussed. 
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Figure  67.  I’liot ograph  of  i neii  probe  within 
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Water  E j ec  t Ion  Test  Results 

Figures  68,  69  and  70  display  the  data  obtained  in  test  series  1, 
2 and  3 of  Table  10.  Referring  first  to  Figure  10,  it  can  be 
noticed  that  for  this  water  diffusion  case,  the  wall  dye  concentra- 
tion has  decreased  to  .1  of  the  ejected  concentration  in  a distance 
of  slightly  less  than  one  body  diameter.  The  diffusion  over  the  next 
nine  inches  of  body  length  near  halves  again  the  concentration  level. 
The  diffusion  process  is  quite  rapid,  characteristic  of  turbulent 
mass  transport.  The  predicted  curve  on  this  figure  applies  the  two 
zone  diffusion  boundary  layer  growth  model  combined  with  the  axisym- 
metric  boundary  layer  model  as  discussed  earlier.  The  concentration 
similarity  profiles  for  water  were  applied  for  the  two-zone  process 
quite  successfully. 

Figure  69  displays  the  data  from  the  series  2 experiment  with 
the  .025  inch  probes  sampling  the  boundary  layer.  Again,  comparison 
with  the  model  for  the  water  case  seems  quite  good.  Although  not 
apparent  in  the  data,  the  computer  solutions  shows  the  concentration 
growth  prior  to  the  sampling  station  and  then  a normal  decay  as  the 
dye  diffuses  outward  in  the  boundary  layer. 

Figure  70  displaying  series  3 tests  with  water  ejection  and  the 
.055  inch  probes,  the  rise  in  the  concentration  levels  between 
the  first  and  second  probe  stations  is  evident  in  the  data  and 
predicted  by  the  theory. 
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BOUNDARY  LAYER  TO  INJECT 
CONCENTRATION  RATIO  C/C, 


Figure  70.  Boundary  layer  to  injected  concentration  ratio  vs  normalized 
axial  length  Test  Series  .‘1 
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Figure  71,  a composite  plot  of  the  calculated  flush,  . 025-inch 
and  .055-inch  probes  displays  an  interesting  result.  First,  the 
increase  in  concentration  initially  for  the  probes  away  from  the 
wall  is  evident.  The  decay  in  concentration  as  the  fluorescent 
dye  diffuses  outward  is  as  expected  In  the  approximately  constant 
diameter  section.  At  x/L  of  approximately  .6,  however,  the  concen- 
tration ratio  begins  to  increase.  The  boundary  layer  in  this  region 
is  growing  very  rapidly  as  both  the  photograph  for  this  case,  Figure 
53  and  the  calculated  and  measured  boundary  layers  indicate,  Figure 
60.  The  cause  of  this  rise  is  the  decreasing  diameter  resulting  in 
a smaller  boundary  layer  flux  area  within  which  the  tracer  resides. 
This  is  a rather  interesting  result  since  in  external  flows  one  would 
expect  a uniform  decrease  in  concentration  with  length  due  to  boundary 
layer  growth,  Figure  72  displays  the  growth  of  the  diffusion  boundary 
layer,  defined  as  height  at  which  the  concentration  is  .5  of  the  wall 
concentration,  within  the  normal  boundary  layer.  The  growth  of 
the  diffusion  boundary  layer  is  reduced  as  compared  to  the  normal 
boundary  layer  in  the  tail  region,  x/L  > .6. 

Figure  73  displays  a plot  of  the  concentration  profile  for  the 
first  two  probe  stations  normalized  to  the  wall  concentration.  The 
remaining  two  stations  were  not  plotted  since  only  the  .055  probe 
heighr  gave  a value  of  C/Cy  different  from  one.  A single  point 
curve  would  have  resulted.  Table  18  provides  a tabulation  of  the 
experimental  data  used  in  generating  this  figure.  As  noted  in  the 
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Figure  72.  Calculated  diffusion  boundary  layer  and  boundary  layer  growth  vs  normalized  axial  length 


figure,  the  diffusion  process  Is  quite  rapid  with  wall  concentrat  Ion 
being  approached  quite  far  out  Into  the  boundary  layer  a short 
distance  back  on  the  body.  By  defining  the  diffusion  boundary  layer 
thickness  as  that  value  of  y at  which  C/Cy  equals  .5,  a normalization 
height  is  determined  for  plotting  this  data  on  a similarity  concen- 
tration profile  plot.  Table  19  presents  these  data,  along  with  data 
from  subsequent  tests  with  polymers.  A comparison  of  these  data 
(water  case)  with  that  from  other  experimentators  is  done  on  the 
similarity  profile  plot  of  Figure  74.  The  dotted  line  on  this  figure 
taken  from  Hsu  (1971)  has  been  generated  from  data  by  Porch  and  Cermak 
(1964)  and  verified  by  Wetzel  and  Ripkin  (1970).  The  equation  for 
this  line  is  as  shown  on  the  figure.  The  limited  data  obtained  in 
this  experiment  verify,  for  the  axisymmetrir  case,  the  applica- 
bility of  this  similarity  profile  for  the  downstream  region.  The 
water  ejection  line,  although  falling  slightly  higher  resulting  in 
a higher  value  of  the  exponent  in  the  C/(\,  relation,  is  in  good 
agreement  with  the  previous  investigation.  The  applicability  of  the 
model  in  predicting  the  wall  concentration  and  concentration  profiles 
in  the  boundary  layer  seems  apparent. 

Polymer  Ejection  Test  Results 

Figures  75,  76  and  77  dinpiay  the  data  obtained  in  test  series 
4,  5 and  6 of  Table  12.  In  this  series,  50  WPPM  polymer  concentra- 
tion was  ejected.  The  results  were  not  as  expected.  At  all  probe 
stations  the  polymer  concentration  remained  high  indicating  no  diffu- 
sion occurring.  As  may  be  seen  in  Figure  75,  the  combined  model 
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Figure  74.  Concentration  profiles  determined  for  water  and  polymer  ejection 
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Figure  75.  Boundary  layer  to  injected  concentration  ratio  vs  normalized 
axial  length  Test  Series  4 
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Figure  7G.  Boundary  layer  to  injected  concentration  ratio  vs  normalized 
axial  length  Test  Series  5 
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Figure  77.  Boundary  layer  to  injected  concentration  ratio  vs  normalized 
axial  length  Test  Series  G 
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predicts  a short  initial  zone.  Figures  76  and  77  show  the  concentra- 
tion measured  with  the  .025-inch  probes  and  the  .055-inch  probes. 

The  experimental  data  display  a slow  build  up  of  the  concentration. 
The  .055-inch  probes  which  should  be  well  outside  the  sublayer, 
display  much  scatter  expecially  at  the  first  station.  Apparently 
there  is  diffusion  from  the  surface  of  the  viscous  layer  otherwise 
the  measured  value  would  be  much  lower.  The  calculations  do  not 
account  for  diffusion  from  the  viscous  layer.  There  may  still  be 
some  bursting  occurring  at  this  concentration  level  as  some  of  the 
polymer  in  the  ejection  process  is  being  hurled  to  regions  further 
from  the  wall.  The  calculated  values  displayed  in  this  figure  show 
the  characteristic  rise  and  fall  off  of  the  rapid  diffusion  process 
as  is  predicted  by  the  intermediate  and  final  zone  models.  Figure 
77  shows  a sudden  rise  in  the  data  at  the  last  measurement  station. 
Quite  possibly,  the  end  of  the  initial  zone  is  being  reached  and 
the  more  rapid  diffusion  is  beginning  to  occur. 

Table  18  summarizes  the  concentration  profile  data  taken  for 
the  50  WPPM  case  from  which  Figure  78  is  plotted.  This  figure 
displays  the  totally  changed  character  of  the  concentration  profile 
from  the  water  case.  Figure  73.  Comparing,  <$<j,  the  value  of  y at 
C/Cy  equal  to  .5,  for  the  water  case  6^  is  .044  inches  at  the  first 
measurement  station,  while  for  the  50  WPPM  case  it  is  .0055  inches, 
(nearly  a factor  of  10  difference)  . Table  19  presents  the  data  taken 
from  Figure  78  from  which  a redefinition  of  the  distance  from  the 
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Figure  78.  Distance  from  wall  vs  normalized  concentration  ratio 
50  WPPM  case 
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wall,  y,  to  a nnnd Imena i orm'l  dlntanco,  In  mailt-.  Tlila  mm- 

dlmensional  height  la  plotted  In  Figure  lb  and  a curve  fit  made. 

The  result,  is  a similarity  profile  with  an  exponent  of  .75  for  all 
measurement  stations.  This  exponent  implies  significantly  reduced 
diffusion  but  not  to  the  level  which  the  data  display.  The  dis- 
crepancy lies  in  the  fitting  of  the  curve  in  Figure  78.  Three 
points  are  available.  These  are  at  y = .025  inches,  y = .055 
inches  and  y = 0 inches.  Clearly,  the  curve  could  be  drawn  a number 
of  ways  which  would  significantly  reduce  the  exponent  value  by 
reducing  5^.  The  information  demonstrating  low  diffusion  rate  is 
contained  in  the  wall  concentration  measurement,  however,  so  no 
attempt  was  made  to  match  the  curves  to  the  data  through  the  com- 
puter routine.  In  fact,  the  variable  K3  model  which  resulted  in  a 
steady  wall  concentration  profile  calculated  an  exponent  value  of 
about  .15.  Let  it  suffice  that  Figure  78  displays  markedly  reduced 
diffusion  to  that  of  water  and  recommend  for  future  work  that  a probe 
be  tested  at  .01-inch  height.  Another  interesting  facet  of  these 
tests  were  that  the  exponent  did  not  change  with  measurement  station. 
This  implies  that  the  diffusion  rate  does  not  change,  in  these 
tests,  with  downstream  distance.  Again  this  may  be  a function  of  the 
curve  fit  only  since,  unless  only  molecular  diffusion  is  occurring, 
it  would  be  expected  that  the  diffusion  rate  would  increase. 

Figures  79,  80  and  81  present  similar  test  results  for  the  test 
series  7,  8 and  9 of  Table  13.  As  can  be  seen  in  Figure  79,  for  the 
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FLUSH  PROBES 
500  WPPM  POLY  MER 
TEST  SERIES  7 
TABLE  12 
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□ TEST  27 
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Figure  79.  Boundary  layer  to  injeeted  concentration  ratio  vs  normalized 
axial  length  Test  Series  7 
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Figure  80.  Boundary  layer  to  injected  concentration  ratio  vs  normalized 
axial  length  Test  Series  8 


500  WPPM  case,  both  the  calculated  and  experimental  data  clearly 
indicate  little  diffusion  is  occurring.  The  .025-inch  probes 
display  the  same  gentle  rise  in  concentration,  Figure  80,  as 
compared  to  the  decay  in  concentration  noted  for  the  water  case, 

Figure  69,  resulting  from  the  higher  diffusion  in  that  case.  The 
.055-inch  probes  again  display  a slow  rise  with  low  levels  indicative 
of  low  diffusion.  It  is  interesting  to  note  the  scatter  at  the  first 
station  implying  some  anomaly  in  the  ejection  process.  No  calculations 
are  evident  for  the  concentration  ratio  in  these  probe  figures  since 
the  molecular  diffusion  portion  of  the  model  does  not  incorporate 
a similarity  concentration  profile.  The  variable  K3  model  would 
account  for  this.  Before  turning  to  a concentration  profile  dis- 
cussion, a test,  test  series  10,  Table  14,  was  performed  with  1000 
WPPM  ejected.  Figure  82  displays  the  results.  Again  a nearly 
constant  concentration  with  body  distance  results  Indicating  little 
or  no  diffusion.  This  is  not  surprising  considering  the  photographs 
taken  showing  the  boundary  layer  structure,  Figures  58  and  59. 

The  data  for  the  500  WPPM  series  have  also  been  normalized  for 
plotting  as  concentration  ratio  vs  height,  Table  18.  Figure  83 
displays  this  plot.  Again  an  extremely  low  diffusion  rate  is 
evident  compared  to  the  water  case,  Figure  73.  The  500  WPPM  case  is 
much  steeper  in  profile  than  that  of  the  50  WPPM  case,  Figure  78, 
implying  an  even  greater  decrease  in  diffusion  rate.  Interestingly, 
wher.  the  height,  y,  is  nondimensionalized  with  6<j,  Table  19,  the 
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Figure  81.  Boundary  layer  to  injected  concentration  ratio  vs  normalized 
axial  length  Test  Series  9 


Q 

LU 

h-  o' 
O 'x 
UJ  o 


Z O 


25 

V « 
CC  I- 

< z 

a Ui 
Z o 
=>  Z 

o o 

CD  O 


— 

- 

FLUSH  PROBES 

1000  WPPM  POLYMER  EJECTION 

TEST  SERIES  10 

- 

TABLE  13 

LEGEND 

O TEST  37 

□ TEST  38 

— 

& TEST  39 

— 

CALCULATED 

/ 

- 

0 / 0 

- 

0 ^ ~~~ ' — 0 

0 r 

O 

A 

- 

MEASURED 

1 1 1 1 

1 

0.1  0.2  0.3  0.4  0.5  0.6 

NORMALIZED  AXIAL  LENGTH  X /L 


Figure  82.  Boundary  layer  to  injected  concentration  ratio  vs  normalized 
axial  length  Test  Series  10 
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Figure  83.  Distance  from  wall  vs  normalized  concentration  ratio 
500  WPPM  case 


2 


the  similarity  profile  is  the  same  as  for  the  50  WPPM  case.  The 
results  are  plotted  in  Figure  74  where  the  polymer  ejection  case 
applies  to  both  the  50  WPPM  and  500  WPPM  tost  series.  Again,  the 
frailties  of  this  plot  are  similar  to  those  of  Figure  78. 

Additional  tests  were  performed  with  reduced  polymer  concentra- 
tions in  an  attempt  to  define  the  position  at  which  the  initial 
diffusion  zone  terminated.  These  test  results  are  displayed  in 
Figure  84,  for  20  WPPM  polymer  solution,  Figure  85,  10  WPPM  polymer 
solution  and  Figure  86,  5 WPPM  polymer  solution  ejected.  For  both 
the  20  WPPM  and  10  WPPM  case,  no  drastic  change  to  a higher  diffu- 
sion rate  is  evident.  The  wall  concentration  to  injected  concentra- 
tion ratio  remaining  about  constant  throughout  the  sampling  region 
of  the  body.  The  5 WPPM  ejected  concentration  did  display  the  result 
being  sought.  As  shown  in  Figure  86,  a definite  reproducible  shift  to 
a lower  concentration  has  occurred  between  the  second  and  third 
sampling  stations.  From  the  flat  plate  data  of  Fruman  and  Tulin 
(1974)  a projected  transition  to  a higher  diffusion  process  would 
have  been  projected  to  occur  at  concentration  below  36  WPPM.  These 
limited  results  indicate,  for  the  axisymmetric  case,  a delay  in 
transition  to  higher  diffusion  occ.rs  at  a significantly  reduced 
concentration  level. 

The  test  series  performed  in  this  study  may  be  compared  directly 
with  those  of  Fruman  and  Tulin  (1974)  if  corrected  for  the  axisymme- 
tric  influence  by  multiplying  the  concentration  ratio  by  D0/Dj_.  The 
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Figure  84.  Boundary  layer  to  injected  concentration  ratio  vs  normalized 
axial  length  Test  Series  11 
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Figure  85.  Boundary  layer  to  injected  concentration  latio  vs  normalized 
axial  length  Test  Series  12 
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Figure  8C>.  Boundary  layer  to  injected  concentration  ratio  vs  normalized 
axial  length  Test  Series  13 


Reynolds  number  ;\mge  for  this  test  series  was  about  5 million 
whereas  the  l'n.avm  and  Tulin  experiments  were  approximately  10 
million.  The  ratio  of  ejection  velocity  to  free  stream  velocity 
(body  velocity  in  this  case)  is  .185  in  this  study  and  ranges  from 
.187  to  .274  in  the  Fruman  and  Tulin  experiments.  Figure  87  displays 
a plot  of  the  concentration  to  injected  concentration  ratio  vs  a 
dimensionless  distance  for  a water  ejected  case  taken  from  Fruman 
and  Tulin  (1974).  The  data  from  this  experiment  compare  well. 
Similar  plots  for  1000  WPPM  Figure  88,  and  500  WPPM,  Figure  89 
show  similar  agreement.  The  data  from  this  experiment  are  believed 
to  be  well  qualified  and  displaying  similar  trends. 

Figure  90,  from  Fruman  and  Tulin  (1974),  displays  the  departure 
from  the  Initial  zone  of  diffusion  as  was  predicted  from  equation 
(96)  from  that  reference.  Data  from  this  experiment  are  plotted 
on  this  figure.  As  is  noted,  all  data  for  concentration  fall  in 
the  initial  zone  region  of  molecular  diffusion  and  have  a corrected 
concentration  ratio  of  about  1.  The  test  series  4 data,  50  WPPM 
case,  falls  to  the  right  of  the  departure  from  the  Initial  zone 
region.  The  reason  is  suspected  to  be  that  equation  (96) , which 
applies  to  slot  ejection,  should  be  modified  to  the  case  of  ejection 
into  the  free  stream  and  becoming  a part  of  the  laminar  boundary 
layer.  The  ejection  velocity  becomes  meaningless  for  this  case 
and  should  be  replaced  with  some  other  velocity  which  may  be 
the  mean  laminar  flow  velocity  or  approximately  .6  of  the  local 
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dimensionless  distance,  x/s,  for  water  injection  (from  Fruman 
and  Tulin  (1974)  ) 


Figure  88.  Ratio  of  concentration  to  injected  concentration,  c/Cj,  vs  dimension- 
less distance,  x/s,  for  1900  ppm  polyox  WSR  301  injection  (from 
Fruman  and  Tulin  (1974)  1 
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free  stream  velocity.  For  this  cast',  equation  (96)  reduces  to 


2.15—  — . (117) 

s Ct 

I 

For  the  50  WPPM  case,  the  reduced  dimensionless  distance  for  the  50 
WPPM  point  would  now  be  placed  at  5.87  rather  than  near  30  where 
it  presently  is.  In  fact,  all  the  points  for  the  present  test  would 
shift  to  the  left  indicating  a large  initial  length  as  the  data 
indicate. 

Figure  91  presents  the  corrected  experimental  data  plotted 
against  the  modified  dimensionless  distance  give  by  equation  (117). 
The  extended  initial  zone  compared  to  the  flat  plate  experiment 
is  apparent  in  this  figure.  The  analytical  model  was  exercised 
for  the  condition  of  the  test  shown  in  Figure  86  with  a transfer 
from  the  initial  to  intermediate  zone  occurring  at  30.  The  dimen- 
sionless distance  given  by  equation  (117)  was  used  to  determine 
the  change  to  the  intermediate  diffusion  zone.  The  calculated 
values  are  displayed  in  Figure  86.  As  may  be  seen,  the  calculated 
intermediate  zone  diffusion  process  is  much  greater  than  the 
experimental  data  indicate  warranting  further  study  of  this  zone. 

The  extended  initial  zone  for  the  ejection  method  applied  is 
significant  in  that  no  special  care  is  needed  In  the  ejection 
process.  Additionally,  this  ejection  process  most  likely  represents 
the  optimal  case  since  higher  ejection  velocities  would  quickly  be 
reduced  to  the  mean  viscous  sublayer  velocity  through  viscous  drag 
and,  for  a slot  ejection  case,  could  result  in  flow  disruption. 
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Figure  91.  Dimensionless  concentration  vs  reduced  dimensionless  distance 
test  results  (from  Fruman  and  Tulin  (1974)  ) 


VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSION 

1.  General 

The  object  of  this  research  was  threefold:  1)  to  develop  methods 

for  wall  concentration  and  drag  reduction  prediction  for  bodies  of  revo- 
lution with  a polymer  ejection  process,  2)  to  qualitatively  define  the 
effects  of  polymer  on  boundary  layers  in  external  flows,  and  3)  to  veri- 
fy the  predictive  techniques  by  experiment.  Drop  tank  testing  was  per- 
formed with  several  especially  designed  bodies  of  revolution.  One  of 
these  bodies  was  capable  of  simultaneous  ejection  and  boundary  layer 
sampling  at  several  heights  in  the  boundary  and  at  several  axial  sta- 
tions. The  ejection  process  for  this  model  was  considered  optimal  In 
that  ejection  occurred  directly  forward  of  the  body  in  a laminar  flow 
region  thus  affecting  the  developing  flow  region  of  turbulent  flow. 

The  experimental  apparatus  is  capable  of  testing  a wide  variety  of  ejec- 
tors at  different  locations  and  over  a range  of  Reynold's  number  to 
about  5 x 10^.  Variable  ejection  rates  and  polymer  concentrations  may 
be  applied  as  well  as  boundary  layer  sampling  at  selected  heights  and 
at  the  wall.  The  drop  tank  facility  provided  a very  low  turbulence 
level  apparatus  in  which  the  external  flow  experiment  could  be  carried  out. 

2.  Analytical  Prediction  Methods 

Much  of  the  past  work  in  the  diffusion  of  polymer  boundary  layers 
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has  centered  in  the  far  downstream  zone  where  the  polymer  attained  the 
same  diffusion  characteristics  as  the  medium.  Postulates  of  an  inter- 
mediate zone  have  been  made  but  not  reasonably  tested  for  polymer  cases 
The  case  of  developing  flow  with  polymer  ejection  has  received  little 
attention  for  external  flows  and  especially  so  for  bodies  of  revolution 
Predictive  methods  are,  therefore,  not  well  developed. 


Several  analytical  modeis  were  developed  for  comparison  with  experi- 
mental data.  The  models  were  based  upon  a proven  integral  boundary 
layer  analysis  for  thick  axisyjmnetric  boundary  layers  combined  with  a 
velocity  profile  relation  accounting  for  pressure  gradients.  The  models 
made  use  of  a Lagrangian  similarity  hypothesis  which  predicts  a diffu- 
sion boundary  layer  growth  within  the  boundary  layer.  Two  versions 
of  the  model  are  proposed  but  only  one  extensively  used  in  the  com- 
parison. The  first  version,  combined  model,  presently  recommended  as 
capable  of  accurate  prediction  of  skin  friction,  polymer  wall  concen- 
tration and  indicating  total  drag  if  provided  with  an  experimental 
pressure  profile,  incorporates  an  initial  diffusion  zone  where  molecular 
diffusion  is  predomanate.  Based  upon  a nondimensional  distance,  given 
by 


(i)1'5  x_L  , 8 

V s ci 


(96) 


a change  to  an  intermediate  zone  diffusion  process  is  made.  Within 
this  region,  determined  by  experimental  data  for  ammonia  diffusion  in 
air  to  extend  to  X/6  < 60,  a similarity  concentration  profile  exists 
having  the  form 
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where  the  exponent  is  adjusted  to  a given  value  dependent  on  the  region 
of  diffusion.  Proceeding  to  X/iS  > 60  places  the  diffusion  process  in 
the  final  zone  where  the  second  value  of  is  applied.  The  Lagrangian 
analysis  provides  for  the  determination  of,  6^,  the  diffusion  boundary 
layer.  The  equation  for  the  diffusion  boundary  layer  development  incor- 
porates the  Karman  constant,  K,  which  is  adjusted  for  the  polymer  flow 
case  resulting  in  suppressed  diffusion. 

The  analytical  model  has  been  compared  to  experimental  data  and 
found  to  predict  proper  levels  of  wAll  and  concentration  profiles  but 
predicts  an  early  increase  in  the  diffusion  process.  It  is  hypothesized 
that  the  ejection  velocity  for  the  forward  ejector,  ejecting  into  the 
laminar  region,  is  not  the  valid  velocity  for  application  into  the  non- 
dimensional  distance  relation  (v^).  A mean  laminar  velocity  is  believed 
to  be  the  correct  velocity  for  this  ejection  process.  When  applied, 
the  predictive  techniques  properly  reproduces  the  experimental  data. 
Unfortunately,  only  limited  data  are  available  displaying  the  initia- 
tion of  the  intermediate  diffusion  zone. 

The  second  postulated  model,  the  variable  model,  eliminates 

the  need  for  the  molecular  diffusion  portion  of  the  previous  model  and 

adjusts,  K , as  a function  of  <$,  x/L,  and  c to  account  for  this 
j d w 

region.  Unfortunately,  insufficient  data  were  available  to  perfect  this 
model.  A simplified  functional  relationship  of  was  found  to  predict 
the  limited  experimental  data 
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The  combined  model  is  considered  valid  for  use  in  prediction  of 
various  flow  parameters  for  axisymmetric  bodies  of  revolution.  The 
model  should  be  exercised  in  several  ways  depending  on  the  ejection 
process.  If  ejection  occurs  prior  to  or  at  the  point  of  transition  to 
turbulent  flows,  the  full  model  should  be  applied.  If  the  ejector  is 
placed  in  a region  of  a developed  turbulent  boundary  layer,  the  initial 
zone  should  be  eliminated  and  the  Karman  constant  reduction  removed. 

3.  Boundary  Layer  Characterization  and  Separation  Tests 

Several  test  bodies  were  applied  in  this  series  of  experiments 
covering  a Reynolds  number  range  of  one  to  five  million.  One  series 
of  tests  were  conducted  in  a polymer  ocean,  where  the  drop  tank  was 
uniformly  mixed  to  a polymer  concentration  level  of  1.23,  2.5,  5,  10, 

20,  50,  and  60  WPPM  using  Polyox-WSR-301  polymer.  In  this  series,  models 
having  6°  tail  cone,  12°  tail  cone,  and  a spherical  tail  were  tested. 

The  photographic  study  displayed  several  interesting  results.  The  tur- 
bulent nature  of  the  boundary  layer  was  apparent  from  the  dye  aspirated 
from  these  bodies  at  the  minimum  pressure  point.  At  very  low  concen- 
tration levels,  2.5  WPPM,  the  fine  structure  turbulence  was  absent  and 
only  the  large  scale  turbulence  remained.  At  higher  concentration,  ex- 
treme boundary  layer  thinning  occurred.  Measurements  of  the  boundary 
layer  thickness  were  made  from  photographs  and  compared  to  computer 
predictions.  The  computer  routines  predicted  the  proper  boundary  layer 
growth  shape  for  these  tests  but  tended  to  overpredict  the  thickness. 

This  was  probably  due  to  the  mean  height  of  the  eddies  being  used  for 


measurement  purposes. 


No  shift  in  the  boundary  layer  separation  region  was  noted  by 
application  of  the  polymer.  This  was  undoubtedly  due  to  the  large  ad- 
verse pressure  gradient  conditions  for  these  tail  configurations. 

Computer  predictions  also  indicated  no  shift  in  separation  point  for 
the  test  conditions. 

The  testing  in  fresh  water  with  the  polymer  ejection  body  showed 
similar  results  but  provided  much  better  visual  representations  at  the 
higher  polymer  concentrations.  The  characteristic  coarse  and  fine  scale 
turbulence  was  noted  on  dyed  water  ejection.  Ejection  of  50  WPPM  dyed 
polymer  solution  eliminated  most  of  the  fine  structure  turbulence  leav- 
ing only  the  coarse  structure  and  somewhat  the  rapid  boundary  layer 
growth  in  the  tail  region.  Application  of  500  WPPM  and  1000  WPPM  of 
ejected  polymer  resulted  in  extremely  thin  boundary  layers,  even  in  the 
tail  region  with  little  evidence  of  turbulent  action.  In  the  tail  region, 
bursting  of  the  eddy  structure  is  evident  but  with  very  low  energy  levels 
as  evidenced  by  their  vectors  being  in  a near  axial  direction. 

At  the  higher  concentration  levels,  dye  streaking  is  apparent.  The 
streaking  spacing  is  consistent  with  laminar  sub  layer  streaks  but  also, 
for  this  case,  agrees  well  with  the  number  of  ejection  holes  in  the 
screen  ejector  in  the  nose  of  the  body.  In  any  event,  turbulent  mixing 
is  not  apparent  at  the  500  WPPM  and  1000  WPPM  concentrations.  Boundary 
layer  measurements  were  made  for  each  case  tested.  A normalization  of 
the  thickness  measurements  shows  an  increased  rate  of  growth  for  all 
concentrations  tested,  in  the  tail  region.  This  suggests  higher  diffu- 
sion rates  lowering  the  wall  concentrations  and  reducing  the  subdued 
growth  in  this  region. 
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4.  Drag  Reduction  Tests 

Drag  reduction  tests  were  performed  with  the  6°,  12°,  and  spherical 
tail  bodies  in  a polymer  ocean.  The  concentration  levels  were  as  dis- 
cussed previously.  A total  drag  reduction  of  33%  was  achieved  with  the 
6°  tail  body  at  polymer  concentrations  of  20  WPPM  and  greater.  This 
compared  favorably,  within  10%  of  calculated  values.  Percent  skin  fric- 
tion reduction  of  near  70%  were  achieved,  comparing  well  with  optimal 
total  shear  drag  reduction  for  flat  plates  as  determined  by  other  in- 
vestigators. The  maximum  total  drag  reduction  achieved  with  the  12° 
tail  and  spherical  tail  bodies  was  less  due  to  the  reduced  surface  area 
and,  therefore,  increased  percentage  of  form  drag.  Maximum  total  drag 
reductions  of  16%  and  10%  were  achieved  respectively  for  the  12°  tail 
and  spherical  tail  configuration.  The  skin  friction  reduction  remained 
near  70%  as  it  should.  Computer  predictions  for  total  drag  with  the 
12°  and  spherical  tail  models  was  considerably  in  error  due  to  the  very 
large  pressure  coefficients  predicted  near  the  separation  region  for 
these  cases  from  potential  flow. 

5.  Polymer  Diffusion  Tests 

The  polymer  ejecting  body  was  applied  for  all  tests  with  a fixed 

3 

ejection  rate  of  20.6  in  /sec  and  nominal  velocity  of  27  ft/sec.  Wall 
samples  and  boundary  layer  samples  at  .025  inches  and  .055  inches  were 
taken  and  analyzed  applying  a fluorometric  technique.  Testing  with 
water  ejection  and  5 WPPM,  10  WPPM,  20  WPPM,  50  WPPM,  500  WPPM,  and 
1000  WPPM  were  performed.  The  water  tests  compared  well  with  available 
predictive  methods  displaying  appropriate  concentration  profiles  as 
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achieved  hy  other  investigators.  The  tests  with  all  polymer  concentra- 
tions displayed  an  unexpected  drastically  reduced  diffusion  resulting 
In  a considerably  extended  initial  mixing  zone.  Application  of  a 
molecular  diffusion  model  in  this  region  agrees  well  with  data.  Boundary 
layer  concentration  profile  data  when  plotted  nondimensional ly  as 
y/6^  vs  c/cw  result  in  similarity  profile  exponent  coefficients,  K^, 
of  .75.  It  is  believed,  however,  the  value  is  below  this  since  data 
were  not  obtained  in  the  critical  region  of  approximately  .01  inches 
from  the  wall  in  these  tests.  These  results  are  significant  in  that 
they  imply  significantly  reduced  polymer  quantities  may  be  required 
with  proper  ejection  techniques.  Apparently  ejection  into  the  laminar, 
and  in  this  case  stagnation  region,  of  developing  flow  substantially 
reduces  the  turbulence  levels  associated  with  turbulent  flow  and  re- 
duces, therefore,  diffusion.  It  is  believed  that  ejection  in  a de- 
veloped flow  would  result  in  an  extremely  short  initial  zone  and  rapid 
mixing. 

6.  Conclusion  Summary 

An  experimental  apparatus  has  been  developed  which  allows  labora- 
tory measurement  of  polymer  drag  reduction  wall  concentration  profiles 
and  ejector  geometries  at  various  ejection  rates  on  axisymmetric  bodies. 
Experiments  have  displayed  reduced  fine  scale  turbulence  structure  with 
polymer  and  a significantly  reduced  diffusion  process  if  ejection  is 
performed  upstream  of  the  turbulent  transition  point.  Maximum  skin 
friction  reductions  were  achieved.  Predictive  methods  were  developed, 
building  on  available  theories,  for  prediction  of  the  wall  concentration 
on  an  axisymmetric  body.  Comparison  with  experiment  has  been  performed. 
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B.  RECOMMENDATIONS 


!.  The  variable  analytical  model  should  be  developed  to  more  ade- 
quately describe  the  diffusion  process.  As  presently  configured,  the 
model  abruptly  changes  from  molecular  diffusion  to  a diffusion  process 
which  is  believed  to  be  too  high  a rate  for  polymer.  Application  of 
the  Lagrangian  approach  with  modified  Karman  constant  should  provide 
an  input  to  the  functional  relationship. 

2.  Additional  data  should  be  obtained  at  boundary  layer  heights 
below  the  .025-inch  level  with  this  experimental  model.  These  data 
are  required  to  adequately  define  the  concentration  curves  and  provide 
the  data  necessary  for  developing  the  model.  The  variation  in  the 
similarity  profile  exponent,  K^,  with  axial  length  should  also  be  de- 
fined by  these  data  and  would  provide  information  on  the  intermediate 
and  final  zone  similarity  profile  exponents. 

3.  Additional  testing  should  be  performed  at  lower  concentration  levels, 
below  50  WPPM,  with  the  ejecting  experimental  model.  This  would  clearly 
establish  the  extent  of  the  initial  mixing  zone  and  verify  the  relation 
used  for  the  extent  of  the  initial  zone  and  the  hypothesis  for  the  ejec- 
tion velocity  equaling  a mean  laminar  velocity  for  the  type  of  ejection 
process  used  here. 

4.  The  experimental  model  should  be  modified  for  ejection  into  the 
established  turbulent  boundary  layer  region  to  clearly  demonstrate  the 
elimination  (or  drastic  reduction)  of  the  initial  zone  length  with  the 
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resulting  high  diffusion  process.  Significantly  higher  concentrations 
of  polymers  are  believed  to  be  required  (higher  flow  rates  and/or  higher 
concentrations)  for  this  case.  Model  modifications  are  extremely  simple. 

5.  Drag  reduction  measurements  should  be  made  with  this  model  during 
the  ejection  process.  Correlations  with  wall  concentration  measure- 
ments and  polymer  ocean  data  would  lend  further  credence  to  the  results. 
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APPENDIX  A 


SKIN  FRICTION  RELATIONS 

The  conservation  of  mass  and  the  momentum  equation  for  axisymmetric 
flow  are  given  by 

(pur)  + (pvr)  = 0 (Al) 

and 

/3u.  , /3u.  /Pe,  , 9 / \ /.on 

pur  (— ) + pvr  (— ) = - r (—)  + — (rr).  (A2) 

White's  derivation  is  reproduced  here  to  show  the  incorporation  of  the 
polymer  concentration  terms. 

Defining  a stream  function  which  satisfies  equation  (Al ) to  eli- 
minate v 
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Solving  (A3)  for  4*  yields 
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Rewriting  in  law  of  the  wall  variables 
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Writing  equation  (A2)  in  wall  variables  and  substituting  in  equation 
(A3)  to  eliminate  v 
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The  x derivatives  must  be  handled  by  the  chain  rule,  since  each  of  the 

parameters  (y"*",  a,  r , c ) in  the  law  of  the  wall  is  a function  of  x. 
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It  is  assumed  throughout  that  p and  v are  constant  in  this  analysis. 
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Terras  (a)  and  (b)  cancel  using  equation  (A6) . Multiplying  out,  cross 
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Where  the  integral  of  the  last  term  is  r x and  the  term 
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Dividing  by  p and  remembering  that  pv  = p yields 
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Defining  dimensionless  parameters 
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